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Introduction 


This document is the second volume in the three volume set of User’s Manuals for the Micromechanics 
Analysis Code with Generalized Method of Cells Version 4.0 (MAC/GMC 4.0). Volume 1 is the Theory 
Manual, this document is the Keywords Manual, and Volume 3 is the Example Problem Manual. 

MAC/GMC 4.0 is a computer code developed at NASA Glenn Research Center that analyzes the thermo- 
inelastic behavior of composite materials and laminates. The code is based on the micromechanics theory 
known as the generalized method of cells (GMC). GMC models the response of composite materials 
using a doubly or triply periodic repeating unit cell (RUC), which is composed of a number of subcells 
(see Figure 3.1 and Figure 3.2). By placing distinct constituent materials within the subcells, a 
heterogeneous (composite) material can be modeled, provided an RUC can be identified in the material’s 
micro scale architecture. The main advantages of GMC over other micromechanics theories and 
approaches are: 1) its fully multi-axial formulation, 2) the availability of local (constituent level) stress 
and strain fields, and 3) its computational efficiency. The availability of the local fields makes the 
method attractive in situations where more detailed composite analysis, beyond simple determination of 
effective thermo-elastic properties, is necessary. In such situations, when matrix inelasticity, constituent 
damage or failure, or fiber-matrix debonding is important, the availability of the local fields is critical in 
order to account for these effects appropriately. It is the purpose of the MAC/GMC 4.0 code to provide a 
convenient and user-friendly “package” for the GMC theories, in addition to providing significant added 
value through a library of inelastic constitutive models, repeating unit cells, thermo-mechanical and yield 
surface loading options, failure and damage analysis capabilities, and results generation options. 

The original method of cells (which allowed only four subcells) was developed by Aboudi (1989, 1991), 
the doubly periodic version of GMC by Paley and Aboudi (1992), and the triply periodic version of GMC 
by Aboudi (1995). The actual GMC theories implemented within MAC/GMC 4.0 were reformulated for 
maximum computational efficiency by Pindera and Bednarcyk (1999) and Bednarcyk and Pindera (2000). 
A recent advancement in the GMC technology has significantly improved the local accuracy of the 
doubly periodic version of the model. This new micromechanics theory, called the high-fidelity 
generalized method of cells (HFGMC), developed by Aboudi et al. (2001, 2002), has also been 
implemented with MAC/GMC 4.0. 

As mentioned, MAC/GMC 4.0 is also capable of analyzing composite laminates. The composite laminate 
analysis capabilities also rely on the GMC composite material model. The doubly and triply periodic 
versions of GMC function within the context of classical lamination theory (Jones, 1975; Herakovich, 
1998) to model the ply level composite material. Thus, the code can analyze the thermo-inelastic 
behavior of arbitrary laminate configurations. 

Finally, a new thermo-electro-magnetic version of GMC that is capable of modeling so-called “smart” 
composites has been implemented within MAC/GMC 4.0. The implementation also enables the analysis 
of smart composite laminates. This specific form of the thermo-electro-magnetic GMC was developed by 
Bednarcyk (2002) based on the work of Aboudi (2000). 


NAS A/TM— 2002-2 1 2077/VOF2 


2 



Introduction 


Current Capabilities and What's New in MAC/GMC 4.0 

The current capabilities of MAC/GMC are listed below, with new or improved features in this version 
highlighted in bold blue type. The superscripts refer to the enumerated details regarding these new 
features given after the list. 


Current Capabilities 

• Usability 

* Intuitive ASCII user interface 1 

* A graphical user interface (GUI) for MAC/GMC 4.0 will be available in March 2002 

* Computer resource management 2 

- Dynamic memory allocation 

- Static problem size limitations removed 

- Improved efficiency of solution for laminate analysis and large problems 

* 43 new example problems are described in the MAC/GMC 4.0 Example Problem Manual 

and are distributed with the code 

• Material Constitutive Models 5 

* Elastic : isotropic, transversely isotropic, completely anisotropic 

* Inelastic : Bodner-Partom, Robinson, isotropic GVIPS, transversely isotropic GVIPS, 

multimechanism GVIPS, SMA, incremental plasticity, Freed-Walker 

* User-defined : user-defined constitutive model subroutine ( usrmat . F90) 

• Library of Material Properties 

* Expanded internal material property library 

* Temperature-dependent or temperature-independent properties specified in input file 

* User-defined subroutine (usrfun.F90) for calculating material properties as a function of 

temperature or other variable 

* External material property database file 6 

• Micromechanics Analysis Models 

■+ Doubly periodic generalized method of cells 

* Triply periodic generalized method of cells 

* Doubly Periodic High-Fidelity Generalized Method of Cells (HFGMC) 4 

• Library of Doubly and Triply Periodic Repeating Unit Cell Architectures 

• Laminate Analysis Model (Classical Lamination Theory) 

* Monolithic layers 

* Doubly periodic composite layers 

* Triply periodic composite layers 7 

* Improved interoperability (doubly and triply periodic layers in the same laminate) 7 

• Smart Material Analysis Capabilities 

Shape memory alloy (SMA) model 5 

* Fully coupled thermo-electro-magneto-elasto-plastic analysis 3 

- Constituent material (arbitrary poling direction) 

- Composite 

- Laminate 
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• Library of Applied Loading History Options 

■ Thermal, uniaxial, biaxial, general (arbitrary) loading, thermo-mechanical 

* Curvature and moment loading for laminate 

* Electromagnetic loading 

• Time Integration Options 

* Forward Euler 

* Self-Adaptive Predictor-Corrector 

• Yield Surface Analysis 8 

* GMC, HFGMC, laminate 

' Global (composite/laminate level), local (subcell level), ply level in laminate 

• Failure, Damage, and Lifing Capabilities 

* Elastic Allowables Estimation 9 

* Static Failure Analysis 9 

- Doubly periodic GMC, triply periodic GMC, laminate 

- Subcell level, repeating unit cell level 

- Maximum stress, maximum strain, and Tsai-Hill criteria 

* Weak Interface Modeling 

- Distinct interfacial phase layer 

- Interfacial displacement discontinuity modeling (GMC and HFGMC) 

* Fiber Breakage Modeling 11 

- Using internal fiber interface 

- Curtin effective fiber breakage model 
Fatigue Damage Analysis 10 

- Doubly periodic GMC, triply periodic GMC, laminate 

- Stiffness reduction, strength reduction 

- Improved interoperability - can be combined with static failure, debonding, and 
fiber breakage 

• Output and Data Visualization 12 

* Unlimited number of macro, micro, and laminate x-y plot files 

* MSC/PATRAN output for use with the MACPOST software add-on fringe plots 

MATLAB output for use with distributed MATLAB source files fringe plots 


Details on New Features 

1) The format of the MAC/GMC 4.0 input file has been completely overhauled. Unneeded keywords 
have been eliminated, more logical keyword and specifier names have been utilized, and the 
placement of specific data has been better organized. The result is a considerably more streamlined 
and easy to use input file format. 

2) The memory allocation of MAC/GMC 4.0 is dynamic. There is no longer a static limit on the size of 
the problem (i.e., number of subcells and laminate plies) that can be analyzed, the number of 
constituent materials that can be employed, or the number of x-y plot result files that can be 
generated. The result is a more robust code that only uses the system resources it needs to execute a 
given problem. 


NAS A/TM— 2002-2 1 2077/VOL2 


4 


Introduction 


3) MAC/GMC 4.0 can analyze smart composites and laminates. A major advancement in the 
generalized method of cells micromechanics models and lamination theory has been implemented 
within MAC/GMC 4.0: the ability to model fully coupled thermo-electro-magneto-visco-elasto- 
plastic composite materials and laminates. These so-called “smart” materials and composites react 
electromagnetically to thermo-mechanical loading and react mechanically to electromagnetic loading. 
These materials are currently finding a growing number of applications in sensors and actuators and 
have the potential to enable the next generation of multi-functional structures that can automatically 
sense and react to various stimuli. See problems 7b and 7c in the MAC/GMC 4.0 Example Manual. 

4) A new high-fidelity doubly periodic micromechanics model has been implemented within 
MAC/GMC 4.0. This new micromechanics model, the high-fidelity generalized method of cells 
(HFGMC) provides much more accurate local stress and strain field predictions than does the 
standard method of cells. MAC/GMC 4.0 can therefore perform composite analyses with a variable 
level of fidelity, maximizing either accuracy or efficiency when appropriate. See problem 3f in the 
MAC/GMC 4.0 Example Manual. 

5) Several new internal material constitutive models have been implemented within MAC/GMC 4.0. 
These include: a shape memory alloy constitutive model, the classical incremental plasticity model, a 
new multimechanism visco-elastic-plastic GVIPS model, the Freed-Walker viscoplasticity model, a 
Robinson isotropic viscoplasticity model for NARloy Z, and a completely anisotropic elastic model. 
The classical incremental plasticity model is formulated to allow the user to simply input stress-total 
strain point pairs from a material’s inelastic stress-strain curve. See problems 2c, 2d, 7a, and 7e in 
the MAC/GMC 4.0 Example Manual. 

6) MAC/GMC 4.0 enables the use of an external material database file. User’s can store their own 
material properties in this external database file, from which the code reads and stores only the 
material properties specified for use in a specific problem. See problem 2h in the MAC/GMC 4.0 
Example Manual. 

7) The laminate analysis capabilities have been significantly expanded and computationally enhanced. 
Discontinuous fiber and particulate reinforced laminate layers can now be analyzed thanks to the 
ability of triply periodic GMC to represent any particular layer within a laminate. Further, all 
laminate data is now kept in core memory during execution (previous versions of MAC/GMC 
performed read/write cycles to disk). This has significantly increased the computational efficiency of 
MAC/GMC 4.0 laminate analyses. See problem 3i in the MAC/GMC 4.0 Example Manual. 

8) The yield surface analysis capabilities of MAC/GMC 4.0 are significantly expanded. Various subcell 
and repeating unit cell level yield surfaces can be automatically generated for composites using 
standard or high-fidelity GMC. Further, yield surfaces can now be automatically generated for 
laminates on the laminate, ply, and micro (subcells within a ply) levels. See problems 4f and 4g the 
MAC/GMC 4.0 Example Manual. 

9) The static failure analysis capabilities of MAC/GMC 4.0 have been significantly enhanced. A simple 
and quick estimation of composite allowable stress and strain quantities can be performed based on 
the constituent allowables. Further, true subcell and repeating unit cell level static failure analysis 
(i.e., checking for failure at various scales during application of simulated thermo-mechanical 
loading) can be performed using several failure criteria. These failure analyses are available for 
composites as well as laminates. See problems 5a, 5b, and 5c in the MAC/GMC 4.0 Example 
Manual. 
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10) The fatigue damage analysis capabilities of MAC/GMC 4.0 have been significantly enhanced. A new 
strength reduction fatigue model has been implemented, and the applicability of the stiffness 
reduction fatigue model has been expanded to include triply periodic repeating unit cells and 
laminates. Further, both fatigue damage models now work in concert with the MAC/GMC 4.0 static 
failure and interfacial debonding capabilities to enable mechanistically realistic life predictions for 
composites and laminates. See problem 5d in the MAC/GMC 4.0 Example Manual. 

11) The local fiber failure analysis capabilities of MAC/GMC 4.0 have been improved via inclusion of 
the Curtin effective fiber breakage model. Based on fiber strength statistical data, the Curtin model is 
a widely used tool for predicting composite stiffness degradation and failure due to longitudinal fiber 
fracture. The Curtin model can be used to model fibers within continuous fiber composites and 
continuous fiber composite layers within laminates. See problem 5f in the MAC/GMC 4.0 Example 
Manual. 

12) The post-processing of MAC/GMC 4.0 results has been significantly improved. The code now 
automatically generates output for use with the popular MATLAB software package. MATLAB 
source files are distributed with MAC/GMC 4.0 that generate fringe plots of the stress, strain, and 
inelastic strain fields predicted by the code. See problem 6c in the MAC/GMC 4.0 Example Manual. 

13) MAC/GMC 4.0 has been extensively tested and debugged. 

^ Note: Attempts have been made to enable a high level of interoperability within MAC/GMC 4.0. 

However, currently not all of the code’s features can work in concert. For example, the 
composite layers within a laminate cannot currently be modeled using the new HFGMC 
micromechanics model. Look for additional improvements in interoperability in future 
versions of MAC/GMC. 

Format of this Manual 


This Keywords Manual first presents a “Getting Started” section that details the installation of the 
MAC/GMC 4.0 software and how the code can be executed. The bulk of this manual presents detailed 
information on each of the keywords that may appear in the MAC/GMC 4.0 input file. Each keyword is 
presented with its purpose, and then the usage of and all specifiers associated with the keyword are 
described. Optional and special case specifiers are shown in blue type. A “Notes” section for each 
keyword is referenced via superscripts in the list of specifiers and their descriptions. For example, the 
specifier t(i) is described as, 

Name Definition 

t(i) 3 Time points in time vs. load magnitude profile 

Here, the superscript “3” would refer to note number 3 in the “Notes” section, which would follow the 
above description. 

One or more examples of the usage of each keyword are given along with any additional information 
needed to use the keyword. Any example problems (appearing in the Example Problem Manual) 
associated with the particular keyword are listed, with example problems that are most relevant to the 
keyword typed in bold. Suggested sources of additional information are also provided when appropriate. 
In some cases, keyword information is broken into subsections in order to facilitate the presentation of the 
needed information. Note that the abbreviation N/A is employed throughout this manual to indicate that a 
particular option is not applicable or unavailable. Finally, - AND - and - OR - are used as logical 
operators. 
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Getting Started 


Distributed Files 


The MAC/GMC 4.0 distribution cd includes the following in its root directory: 


Of mac4 . exe MAC/GMC 4.0 executable file 


□ library 
ra examples 
J23 manuals 

□ documents 
jfif * macpost 

□ matlab 


File folder containing the MAC/GMC 4.0 library file and user-defined 
subroutines 

File folder containing the MAC/GMC 4.0 Example Problem input files and the 
sample external material database file 

File folder containing the MAC/GMC 4.0 Theory Manual, Keywords 
Manual, and Example Problem Manual in PDF format 
File folder containing many of the technical references cited in the 
MAC/GMC 4.0 manuals in PDF format 

File folder containing MACPOST files required for MAC/GMC 4.0 results 
data post-processing using the MSC/PATRAN software package 
File folder containing the MATLAB source files required for MAC/GMC 
4.0 results data post-processing using the MATLAB software package 


It is recommended that the user create the directory “mac 4” on their personal computer running the 
Microsoft Windows® operating system and copy the entire contents of the cd to this directory. If disk 
drive space is at a premium, the contents of the “manuals” and “documents” folders can be omitted and 
the user can employ the cd when viewing the documentation contained in these folders. Note that, in 
order to read the documentation files provided in the “manuals” and “documents” folders, the user will 
need the Adobe Acrobat® Reader software installed on his/her computer. This software can be 
downloaded for free from: 


http ://www . adobe . com/products/acrobat/readstep .html 


Executing MAC/GMC 4.0 Problems 

The MAC/GMC 4.0 executable file (mac4.exe) is a DOS executable file intended to be executed 
through Microsoft Windows® (NT, 2000, XP) at the command prompt. In order to execute a MAC/GMC 
4.0 problem, at the command prompt, in the directory containing the executable, type: 

mac 4 in_name out_name 

where, 

in_name MAC/GMC 4.0 input file name 
outjiame MAC/GMC 4.0 output file name (Optional) 
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The code will look for the input file with the specified injiame and a “.mac” extension, so the specified 
injame should not include an extension. The outjame specification is optional. If outjame is 
omitted, the code will use the specified input file name and write output to the file injiame. out. If 
outjame is specified and does not have an extension, the code will write output to the file 
out jame . out. Finally, if outjame is specified and does have an extension, the code will write output 
to the file outjame with the specified extension. If the desired input and output file locations are 
different than the current working directory, the directory path can be added to the specified injame and 
outjame. 

Example : Execute MAC/GMC 4.0 using the input file “example_la.mac” located in the “mac4” 

directory, write output to the file “example_la . out” located in the “mac4” directory. 

C:\mac4>mac4 example_la 

Example : Execute MAC/GMC 4.0 using the input file “example_la.mac” located in the “mac4” 

directory, write output to the file “test . out” located in the “mac4” directory. 

C:\mac4>mac4 example_la test 


In the two example cases above, the following is written to the screen: 


+ + 

+ Micromechanics Analysis Code + 
+ with + 
+ Generalized Method of Cells + 
+ + 
+ MAC/GMC + 
+ + 
+ Version 4.0 + 
+ + 
+ *** WARNING *** + 
+ No warranties are granted regarding the results of this code + 
+ This code is export controlled by the U.S. government + 
+ + 


**** Checking Input File **** 

> Input file name: c:\mac4\example_la.mac 

> 0 Error (s) Detected < 


Determining Effective Properties 


* * * 


*************************** 

*** EXECUTION COMPLETED *** 

*************************** 

Execution Time: 0.016 second(s) 

As the above text indicates, MAC/GMC 4.0 performs a self-diagnostic “prescan” step prior to executing a 
given problem. During this prescan, the input file is checked for errors. If errors are detected, execution 
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of the problem will not be completed, and a description of the detected errors will be written to the output 
file. 

The MAC/GMC 4.0 input, output, and results data files can be viewed and edited with any standard 
ASCII file editor (such as Microsoft NotePad or Microsoft WordPad). Alternatively, MAC/GMC 4.0 can 
be set up to execute within a programming environment, such as Microsoft Visual Workspace. The best 
way to do this is to compile and link the MAC/GMC 4.0 library file and the user-defined subroutines 
within the programming environment, as described in the next section. The MAC/GMC 4.0 files can then 
be conveniently manipulated within the utilized programming environment. 


Using User-Defined Subroutines with MAC/GMC 4.0 

The MAC/GMC 4.0 static library file, along with the three MAC/GMC 4.0 user-defined subroutines, are 
included on the MAC/GMC 4.0 distribution cd in the “library” file folder. The contents of the 
“library” folder are: 


mac4 . lib MAC/GMC 4.0 static library file 

usrf ormde . F90 MAC/GMC 4.0 user-defined FORTRAN 90 source file for determining 
constituent material elastic stiffness matrices from user-defined elastic 
constants (called when usrmat . F90 is employed) 
usrf un . F90 MAC/GMC 4.0 user-defined FORTRAN 90 source file for determining 

constituent material properties from a user-defined function of temperature 
or other variable 

usrmat . F90 MAC/GMC 4.0 user-defined FORTRAN 90 source file for determining 
inelastic strain and state variable increments from a user-defined inelastic 
constitutive model (material properties can be taken from the input file, the 
external material database, or the usrfun . F90 routine) 

In order to generate a MAC/GMC 4.0 executable file from these four files, the . F90 source files must be 
compiled and linked with the mac4 . lib static library file using an appropriate Microsoft Windows® (NT, 
2000, XP) FORTRAN 90 compiler. MAC/GMC 4.0 was developed using the Compaq Visual Fortran 
product (v. 6.6), and it is thus recommended that this product be used. For information on this product, 
see: http://www.compaq.com/fortran/ 

The procedure for creating the MAC/GMC 4.0 executable from the four files above within Compaq 
Visual Fortran involves: 

1) Create new fortran console application project: 

a) Click “File” — » “New” — » “Projects” Tab — » select “Fortran Console Application” 

— » specify Location (directory path) and Project name and click “OK” 

b) Select “An empty project” to create, click “Finish” and click “OK” 

2) Add files to the project: 

a) Click “Project” -> “Add to Project” -> “Files” 

b) Select the MAC/GMC 4.0 static library file and the three user defined source files 

listed above from the appropriate directory and click “OK” 
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3) Choose “Release” configuration to maximize execution speed 

a) Click “Build” — » “Set Active Configuration ...” 

b) Select “. . . Win32 Release” and click “OK” 

4) Build (compile and link) the executable 

a) Click “Build” — » “Build projectjiame ” 

b) Alternatively, just hit F7 on the keyboard 

This procedure will build an executable MAC/GMC 4.0 file in the “Release” directory. If a debug 
version is desired (which will allow the debugging of the user-defined subroutines only), the active 
configuration can be set to “Debug”. 

To execute MAC/GMC 4.0 problems using this executable file from within the Compaq Visual Fortran 
environment: 

1) Click “Project” — » “Settings” 

2) Make sure that on the left the “Settings For:” box has the appropriate configuration listed and that the 
project name (rather than one of the source files or the library file) is highlighted. 

3) On the right, choose the “Debug” tab 

4) Under “Program arguments:”, type: in_name outjiame (as you would after typing “mac 4” at for 

execution via the DOS command prompt) and click “OK” 

5) Select “Build” — » “Execute ...” to run the code. (Alternatively hit the ! toolbar button) 

Note that the MAC/GMC 4.0 input, output, and data results files can be viewed and edited within the 
Compaq Visual Fortran environment. The user-defined FORTRAN 90 source files can also be edited 
within this environment. It is important for the user to realize that the three provided user-defined 
FORTRAN 90 source files contain code that is used to demonstrate the MAC/GMC 4.0 user-defined 
capabilities. These same three subroutines are complied and linked into the MAC/GMC 4.0 distribution 
executable file (mac4 . exe). If users wish to employ their own code in the user-defined subroutines, this 
code must be placed in the three user-defined FORTRAN 90 source files, and these files must be 
compiled and linked to the MAC/GMC 4.0 static library file (see procedure above). It is the user’s 
responsibility to ensure that his code is functioning properly and that his code is actually being utilized 
during execution rather than the demonstration code distributed in the three user-defined FORTRAN 90 
source files. Note that a message is written to the screen when the distributed (default) versions of the 
user-defined subroutines are called. It is recommended that the user write a similar message to the screen 
to serve as notification that the intended subroutine is being called. 


Format of the MAC/GMC 4.0 Input File 

The remainder of this Keywords Manual deals with the specific input requirements of the MAC/GMC 4.0 
input file to execute analyses of composite materials and laminates. This input file is an ASCII file that 
can be edited with any standard ASCII file editor. The first line of the input file is always interpreted as a 
title line. The remainder of the input file is divided into a number of keywords, which are placed on their 
own line in the file and are preceded by the asterisk character (*). A number of input lines follows each 
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keyword, and each line contains specifier words, followed by an “=” sign. The input data read by the 
code then follows the “=” sign. An example of a MAC/GMC 4.0 input file is: 

MAC/GMC 4.0 Example la - graphite/epoxy effective properties 

* CONDUCTIVITY 

NTEMP=1 TEMP=2 1 . 

^^CONSTITUENTS 

NMATS=2 

# -- Graphite fiber 

M=1 CMOD=6 MATID=U MATDB=1 & 

EL=388 . 2E9,7.6E9,0.41,0.45,14.9E9,-0. 68E-6, 9 . 74E-6 
K=500. , 10. 

# -- Epoxy matrix 

M=2 CMOD=6 MATID=U MATDB=1 & 

EL=3 . 45E9, 3 . 45E9, 0.35, 0.35, 1 . 278E9, 45 . E-6 , 45 . E-6 
K=0 .19,0.19 
*RUC 

MOD=2 ARCHID=1 VF=0 . 65 F=1 M=2 

* PRINT 

NPL=-1 

*END 

This MAC/GMC 4.0 input file contains five keywords (highlighted in bold). In addition, the MAC/GMC 
4.0 input file recognizes the following two special characters: 

# Placed at the start of a line, makes the line a comment (which is ignored by the code) 

& Placed at the end of a line, serves as a continuation character and causes the code to treat 

the next line as part of the same line. 

These two special characters are very important as the the line on which the data specifiers are located 
within the MAC/GMC 4.0 input file is important. That is, examining a portion of the example input file 
given above, 

# -- Graphite fiber 

M=1 CMOD=6 MATID=U MATDB=1 & 

EL=388.2E9,7.6E9,0.41,0.45,14.9E9,-0.68E-6, 9.74E-6 
K=500 . , 10. 


The specifiers matid=, matdb=, and el= are all on the same line (due to the presence of the 66 &” 
continuation character), but the specifier K= is on its own line. If the correct input file format is not 
followed, the code will not execute and the errors detected within the input file will be described in the 
output file. 

Finally, the order of the keywords in the input file is unimportant. The order of the specifiers on a 
particular line is also unimportant. As mentioned above, however, the order of the specifiers that are on 
different lines under a particular keyword is important. Despite the freedom to place the keywords and 
specifiers on a particular line in any order, this manual presents a logical suggested order for both. It is 
recommended that the user follow this order for the sake of logic and clarity. 
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Section 1 : Flag-Type Keywords 


The first section of the MAC/GMC 4.0 keywords involves the flag-type keywords. These three keywords 
are like switches that tell the code whether or not to perform certain analyses or do certain things. Like 
all keywords, they may be placed anywhere within the MAC/GMC 4.0 input file, but it is recommended 
that they be placed at the beginning of the file. 

The three MAC/GMC input file keywords covered in this section are: 

* CHECK — » Causes code only to check the input file without executing. 

* CONDUCTIVITY — » Causes code to calculate effective thermal conductivities. 

* ELECTROMAG —> Indicates that code execution will involve electromagnetic (smart) analysis. 


It should be noted that if * CONDUCTIVITY is specified, there are additional input requirements for and 
restrictions on the constituent materials (specified with CONSTITUENTS). Similarly, if 

* ELECTROMAG is specified, additional input requirements are associated with the constituent materials 
(specified with CONSTITUENTS), the repeating unit cell or laminate (specified with *RUC or 

* LAMINATE), and the applied mechanical loading (specified with *MECH). 
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* check: Check input file for errors 


* CHECK: Check input file for errors 


Purpose: 

Specify that MAC/GMC 4.0 should only check the input file for errors and then stop. The code will not 
execute an analysis or generate any results. 


Usage: 

♦CHECK 


Example: 

♦CHECK 


Associated MAC/GMC 4.0 Example Problemfsn 

None 


Notes: 

1) MAC/GMC 4.0 always performs a “prescan” of the input file in order to check the data for errors and 
determine array dimensions for dynamic memory allocation. The * CHECK option stops the code 
after the prescan is complete. 

2) There are some input data inconsistencies that cannot be detected during the prescan. Such errors will 
not be detected using the * CHECK option, but will be detected during analysis execution, causing the 
code to stop. 
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* CONDUCT I VI T Y : Effective thermal conductivities 


Purpose: 

Specify that MAC/GMC 4.0 should determine effective thermal conductivity components at the specified 
temperatures. 


Usage: 

* CONDUCTIVITY 

NTEMP =ntemp TFMP=T(l),T(2),...,T(ntemp) 

Name Definition 

ntemp Number of temperatures at which to calculate effective thermal conductivities 
T(i) Temperatures at which to calculate effective thermal conductivities 


Example: 

♦CONDUCTIVITY 

NTEMP=3 TEMP=21 . , 200 . , 500 . 

Associated MAC/GMC 4.0 Example Problemfs): 

Example la 


Notes: 

1) MAC/GMC 4.0 determines the three effective thermal conductivity components (in the three 
Cartesian coordinate directions) of a GMC repeating unit cell (RUC) at the specified temperatures as 
a separate pre-step. This is done prior to determination of the effective thermo-elastic properties and 
application of the simulated loading. The effective thermal conductivities are printed to the output 
file and then no longer stored or used by the code. 

2) If the * CONDUCTIVITY keyword is specified, two thermal conductivity components (axial and 
transverse) must be specified for each constituent material read by MAC/GMC 4.0 (see Section 2, 
* CONSTITUENTS). Further, the material properties of the constituents must be user-defined 

(MATID=U). 

3) The * CONDUCT I VI T Y keyword is only compatible with repeating unit cell analysis (*RUC), not 
laminate analysis (*LAMINATE). 
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Section 1: Flag-Type Keywords 


* conduct I vi T Y : Effective thermal conductivities 


4) The effective thermal conductivities are calculated based on the assumption that all constituent 
materials are (at most) transversely isotropic (in terms of their local thermal conductivities) and that 
their plane of isotropy is the x 2 -x 3 plane in global RUC coordinates. 


Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 2 (*CONSTITUENTS) 
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* ELECTROMAG : Electromagnetic analysis 


Purpose: 

Specify that MAC/GMC 4.0 should perform coupled thermo-electro-magneto-elasto-visco-plastic 
analysis. 


Usage: 

*ELECTROMAG 


Example: 

*ELECTROMAG 


Associated MAC/GMC 4.0 Example ProblemfsU 

Example 7b, Example 7c 


Notes: 

1) Unlike determination of the effective thermal conductivities, electromagnetic analysis is tightly 
coupled within the MAC/GMC 4.0 formulation. Field variables associated with the thermo-electro- 
magnetic effects are stored and tracked during the MAC/GMC 4.0 execution and simulated electro- 
magnetic loading must be applied. 

2) If the * ELECTROMAG keyword is specified, each constituent material read by MAC/GMC 4.0 must 

include the electromagnetic specifier: EM=1 indicates that the constituent material has 

electromagnetic material parameters associated with it; EM=0 indicates that the constituent material 
does not have electromagnetic material parameters associated with it (see Section 2, 
* CONSTITUENTS). For constituent materials with EM=1, electromagnetic material parameters must 
be specified, the material’s thermo-mechanical material properties must be user-defined (MATID=U), 
and the arbitrary transversely isotropic elastic constitutive model must be employed (CMOD=9). 

3) Electromagnetic analysis may be performed for both repeating unit cell analysis (*RUC) and laminate 
analysis (* LAMINATE). However, the electromagnetic materials may only be part of triply periodic 
GMC RUCs. That is, if RUC analysis is specified, triply periodic GMC must be employed (MOD=3). 
If laminate analysis is specified, any layer that is electromagnetic (i.e., contains electromagnetic 
constituent materials) must employ triply periodic GMC (MOD=3). 
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Section 1: Flag-Type Keywords 


*electromag: Electromagnetic analysis 


4) In the case of electromagnetic laminate analysis (* LAMINATE), each layer must include the 
electromagnetic specifier: EM=1 indicates that the layer is electromagnetic (see Section 3, 

* LAMINATE). 

5) Because electric and magnetic loading terms must be specified in the case of electromagnetic 
analysis, the arbitrary loading option (LOP=9 9) must be employed. The electromagnetic loading 
history is specified as mechanical loading components 7-12 (see Section 4, *MECH). 


Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 2 (^CONSTITUENTS) 
MAC/GMC 4.0 Keywords Manual Section 3 (*LAMINATE) 
MAC/GMC 4.0 Keywords Manual Section 4 (*MECH) 

MAC/GMC 4.0 Theory Manual Section 2.2.2 
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Section 2 : Material Keywords 


This section deals with two keywords related to the constituent materials that will form the repeating unit 
cell or laminate to be analyzed by MAC/GMC 4.0. These keywords are used to provide the code with the 
constituent material properties that will occupy the subcells in the analyzed repeating unit cell (RUC) or 
repeating unit cells that represent each ply of the laminate. Then, using keywords described in Section 3, 
the RUC is constructed using subcells, which contain some or all of the constituent materials. If the 
MAC/GMC 4.0 problem involves a laminate, the laminate is constructed of layers, each of which is 
represented by an RUC. The RUC for each layer is then constructed of subcells containing the 
constituent materials. Because the constituent materials always occupy subcells within an RUC, they can 
also be thought of as subcell materials. However, the subcell materials are not necessarily elastic or 
isotropic. MAC/GMC 4.0 handles a wide range of material visco-elasto-plastic material behavior and 
material anisotropy on the level of the subcell. 

The two MAC/GMC 4.0 input file keywords covered in this section are: 

★MDBPATH — » Specifies name and location of the external material database file. 

* CON S T I TUENT S — » Specifies the properties of constituent materials that will occupy subcells. 

As one might expect, there are several ways in which MAC/GMC 4.0 can obtain the material properties 
for the constituent materials. These include: 

1) Reading the properties directly from the MAC/GMC 4.0 input file. 

2) Obtaining the properties from the internal MAC/GMC 4.0 material database. 

3) Reading the properties from a user’s personal external material database file. 

4) Obtaining the properties from a user-defined subroutine that calculates them based on a specified 
function. 

In addition, different material properties are required for different analysis types and constitutive models. 
For example, if the user indicates that effective thermal conductivities should be calculated by the code 
(by including *CONDUCTIVITY in the input file), the thermal conductivities of the constituent materials 
must be specified in *CONSTITUENTS. Similarly, if the user chooses an inelastic constitutive model for 
a particular constituent material rather than a linear elastic constitutive model, material parameters for the 
particular inelastic constitutive model must be specified in addition to the material elastic constants. 
Because of this specialization of the constitutive material properties for many of the different problems 
that MAC/GMC 4.0 can solve, the data specification associated with the CONSTITUENTS keyword is 
also quite specific compared to the other input parameters associated with the problem. As such, the part 
of this section dealing with the CONSTITUENTS keyword is the longest and most involved component 
of this manual. 

It is important to realize that the employed constituent material property units must be consistent. That is, 
all constituent materials that comprise an RUC or laminate must employ the same units. The units of the 
material properties within the MAC/GMC 4.0 internal material database are given in Table 2.2. Both ksi 
and MPa are used for the material property stress (and stiffness) unit. The time unit for all material 
properties is seconds, the temperature unit is °C, and the strain values are unitless (as opposed to %). 
Note also that any applied loading should be specified in units consistent with the material property units. 
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*mdbpath : External material database file & location 


*MDBPATH : External material database file & location 


Purpose: 

Specify the name and path to an external material database file. 


Usage: 

*MDBPATH 

PATH =mat _path NAME =mat_name 

Name Definition 

mat jpath 1 Computer directory path to the external material database file (Optional) 
mat name Name of the external material database file with extension 


Examples: 

*MDBPATH 

NAME=sample_material_database . mat 

*MDBPATH 

PATH=c : \mac4 \mat_db\ NAME=sample_material_database . mat 

Associated MAC/GMC 4.0 Example Problemfsb 

Example 2h 


Notes: 

1) If the optional PATE=mat _path is omitted, the code will look for the external material database file in 
the same directory as the input file. 

2) A sample external material database file, sample_material_database .mat, is distributed with 
MAC/GMC 4.0. This sample file is used in Example Problem 2h. A print out of this sample file is 
given in the Appendix of this manual. 

3) The format of the external material database file is similar to the user-defined material property input 
format employed when the properties are specified directly in the input file. The specifics of the 
external material database file format are discussed under the CONSTITUENTS keyword in the 
present section of this manual. 

4) To choose a particular material from the external material database file, specific input is required 
under CONSTITUENTS in order to specify the name of the material from the external material 
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database file and that the code should look for the particular material in the external database file (see 
* CON S T I TUENT S in the present section of this manual). 

5) MAC/GMC 4.0 does not read in and store all constituent materials in the external material database 
file. The code only reads and stores the materials specified under *CONSTITUENTS. Thus, 
requirements that apply to all constituent materials only apply to those that are actually read from the 
external material database file. For example, electromagnetic analysis requires that all constituent 
materials have the electromagnetic specifier (em=) included. If using an external material database, 
this requirement would not apply to every material appearing in the external material database file. It 
would only apply to those read by MAC/GMC 4.0 as specified under *CONSTITUENTS in the 
particular input file being employed. 


Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 2 ^CONSTITUENTS) 
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Section 2: Material Keywords 


♦constituents : Constituent (subcell) materials 


* CONSTITUENTS : Constituent (subcell) materials 


^ * CONSTITUENTS is a required keyword ^ 

Purpose: 

Specify the constituent material constitutive models and material properties that will be read and stored 
by MAC/GMC 4.0. The local constitutive models supported by MAC/GMC 4.0 are listed in Table 2.1. 
Each is described in detail in the MAC/GMC 4.0 Theory Manual Section 4. 


Table 2.1 Constituent material constitutive models supported by MAC/GMC 4.0. 


cmod 

Model 

i 

Bodner-Partom viscoplastic model 

2 

Modified (non-isothermal) Bodner-Partom viscoplastic model 

3 

Robinson transversely isotropic viscoplastic model 

4 

Isotropic GVIPS model 

6 

Standard elastic model (transversely isotropic with RUC x 2 -x 3 plane of isotropy) 

7 

Transversely isotropic GVIPS model 

9 

Transversely isotropic elastic model (arbitrary specified plane of isotropy) 

10 

Freed-Walker viscoplastic model 

11 

Robinson isotropic viscoplastic model for NARloy Z 

15 

Anisotropic elastic model 

21 

Classical incremental (rate-independent) plasticity model 

22 + 

Multimechanism visco-elasto-plastic GVIPS model 

30 

Graesser - Cozzarelli - Witting shape memory alloy (SMA) model 

99 

User-defined constitutive model (requires compilation of user code in usrmat.F90 and 
usrformde . F90 and linking with the MAC/GMC 4.0 library file mac4 . lib) 


+ Non-U. S. government users may require an additional software license to use this constitutive model 


The material information specified under * CONSTITUENTS has many special cases based on the options 
chosen by the user and can become somewhat complicated. As such, the description for this keyword has 
been written in outline form in order to present the input requirements as simply as possible. Each outline 
heading has its own Usage description and Examples. The Notes have been consolidated and placed at 
the end of each capital lettered section. The outline format for the CONSTITUENTS keyword is: 


I. Properties taken from the internal material properties database ( matid ^ u) 

II. User-defined material properties ( matid = u) 

A. Material properties read directly from the input file ( matdb =1) 

1. Standard constitutive models (cmod^ 15, cmod^ 9 9) 

a) Temperature-independent material properties 

b) Temperature-dependent material properties 

2. Special cases 

a) Internal anisotropic constitutive model (cmod= 15) 

i) Temperature-independent material properties 

ii) Temperature-dependent material properties 
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b) User-defined constitutive model (cmod= 9 9) 

i) Temperature-independent material properties 

ii) Temperature-dependent material properties 

B. Material properties determined from a user-defined function ( matdb = 2) 

C. Material properties read from external material database file ( matdb = 3) 

III. Electromagnetic (smart) material properties 

A. Temperature-independent material properties 

B. Temperature-dependent material properties 


I. Internal material properties database (matid ^ U) 


Usage: 

^^CONSTITUENTS 

NMATS =nmats 

n=matnum cmod =cmod tref =tref matid =matid D=d 1 ,d 2 ,d 3 
M=matnum CMOD =cmod TREF =tref MATID =matid D=d 1 ,d 2 ,d 3 


—> Repeated nmats times 


Name 

nmats 

matnum 

cmod 

tref 1 ' 2 

matid 


d i, d 2 , d 3 


3 


Definition 

Number of constituent materials specified 

Material number (data following M =matnum is associated with material number matnum) 
numbered sequentially 

Constitutive model identification number (see Table 2.1) 

Reference temperature for material (Special Case) 

Material i.d. letter, which indicates location of material parameters 
= a-h — » Internal material database (see Table 2.2) 

= U —> User-defined material properties 

Components of vector normal to the plane of transverse isotropy (Special Case) (see 
Figure 2.1) 


Example: 

♦CONSTITUENTS 

NMATS=4 

M=1 CMOD=l MATID=C 
M=2 CMOD=6 TREF=650 . MATID=B 
M=3 CMOD=9 MATID=A D=0 . , 0 . , 1 . 
M=4 CMOD=21 MATID=F 
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♦constituents : Constituent (subcell) materials 



Figure 2.1 Definition of the d vector for local transverse isotropy. 


Associated MAC/GMC 4.0 Example Problem(s): 

Example lc-d, Example 2a-d, Example 3a-i, Example 4a-d, Example 4f-h, Example 5e-f, Example 6a-c, 
Example 7a (see Table 2.2) 


Notes: 

1) TREF=£re/is only required in the special case when the code cannot determine at what temperature 
the material properties should be taken. This special case occurs when: 

• Thermal loading is not specified (* THERM omitted) 

-AND- 

• Material properties are taken from the internal material database ( matid ^ u) 

-OR- 

Material properties are user-defined ( matid = u) - AND - temperature-dependent 

2) In all other cases TREF=£re/is optional. 

Warning : Specifying TREF =tref will override the temperature dependence and cause the code to 
employ temperature-independent material properties for the material taken at the temperature tref. 

3) D=d 1 ,d 2 ,d 3 is only required in the special case when the specified material constitutive model is one 
that admits arbitrary local transverse isotropy. This occurs when: 

cmod= 3 - OR - cmod= 7 - OR - cmod= 9 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 4 (*SOLVER) 
MAC/GMC 4.0 Theory Manual Section 4 
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Table 2.2 MAC/GMC 4.0 constitutive models and internal material property database. Numerical 
values for the material properties can be obtained from the MAC/GMC output file. 


Model 

matid 

Material 

Stress 

Units* 

Temperature 

Dependent? 

Example 

Problem(s) 

Bodner-Partom 

(CMOD=l) 

A 

Aluminum (2024-T4) 

MPa 

Yes 

2a 

B 

Aluminum (2024-0) 

MPa 

Yes 

2a 

C 

Aluminum (6061-0) a 

MPa 

Yes 

2a, 3c 

D 

Aluminum (6061-0) b 

MPa 

Yes 

2a 

E 

Aluminum (pure) 

MPa 

Yes 

2a 

F 

Titanium (pure) 

MPa 

No 

2a 

G 

Copper 

MPa 

No 

2a 

U 

User-Defined 

- 

- 

2e,2f,2h 

Modified 

Bodner-Partom 

(CMOD=2) 

A 

Ti-218 

MPa 

Yes 

2b 

U 

User-Defined 

— 

— 

— 

Robinson 

Transversely 

Isotropic 

Viscoplastic 

(CMOD=3) 

A 

Kanthal 

ksi 

No (600° C) 


B 

FeCrAlY 

ksi 

Yes 


C 

35% Tungsten/Kanthal Composite 

ksi 

No (600° C) 

- 

U 

User-Defined 

— 

— 

— 

Isotropic GVIPS 

(CMOD=4) 

A 

Ti-21S 

ksi 

Yes 

lc, Id, 2b, 3a, 
3b,3d-i,4a- 
d,4f-h,5e, 
5f,6a-c,7a 

U 

User-Defined 

- 

- 

2e,3b 

Standard Linear 
Elastic 

(CMOD=6) 

A 

Boron Fiber 

MPa 

No 

3c 

B 

SiC (SCS-6) Fiber 

MPa 

Yes 

- 

C 

Tungsten Fiber 

MPa 

No 


D 

Boron Fiber 

ksi 

No 

3a, 3b 

E 

SiC (SCS-6) Fiber 

ksi 

Yes 

Id, 3a, 3b, 3d- 

i,4a-d,4f-h, 

5e,5f,6-c,7a 

F 

Tungsten Fiber 

ksi 

No 


U 

User-Defined 



la,lb,2e,2f, 

2h,4e,5a-d, 

7d,7e 

Transversely 
Isotropic GVIPS 

(CMOD=7) 

A 

Ti-6-4 

ksi 

Yes 

— 

U 

User-Defined 

— 

— 

— 

Transversely 
Isotropic Elastic 

(CMOD=9) 

A 

T50 Graphite 

MPa 

No 

— 

B 

T300 Graphite 

MPa 

No 


C 

PI 00 Graphite 

MPa 

Yes 

- 

D 

T50 Graphite 

ksi 

No 


E 

T300 Graphite 

ksi 

No 


F 

PI 00 Graphite 

ksi 

Yes 


U 

User-Defined 

- 

- 

7b-e 
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* constituents : Constituent (subcell) materials 


Table 2.2 (continued) 


Freed-Walker 

(CMOD=l 0) 

A 

Copper 

MPa 

Yes 

— 

B 

NARloy Z 

MPa 

Yes 


U 

User-Defined 

- 

- 

- 

Robinson 

Isotropic 

(CMOD=l l) 

A 

NARloy Z 

MPa 

Yes 

— 

U 

User-Defined 

— 

— 

— 

Anisotropic 

Elastic 

(CMOD=15) 

U 

User-Defined 



7e 

Classical 

Incremental 

Plasticity 

(CMOD=2 1) 

A 

Copper (bilinear) 

MPa 

Yes 

2c 

B 

Ti-24-11 (bilinear) 

MPa 

Yes 

2c 

C 

Ti-15-3 (point-wise) 

MPa 

Yes 

2c 

D 

Ti-24-11 (point-wise) 

MPa 

Yes 

2c 

E 

Copper (bilinear) 

ksi 

Yes 


F 

Ti-24-11 (bilinear) 

ksi 

Yes 


G 

Ti-15-3 (point-wise) 

ksi 

Yes 

- 

H 

Ti-24-11 (point-wise) 

ksi 

Yes 


U 

User-Defined 

- 

- 

2e,2f 

Multimechanism 

GVIPS + 

(CMOD=22) 

A 

Ti-21S 

ksi 

Yes 

7a 

Shape Memory 

(CMOD=30) 

A 

NiTi 

MPa 

Yes 

2d 

B 

NiTi 

ksi 

Yes 


U 

User-Defined 

- 

- 


User-Defined 

(CMOD=99) 

U 

User-Defined 

- 

- 

2g,2h 


* Time units are all seconds, temperature units are all °C, strain is unitless 

+ Non-U. S. government users may require an additional software license to use this constitutive model 


NAS A/TM— 2002-2 1 2077 /VOL2 


25 




MAC/GMC 4.0 Keywords Manual 


II. User-defined material properties ( matid = U) 

A. Material properties read directly from the input file ( matdb = 1) 

1. Standard constitutive models (cmod* 15, cmod* 99) 


Usage: 

a) Temperature-independent material properties 

*CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid matdb =matdb & 
EL =E A ,E T ,v A ,v T ,G A ,a T ,a A np =np vi=vi(l),vi(2),...,vi(n) D =d 1 ,d 2 ,d 3 
k=k a ,k t 


—> Repeated nmats times 

b) Temperature-dependent material properties 

CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid MATDB =matdb 
NTP =ntp 

tem =T(l),T(2),...,T(ntp) 

EA =E a (1), E a ( 2),..., E A (ntp) 

ET= E t (1), E t ( 2),..., E T (ntp) 

NUA= v a (1), v a (2), ..., v A (ntp) 
nut= v T ( 1 ), v T ( 2), v T (ntp) 

GA=G a (1),G a ( 2), ...,G A (ntp) 

RLPR=a A (l), a A (2), a A (ntp) 
alpt =cct(1), a T (2 ), .... a T (ntp) 

NP=«p 

vi =Vj(l), V,(2),..„ V,(ntp) 

V2=V 2 (1), V 2 (2),..„ V 2 (ntp) 

. — > Repeated n times, where n = number of visco-inelastic properties (see Table 2.3) 

Vn=V n (l), VJ2),..., V„(ntp) 
ka= k a ( ]), k a (2) ,...,k a ( ntp) 
ki=k t (1),k t ( 2) ,...,K T (ntp) 

'iJ = d i,(l 2 , f/j 


— > Repeated nmats times 
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♦constituents : Constituent (subcell) materials 


Name 

nmats 

matnum 

Definition 

Number of constituent materials specified 

Material number (data following M =matnum is associated with material number 
matnum) 

numbered sequentially 

cmod 

tref 1 ’ 2 

matid 

[=U] 

Constitutive model identification number 
Reference temperature for material (Special Case) 

Material i.d. letter, which indicates location of material parameters 
= a-h — » Internal material database (see Table 2.2) 

matdb 

[=1] 

= U —> User-defined material properties 
Material database location: 

= 1 — » Material properties read from input file 

= 2 — » Material properties determined from user-defined function (usrfun . F90) 

E a ,E r 3 
Va, Vt 3 

G a 3 

ccj, a a 
np 4 

= 3 — » Material properties read from external material database file 
Axial and transverse elastic (i.e., “Young’s”) modulus 
Axial and transverse Poisson ratios 
Axial shear modulus 

Axial and transverse coefficients of thermal expansion (CTEs) 

Number of stress-strain points employed to characterize post-yield portion of the stress- 
strain curve for classical plasticity (Special Case) 

vi(l-n) 5 

Visco-inelastic material properties (Special Case) 

n = number of visco-inelastic material properties and is cmod dependent (Table 2.3) 

d i, d 2 , d 2 

Components of vector normal to the plane of transverse isotropy (Special Case) (see 
Figure 2.1) 

Ka^t 1 

ntp 

Tfi) 

E A (i) 3 
Et d) 3 

VA0) 3 

v T (i) 3 
G a ( i) 3 
a A (i) 3 
ocr(i) 3 
E/0 5 

K/0 6 

k/0 6 

Axial and transverse thermal conductivities (Special Case) 

Number of temperatures at which the material properties are specified 

Temperatures at which the material properties are specified 

Axial elastic (i.e., “Young’s”) modulus at input temperature number i 

Transverse elastic (i.e., “Young’s”) modulus at input temperature number i 

Axial Poisson ratio at input temperature number i 

Transverse Poisson ratios at input temperature number i 

Axial shear modulus at input temperature number i 

Axial coefficient of thermal expansion (CTEs) at input temperature number i 
Transverse coefficients of thermal expansion (CTEs) at input temperature number i 
Visco-inelastic material property number j at input temperature number i (Special Case) 
Axial thermal conductivity at input temperature number i (Special Case) 

Transverse thermal conductivity at input temperature number i (Special Case) 
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Table 2.3 Temperature-independent user-defined inelastic material property specification. See the 


MAC/GMC 4.0 Theory manual for details on the constitutive models and their parameters. 



n 

Inelastic material property specification 

Plane of 
transverse 
isotropy 
specification 

Bodner-Partom 

(CMOD=l) 

6 

VI = D 0 , Zo, Z b m, n, q 

N/A 

Modified 

Bodner-Partom 

(CMOD=2) 

17 

VI = D 0 , Z 0 , Z b Z 2 , Z 3 , m ls m 2 , n, a u a 2 , r l5 r 2 , Dm,, 
Dm 2 , Dzi, Dz 2 , Dz 3 

N/A 

Robinson 

Transversely 

Isotropic 

Viscoplastic 

(CMOD=3) 

10 

VI = p, k t , R, H, n, m, P, q, co, G 0 

D=d],d 2 ,d 3 

Isotropic GVIPS 

(CMOD=4) 

11 

VI = p, k, R„, B 0 , B b n, p, q, Ko, B 0 ’, P 

N/A 

Standard Linear 
Elastic 

(CMOD=6) 

N/A 

N/A 

Elastic Only 

N/A 

Transversely 
Isotropic GVIPS 

(CMOD=7) 

10 

VI = k, n, p, m, P, R, H, G 0 , eo, q 

D =d h d 2 ,d 3 

Transversely 
Isotropic Elastic 

(CMOD=9) 

N/A 

N/A 

Elastic Only 

D=d],d 2 ,d 3 

Freed-Walker 

(CMOD=l 0) 

14 

VI = A, C, n, Q, T m , T t , 5, f, h, 1, D 0 , s p \ , ^ll , 

11/ Us 

D L 

N/A 

Robinson Isotropic 

(CMOD=l l) 

9 

VI = A, n, m, P, H, R, K 2 , Ko, Go 

N/A 

Anisotropic 

Elastic 

(CMOD=15) 

N/A 

N/A 

Elastic Only, Special Case - see below 

N/A 

Classical Plasticity 

(CMOD=2l) 

np + 1 

NP —tip VI — CJy, C7i, O/? ••• ? ^Tnp? £l? £*2? ••• ? £np 

(a y is the yield stress and a i? 8i are stress-total strain 
pairs that characterize the post-yield response) 

N/A 

Multimechanism 

GVIPS 

(CMOD=22) 

N/A 

N/A 

Model not currently supported 
with user-defined properties 

N/A 

Shape Memory 

(CMOD=30) 

5 

vi = f T , Y, a, a, n 

N/A 

User-Defined 

(CMOD=99) 

N/A 

N/A 

Special Case - see below 

N/A 
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♦constituents : Constituent (subcell) materials 


Examples: 

a) Temperature-independent material properties 

*CONSTITUENTS 

NMATS=2 

M=1 CMOD=6 MATID=U MATDB=1 & 

EL=388.2E9,7.6E9,0.41,0.45,14.9E9,-0.68E-6, 9.74E-6 
K=500 . , 10. 

M=2 CMOD=6 MATID=U MATDB=1 & 

EL=3 . 45E9, 3 . 45E9, 0.35, 0.35, 1 . 278E9, 45 . E-6, 45.E-6 
K=0 .19,0.19 


CONSTITUENTS 

NMATS=5 

M=1 CMOD=l MATID=U MATDB=1 & 

EL=9 . 53E3, 9. 53E3, 0 .33,0.33,3. 58E3, 21 . 06E-6, 21 . 06E-6 & 

VI=1.E4,49.,63.,300.,4.,1. 

M=2 CMOD=9 MATID=U MATDB=1 & 

EL=2 53 . 5E 9 , 6 . 05E 9 , 0 . 3 9 , 0 . 4 7 , 4 . 1 67E9 , -0 . 4 72 4E- 6 , 2 6 . 63E- 6 D=0.5,l.,0. 
M=3 CMOD=4 MATID=U MATDB=1 & 

EL=14009 . , 14009. , 0.365, 0.365, 5131 .5, 5. 862E-6, 5. 862E-6 & 

VI =0 . 000999275, 44 . 960, 1. 679E-07, 2 . 494561E-05, 0 . 05, & 

3. 3. 1. 8. 1. 35. 0. 85. 3. 018 3E -7, 0.001 
M=4 CMOD=l 0 MATID=U MATDB=1 & 

EL=0 . 1305E + 0 6 , 0 . 1305E+0 6, 0.36,0.36,0. 4 7 98E+05, 0 . 165E-04, 0 . 165E-04 & 
VI=0 . 2000E+08, 0 . 1300E+02 , 0 . 4500E+01, 0 . 2000E+06, 1082.85, & 

404 .85. 0. 3500E-01, 0.5000E+00, 0.4000E+02, 0. 6500E+01, & 

0 . 1300E+00, 0 . 1000E-01 , 0 . 1000E-06, 0 . 7000E+00 

M=5 CMOD=2 1 MATID=U MATDB=1 & 

EL=10000 . , 10000. , 0.326, 0.326, 3770. , 12.00E-6, 12.00E-6 & 

NP=1 VI =5 .08,13.1,0.15 

b) Temperature-dependent material properties 

CONSTITUENTS 

NMATS=2 

M=1 CMOD=2 1 MATID=U MATDB=1 
NTP=3 

TEM=2 1 . , 400 . , 800 . 

E A= 10000. , 9000. ,7300. 

ET=1 0 0 0 0 . , 9000 . , 7300 . 

NUA= 0.326, 0.351, 0.345 
NUT=0 .326, 0.351, 0.345 
GA=3771 .,3331., 2714. 

ALPA=12 . 00E-6, 13 . 50E-6, 22 . 72E-6 
ALPT=12 . 00E-6, 13 . 50E-6, 22 . 72E-6 
NP=3 

VI =4. 20, 3. 81, 3. 05 
V2=6 . 13 , 4.77, 3. 69 
V3=6.60,5.10,3.90 
V4=13.1, 9.70,5.10 
V5=0 . 004 , 0.004, 0.004 
V6=0 . 01 , 0.01, 0.01 
V7=0 .15,0.15,0.15 
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M=2 CMOD=l MATID=U MATDB=1 
NTP=2 

TEM=18 . , 700 . 

EA=9 . 53E3 , 4 . 12E3 
ET=9.53E3,4.12E3 
NUA=0 .41,0.41 
NUT=0 .41,0.41 
GA=3 . 58E3 , 1 . 46E3 
ALPA=21 . 06E-6, 28 . 92E-6 
ALPT=2 1 . 06E-6, 28 . 92E-6 
Vl=l ,E4, 1 ,E4 
V2=49. , 49. 

V3=63 . , 63. 

V4=300 . , 300 . 

V5=4 .,2.5 
V6=l . , 1 . 

Associated MAC/GMC 4.0 Example ProblemfsR 

a) Temperature-independent material properties 

Example la, Example lb, Example 2e, Example 3b, Example 7e 

b) Temperature-dependent material properties 

Example 2e, Example 4 e, Example 5a, Example 5b, Example 5c, Example 5d 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 4 (*solver) 
MAC/GMC 4.0 Theory Manual Section 4 


2. Special cases 

a) Interna] anisotropic constitutive model (cmod= 15) 


Usage: 

i) Temperature-independent material properties 

*CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid matdb =matdb & 

EL=C//, Cj2, C]3, C 14, Ci 5, Ci & C22, C23, C24, C25, C26, C33, C34, C35, C36, C44, C45, C46, C55, C56, Cm, & 

ocn, CC22, OC33, a 2 3, ax, a n 


— » Repeated nmats times 
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* constituents : Constituent (subcell) materials 


ii) Temperature-independent material properties 

^^CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid matdb =matdb 
NTP =ntp 

tem=T( 1),T(2), ...,T(ntp) 

Cll =C 11 (1),C]](2), ...,C]j(ntp) 

C 1 2=C 12 (1), C 12 (2), ..., C 12 (ntp) 

C 1 3=Cis( 1), Cj 3 (2), ...,C 12 (ntp) 

C 1 4 =CuC 1), c 14 (2), ..., C 14 (ntp) 

C15 =C 15 (1),C 15 (2), ...,C 15 (ntp) 

Cl <o=C i6(l), C 16(2), ...,C]6(ntp) 

C 2 2=C 22 ( 1), C 22 (2), ..., C 22 ( ntp) 

C 2 3=C 23 ( 1), C 23 (2), ...,C 23 (ntp) 

C 2 4 =C 24 ( 1 ), C 24 (2), ..., C 24 (ntp) 

C 2 5 = C 2 j( 1), C 25 ( 2),..., C 23 ( ntp) 

C2 6 = C 26 ( 1), C 2 $( 2), ...,C 2 6(ntp) 

C 3 3 = C33 (1), C33 (2), ..., C33 (n tp) 

C 3 4=C 34 (1 ), C 34 (2), ..., C 34 (ntp) 

C3S=C 35 (1),C 35 (2), ... ,C 35 (fitp) 

C3 6=C 3 6(1),C 3(3 (2), ...,C 3 ^(ntp) 

C 4 4 =C 44 ( 1 ), C 44 (2), ..., C 44 (ntp) 

C4 S=C 43 ( l ), C 45 (2), ... ,C 43 (ntp ) 

C4 6 =C 46 (1),C 46 (2), ...,C 46 (ntp) 

C55= C 55 ( 1), C 55 ( 2),..., C 55 ( ntp) 

Cb6=C 56 (l),C 56 (2), ...,C 3 ^(ntp) 

C 6 6 =C 66 ( l), C 66 (2), ..., C 66 (ntp) 

ALFl^a 11 (l),a 11 (2), ...,a n (ntp) 

ALF2 =a 22 (l),a 22 (2),...,a 22 (ntp) 

KLE2=a 33 (l), a 33 (2), ..., a 33 (ntp) 
al F 4 = a 23 (l), a 23 (2), ..., a 23 (ntp) 

KLEb=a 13 (l), a 13 (2), ..., a 13 (ntp) 

ALF6 =a 12 (l), a 12 (2), ..., a 12 (ntp) 


— » Repeated nmats times 


Name Definition 

nmats Number of constituent materials specified 

matnum Material number (data following M =matnum is associated with material number 
matnum) 

numbered sequentially 

cmod Constitutive model identification number 

[=15] 

tref 1,2 Reference temperature for material (Special Case) 
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Name Definition 

matid Material i.d. letter, which indicates location of material parameters 
= a-h — » Internal material database (see Table 2.2) 

= U — » User-defined material properties 

matdb Material database location: 

[ =1 ] = 1 — » Material properties read from input file 

= 2 Material properties determined from user-defined function (usrfun . F90) 
= 3 — » Material properties read from external material database file 
Cy Stiffness matrix coefficients 

ay 3 CTE vector components (anisotropic materials can have 6 CTE components) 

ntp Number of temperatures at which the material properties are specified 

T(i) Temperatures at which the material properties are specified 

Cjk(i) Stiffness matrix component jh at input temperature number i 

a Jk (i) 3 CTE component jh at input temperature number i 


Examples: 

i) Temperature-independent material properties 


CONSTITUENTS 

NMATS=1 

M=1 CMOD=15 MATID=U MATDB=1 
EL=6. 677E9, 3.274E9, 


4.843E9, 

0 ., 

-0.1887E9, 

0 . , & 

3.384E9, 

0 ., 

0.006711E9, 

0 . , & 

55.80E9, 

0 

-6.96 9E 9 , 

0 . , & 


2 . 625E9, 

o.. 

-0 . 1830E9, 



2.269E9, 

0 . , & 


1.613E9 & 


4 3 . 7 4E- 6 , 3 9 . 04E- 6 , 4.043E-6, 0., 12.19E-6, 0. 


ii) Temperature-dependent material properties 

CONSTITUENTS 

NMATS=1 

M=1 CMOD=15 MATID=U MATDB=1 
NTP=2 

TEM=2 3 . f 250 . 

C11=6.677E9,5.973E9 
C12=3.274E9, 2 . 931E9 
C13=4 . 843E9, 4 . 517E9 
C14=0. ,0. 

C15=0 . 1887E9, 0 . 1932E9 
Cl 6=0 ., 0 . 

C22=6.662E9,5.867E9 
C23=3.384E9,3.023E9 
C24=0. ,0. 

C25=-0 . 006711E9, 0 . 008475E9 
C26=0. f 0. 

C33=55.80E9 f 48.22E9 
C34=0. , 0. 

C35=6.969E9,7.105E9 
C3 6=0 .,0. 

C4 4=2 . 625E9, 2 . 124E9 
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C45=0. , 0. 

C46=0.1830E9, 0.1363E9 
C55=2.269E9,2.123E9 
C56=0 . , 0 . 

C66=l . 613E9, 1 . 444E9 
ALFl=43.74E-6, 46.75E-6 
ALF2=39 . 04E-6, 43.81E-6 
ALF3=4 . 043E-6, 5 .241E-6 
ALF4=0 . , 0 . 

ALF5=-12.19E-6,-16.45E-6 
ALF6=0 . , 0 . 


Associated MAC/GMC 4.0 Example ProblemfsL 

Example 7e 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 4 (*SOLVER) 


b) User-defined constitutive model (cmod= 99) 


Usage: 

i) Temperature-independent material properties 

CONSTITUENTS 

nmats =nmats 

n=matnum cmod =cmod TREE=tref matid =matid matdb =matdb npe =npe np v=npv & 
el =E(l),E(2),...,E(npe) alp a A ,a T vi=V(l),V(2),...,V(npv) 
k=k a ,k t 


—> Repeated nmats times 

ii) Temperature-dependent material properties 

^^CONSTITUENTS 

NMATS =nmats 

n=matnum cmod =cmod tref =tref matid =matid matdb =matdb 
NPE =npe NP v=npv 
ntp =ntp 

tem =T(l),T(2),...,T(ntp) 
e^E^I), Ei( 2), Ei(ntp) 

E2= E 2 (l), E 2 (2), E 2 (ntp) 

. — > Repeated npe times 
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E nvp E nvp (l), E nvp (2), ... , E nvp (ntp) 
ALPA=au(X), &a(2), • ••, oc A (ntp) 
ALPT=<2 r (7), 6^(2), a T (ntp) 
vi=V!(l), Vj(2),..., Vj(ntp) 
V2=V 2 (1), V 2 (2),..„ V 2 (ntp) 

. — » Repeated npv times 

v 12 vp- V nvp (l), V nvp (2), ... , V nvp (ntp) 
KA =k a (1),k a ( 2) ,...,K A (ntp) 

KT =k t (1),k t (2) , ...,K T (ntp) 


— > Repeated nmats times 


Name 

nmats 

matnum 


cmod 
[=99] 9 ’ 10 
tref U2 
matid 
[=U] 


matdb 

[= 1 ] 


npe 

4 

npv 

E(l-npe) 
a A a T 
V(1 -npv) 4 

Ka^t 1 

ntp 

T(i) 

E/i ) 
a A (i) 
a T (i) 

Vj(i) 11 
K A (i) 1 
KtO) 1 


Definition 

Number of constituent materials specified 

Material number (data following M =matnum is associated with material number 
matnum) 

numbered sequentially 

Constitutive model identification number 

Reference temperature for material (Special Case) 

Material i.d. letter, which indicates location of material parameters 
= a-h — » Internal material database (see Table 2.2) 

= U — » User-defined material properties 

Material database location: 

= 1 — » Material properties read from input file 

= 2 — » Material properties determined from user-defined function (usrfun . F90) 

= 3 — » Material properties read from external material database file 
Number of elastic constants employed in user-defined constitutive model 
Number of visco-inelastic constants employed in user-defined constitutive model 
Elastic constants 

Axial and transverse coefficients of thermal expansion (CTEs) 

Visco-inelastic constants 

Axial and transverse thermal conductivities (Special Case) 

Number of temperatures at which the material properties are specified 

Temperatures at which the material properties are specified 

Elastic material property number j at input temperature number i 

Axial coefficient of thermal expansion (CTEs) at input temperature number i 

Transverse coefficients of thermal expansion (CTEs) at input temperature number i 

Visco-inelastic material property number j at input temperature number i 

Axial thermal conductivity at input temperature number i (Special Case) 

Transverse thermal conductivity at input temperature number i (Special Case) 
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♦constituents : Constituent (subcell) materials 


Examples: 

i) Temperature- independent material properties 

♦CONSTITUENTS 

NMATS=2 

M=1 CMOD=99 MATID=U MATDB=1 NPE=2 NPV=2 & 

EL=55.2E9, 0.30 ALP=22 . 5E-6, 22 . 5E-6 VI=1 . 5E-28, 3 . 0 
K=2 7 3 . , 273 . 

M=2 CMOD=99 MATID=U MATDB=1 NPE=2 NPV=6 & 

EL=55.2E9, 0.30 ALP=22 . 5E-6, 22 . 5E-6 & 

VI=1000. , 103 . 42E6, 103 . 42E6, 1700 . , 10 . , 1 . 0 
K=273 . , 273 . 

ii) Temperature-dependent material properties 

♦CONSTITUENTS 

NMATS=1 

M=1 CMOD=99 MATID=U MATDB=1 

NPE=2 NPV=2 

NTP=3 

TEM=2 1 . , 350 . , 650 . 

E1=55.2E9,52.7E9,47.1E9 
E2=0. 30, 0.30, 0.30 
ALPA=22 . 5E-6, 23 . 4E-6, 27 . OE-6 
ALPT=22 . 5E-6, 23 . 4E-6, 27 . OE-6 
VI =1 . 5E-28, 1 . 5E-28, 1 . 5E-28 
V2=3.0, 3.017, 3.052 


Associated MAC/GMC 4.0 Example Problem (s): 

Example 2g 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Appendix 


Notes: 

1) TREF=tref is only required in the special case when the code cannot determine at what temperature 
the material properties should be taken. This special case occurs when: 

• Thermal loading is not specified (* THERM omitted) 

-AND- 

• Material properties are taken from the internal material database ( matid ^ u) 

-OR- 

Material properties are user-defined ( matid = u) - AND - temperature-dependent 
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2) In all other cases TREF=*ref is optional. 

Warning : Specifying TREF =tref will override the temperature dependence and cause the code to 
employ temperature-independent material properties for the material taken at the temperature tref. 

3) To obtain isotropic material behavior, set, E A = E T} v A = v T , G A = E A / 2(l+v A ), and a A = a T . 

4) NP =np is only required in the special case that the material constitutive model is classical incremental 
plasticity (cmod= 2 1). A maximum value of NP=9 may be specified. 

5) The lines associated with the visco-inelastic material properties (vi= vi(l Vj= V/i),...) are only 
required in the special case when the constitutive model (cmod) chosen is an inelastic model. This 
occurs when: 

cmod^ 6 AND cmod^ 9 AND cmod^l 5 

6) D=dj,d 2 ,d 3 is only required in the special case when the specified material constitutive model is on e 
that admits local transverse isotropy. This occurs when: 

cmod= 3 - OR - cmod- 7 - AND - cmod= 9 

7) The thermal conductivity specifiers (k =k a ,k t , ka=k a ,..., kt=k t , ...) are only required in the special 
case when the code is calculating effective thermal conductivities of the RUC. This occurs when the 
* CONDUCTIVITY keyword is specified. The thermal conductivity specification is only valid for 
materials that are isotropic or transversely isotropic with an x 2 -x 3 plane of isotropy (in the global RUC 
coordinates). 

8) Shear CTEs are engineering quantities. That is, for example, = a 23 AT . 

9) When a user-defined constitutive model is selected (cmod=9 9), MAC/GMC 4.0 will use the code in 
the file usrmat . F90 to determine the inelastic strain increments for any subcells contain that 
particular material. The MAC/GMC 4.0 executable that is distributed has been linked with the 
demonstration usrmat subroutine, the source code of which is also distributed with MAC/GMC 4.0 
(see Appendix). In order for users to employ their own usrmat subroutine, it is necessary to place 
the appropriate code in the usrmat. F90 file and to compile the subroutine and link with the 
MAC/GMC 4.0 library file (mac4.1ib), which is also distributed with MAC/GMC 4.0. This 
procedure results in a new MAC/GMC 4.0 executable file that is now linked to the user’s own 
usrmat subroutine rather than the demonstration version of usrmat distributed with the code. It is 
the user’s responsibility to ensure that the coding, compilation, and linking of the user’s own usrmat 
subroutine is done correctly and that the desired user-defined constitutive model is actually being 
employed when the code is executed. A printout of the demonstration usrmat subroutine that is 
distributed with MAC/GMC 4.0 is given in the Appendix of this manual. Note that when the code 
executes the demonstration version of the usrmat subroutine, a message is printed to the screen. 

10) When a user-defined constitutive model is selected (CMOD=9 9), MAC/GMC 4.0 will use the code in 
the file usrformde . F90 to calculate the material elastic stiffness matrix from the user-defined 
elastic properties. This is necessary because the user is free to employ any elastic constants desired 
(e.g., Lame constants, compliance components). Thus the code does not (in general) know how to 
form the stiffness matrix from the user’s elastic constants. As with the usrmat file, the MAC/GMC 
4.0 executable that is distributed has been linked with a demonstration usrformde subroutine, the 
source code of which is also distributed with MAC/GMC 4.0 (see Appendix). It is the user’s 
responsibility to ensure that the coding, compilation, and linking of the user’s own usrformde 
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* constituents : Constituent (subcell) materials 


subroutine is done correctly and that the stiffness matrix of the user’s material is being calculated 
correctly when the code is executed. A printout of the demonstration usrformde subroutine that is 
distributed with MAC/GMC 4.0 is given in the Appendix of this manual. Note that when the code 
executes the demonstration version of the usrformde subroutine, a message is printed to the screen. 

1 1) For a user-defined elastic constitutive model, set npv = 0 and omit the VI - Vnvp lines. 


B. Material properties determined from a user-defined function ( matdb = 2) 


Usage: 

♦CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid MATDB =matdb 


—> Repeated nmats times 

Name Definition 

nmats Number of constituent materials specified 

matnum Material number (data following M =matnum is associated with material number 
matnum) 

numbered sequentially 

cmod Constitutive model identification number 

tref 1,2 Reference temperature for material (Special Case) 

matid Material i.d. letter, which indicates location of material parameters 

= a-h Internal material database (see Table 2.2) 

= U — » User-defined material properties 

matdb Material database location: 

[ =2 ] 3,4 = 1 Material properties read from input file 

= 2 —> Material properties determined from user-defined function (usrfun.F90) 

= 3 — » Material properties read from external material database file 

Example: 

♦CONSTITUENTS 

NMATS=5 

M=1 CMOD=6 MATID=U MATDB=2 
M=2 CMOD=21 MATID=U MATDB=2 
M=3 CMOD=l TREF=650 . MATID=U MATDB=2 
M=4 CMOD=l TREF=650 . MATID=U MATDB=2 
M=5 CMOD=99 MATID=U MATDB=2 


Associated MAC/GMC 4.0 Example ProblemfsU 

Example 2f 
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Notes: 

1) TREF=tref is only required in the special case when the code cannot determine at what temperature 
the material properties should be taken. This special case occurs when: 

• Thermal loading is not specified (* THERM omitted) 

-AND- 

• Material properties are taken from the internal material database ( matid ^ u) 

-OR- 

Material properties are user-defined ( matid = u) - AND - temperature-dependent 

2) In all other cases TREF=£re/is optional. 

Warning : Specifying TREF =tref will override the temperature dependence and cause the code to 
employ temperature-independent material properties for the material taken at the temperature tref. 

3) When a user-defined function is selected to determine a material’s properties (matdb=2), MAC/GMC 
4.0 will use the code in the file usrfun.F90 to determine the material properties at every increment 
of the simulated applied loading. The MAC/GMC 4.0 executable that is distributed has been linked 
with the demonstration us r fun subroutine, the source code of which is also distributed with 
MAC/GMC 4.0 (see Appendix). In order for users to employ their own usrfun subroutine, it is 
necessary to place the appropriate code in the usrfun. F90 file and to compile the subroutine and 
link with the MAC/GMC 4.0 library file (mac4 . lib), which is also distributed with MAC/GMC 4.0. 
This procedure results in a new MAC/GMC 4.0 executable file that is now linked to the user’s own 
usrfun subroutine rather than the demonstration version of usrfun distributed with the code. It is 
the user’s responsibility to ensure that the coding, compilation, and linking of the user’s own usrfun 
subroutine is done correctly and that the desired user-defined function is actually being employed 
when the code is executed. A printout of the demonstration usrfun subroutine that is distributed 
with MAC/GMC 4.0 is given in the Appendix of this manual. Note that when the code executes the 
demonstration version of the usrfun subroutine, a message is printed to the screen. The usrfun 
subroutine can be used to allow the elastic and inelastic material properties to be a function of 
temperature or another variable. It works in concert with the MAC/GMC internal constitutive models 
as well as the user-defined constitutive model subroutine (usrmat). 

4) In MAC/GMC 4.0, the material properties are associated with constituent materials as opposed to 
subcells. That is, several subcells may contain material number 1, and thus all of those subcells must 
have the same material properties. This is relevant when employing the usrfun . F90 subroutine to 
make material properties a function of field variables (as is done in Example Problem 2f). If, for 
instance, the properties of material number 1 are determined from a function of strain and material 
number 1 is placed in two subcells, the subcells will in general experience different strain states. It is 
thus impossible for material number 1 to have the correct material properties to correspond to the two 
different states of strain in each subcell. For this reason, when material properties are determined 
from a function of field variables (other than temperature, which does not vary within the repeating 
unit cell), each subcell should be given its own material. 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Appendix 
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Section 2: Material Keywords 


♦constituents : Constituent (subcell) materials 


C. Material properties read from external material database file ( matdb = 3) 


Usage: 

♦CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF=tref MATID =matid MATDB =matdb EX r mAT=matname 


Repeated nmats times 


Name 

nmats 

matnum 


cmod 
tref x ' 2 
matid 

[=u] 


matid 


[=u] 


matdb 
[=3] 3 


matname 4 


Definition 

Number of constituent materials specified 

Material number (data following M =matnum is associated with material number 
matnum) 

numbered sequentially 

Constitutive model identification number 

Reference temperature for material (Special Case) 

Material i.d. letter, which indicates location of material parameters 
= a-h — » Internal material database (see Table 2.2) 

= U — » User-defined material properties 
Material i.d. letter from internal material database 

Material database location: 

= 1 —> Material properties read from input file 

= 2 — » Material properties determined from user-defined function (usrfun.F90) 

= 3 — » Material properties read from external material database file 
Name of material in external material database file 


Examples: 

CONSTITUENTS 

NMATS=3 

M=1 CM0D=1 MATID=U MATDB=3 EXTMAT=copper 
M=2 CM0D=9 9 MATID=U MATDB=3 EXTMAT=A1 (6061-0) 

M=3 CMOD=6 MATID=U MATDB=3 EXTMAT=SCS-6 

Example portion of the external material database file : 

# 

# Bodner-Partom Model Material 

MAT=copper 

CMOD=l 

EL=120.E9,120.E9,0.33,0.33,45.11E9,14.7E-6,14.7E-6 & 
VI=1 . E4 , 63.E6,250.E6, 8.19,7.5, 0.55 
K=1 ., 2 . 
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Associated MAC/GMC 4.0 Example ProblemfsR 

Example 2h 


Notes: 

1) TREF=tref is only required in the special case when the code cannot determine at what temperature 
the material properties should be taken. This special case occurs when: 

• Thermal loading is not specified (* THERM omitted) 

-AND- 

• Material properties are taken from the internal material database ( matid ^ u) 

-OR- 

Material properties are user-defined ( matid = u) - AND - temperature-dependent 

2) In all other cases TREF=£re/is optional. 

Warning : Specifying TREF =tref will override the temperature dependence and cause the code to 
employ temperature-independent material properties for the material taken at the temperature tref. 

3) In order to use the external material database feature of MAC/GMC 4.0, the *mdbpath keyword 
must be specified (see Section 1 of this manual). It is there that the name of the external material 
database file is specified. 

4) The format of the material property specification in the external material database file is very similar 
to that employed in the MAC/GMC 4.0 input file for user-defined material properties. A printout of 
the sample external material database file distributed with MAC/GMC 4.0 is given in the appendix of 
this manual. The same special characters (“#” for a comment, for a continuation) that can be 
used in the input file can be used in the external material database file, and the first line of the file is 
interpreted as a title. Each material in the file is specified by a name using the line: 

MAT =matname 

where matname is the name of the material. The matname specified in the input file should 
correspond to one (and only one) of the material names in the external material database file. The 
next line for each material re-specifies the material constitutive model: 

CMOD=cmod 

The code uses this as a consistency check as it must match the value of cmod specified for the 
material in the input file. The remainder of the material property input data in the external material 
database file follows the format of the user-defined material properties that read directly from the 
input file (see outline section II. A. above) 

Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 2 (*mdbpath) 
MAC/GMC 4.0 Keywords Manual Appendix 
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♦constituents : Constituent (subcell) materials 


III. Electromagnetic (smart) material properties 


Usage: 

A) Temperature-independent electromagnetic material properties 

^^CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid MATDB =matdb EM =em & 
E L —Ea> Ej 1 , Va> Vt> Ga, &a d = d i , d 2 , d 3 
ES=ejj,ej2,e 2 6 
QS=qn,qi2,q 2 6 

KS = fCjj, k 22 
AS=( 2 1 1 , U 22 
MS fA 22 
PELS=^/, £ 

pmgs=^, y / > 


—> Repeated nmats times 

B) Temperature-dependent electromagnetic material properties 

^^CONSTITUENTS 

NMATS =nmats 

M=matnum CMOD =cmod TREF =tref MATID =matid MATDB =matdb EM =em 
NTP =ntp 

tem =T(l),T(2),...,T(ntp) 
ea=E a (1), E a ( 2), ..., E A (ntp) 

ET= E t (1), E t ( 2), E T (ntp) 

NUA= v A (l), v a ( 2), ..., v A (ntp) 

NUT= v T (l), v t (2), ..., v T (ntp) 

GA=G a (1),G a ( 2), ... ,G A (ntp j 
A’LVA=a A (l), a A (2), ..., a A (ntp) 
alpt =a T (l), gct(2), .... a T (ntp) 

D~cl j.clj, d] 

ES11 =e 1I (l),e 1I (2),...,e 1I (ntp) 

E s 1 2 =e I2 ( 1), e I2 (2), ..., e I2 (ntp) 

ES2 6=e 26 (l),e 26 (2), ...,e 26 (ntp) 

QS11 =q I1 (l),qn(2),...,q 1 i(ntp) 

QS12 ^q 12 (l),q 12 (2),...,q I2 (ntp) 

Q S 2 6 =q 26 ( l),q 26 (2),...,q 26 (ntp) 

KS1 1 = hCi [( 1 ), K[ [(2), .... K[ [(ntp) 

KS 2 2 =k 22 (1), k 22 (2), ...,K 22 (ntp) 

AS 11 =a 1I (l),a I i(2), ...,au(ntp) 

AS22 =a 22 (l),a 22 (2), ...,a 22 (ntp) 

MS 1 1 =p, i(l),p n (2),..., p, , (ntp) 

MS 2 2 =a 22 (l), p 22 (2), ...,p 22 (ntp) 
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PELSl =£ I (l),£i(2),...,£ I (ntp) 
pels 2=C 2 (1),C 2 (2), ...,Q(ntp) 
PMGS1 =y/i(l), y/i(2), i// i (nip) 
PMGS2 = y/ 2 (l), V2(2), y/ 2 (ntp) 


— » Repeated nmats times 


Name 

nmats 

matnum 


cmod 

[ = 9] 

tre/ 2 ’ 3 

matid 

[=U] 


matdb 
[=i] 1 


em 

[ = 1 ] 4 


Ea,Et 

Va,Vt 

G a 

CCj, CC A 
djf d 2 , d 2 

&kj 

Vkj 


Mu 

Ck 

Vk 

Ea(p) 
Et(p) 
Va(p) 
v T (p) 
Ga(p) 
cc A (p) 
a T (p) 
d i, d 2} d 2 

e kj(p) 

<lkj(p) 


Definition 

Number of constituent materials specified 

Material number (data following M =matnum is associated with material number 
matnum) 

numbered sequentially 

Constitutive model identification number 

Reference temperature for material (Special Case) 

Material i.d. letter, which indicates location of material parameters 
= a-h — » Internal material database (see Table 2.2) 

= U — » User-defined material properties 

Material database location: 

= 1 —> Material properties read from input file 

= 2 — » Material properties determined from user-defined function (usrfun . F90) 

= 3 — » Material properties read from external material database file 
Electromagnetic specifier for material: 

= 0 — » Material is not electromagnetic 
= 1 — » Material is electromagnetic 
Axial and transverse elastic (i.e., “Young’s”) modulus 
Axial and transverse Poisson ratios 
Axial shear modulus 

Axial and transverse coefficients of thermal expansion (CTEs) 

Components of vector normal to the plane of transverse isotropy (see Figure 2.1) 

Material piezoelectric component kj 

Material piezomagnetic component kj 

Material dielectric component ij 

Material magnetoelectric component ij 

Material magnetic permeability component ij 

Material pyroelectric constant component k 

Material pyromagnetic constant component k 

Axial elastic (i.e., “Young’s”) modulus at input temperature number p 
Transverse elastic (i.e., “Young’s”) modulus at input temperature number p 
Axial Poisson ratio at input temperature number p 
Transverse Poisson ratios at input temperature number p 
Axial shear modulus at input temperature number p 

Axial coefficient of thermal expansion (CTEs) at input temperature number p 
Transverse coefficients of thermal expansion (CTEs) at input temperature number p 
Components of vector normal to the plane of transverse isotropy (see Figure 2.1) 
Material piezoelectric component kj at input temperature number p 
Material piezomagnetic component kj at input temperature number p 
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* constituents : Constituent (subcell) materials 


Name Definition 

Ky(p) Material dielectric component ij at input temperature number p 
cii/p) Material magnetoelectric component ij at input temperature number p 

/Uij(p) Material magnetic permeability component ij at input temperature number p 
dj k (p) Material pyroelectric constant component k at input temperature number p 

Yk(p) Material pyromagnetic constant component k at input temperature number p 

Examples: 

A) Temperature-independent electromagnetic material properties 

♦CONSTITUENTS 

NMATS=1 

M=1 CMOD=9 MATID=U MATDB=1 EM=1 & 

EL=143.57E9,154.57E9,0.37,0.368,45.3E9,0.00E-6,0.00E-6 & 

D=0 . , 0 . , 1 . 

ES=0 . 0, 0 . 0, 0 . 0 
QS=699.7, 580.3, 550. 

KS=0 . 93E-10, 0 . 08E-9 
AS=0 .0,0.0 
MS=157 . E-6, -590 .E-6 
PELS=0 .0,0.0 
PMGS=0 .0,0.0 

B) Temperature-dependent electromagnetic material properties 

♦CONSTITUENTS 

NMATS=2 

M=1 CMOD=9 MATID=U MATDB=1 EM=1 
NTP=2 

TEM=24 . , 600 . 

EA=1 1 1 . 93E9, 101 . 43E9 
ET=116.33E9,107.10E9 
NUA=0 . 321, 0 . 321 
NUT=0 . 307,0 . 307 
GA=43 . 0E9, 41 . 5E9 
ALPA=1 . 99E-6, 2 . 55E-6 
ALPT=8 . 53E-6, 10 . 80E-6 
D=0 . , 0 . , 1 . 

ESI 1=1 8 .6,18.6 
ES12=-4 . 4, -4 . 4 
ES26=11 .6,11.6 
QS11=0 .0,0.0 
QS12=0 .0,0.0 
QS26=0.0, 0.0 
KS11=12 . 6E-9, 12 . 6E-9 
KS 22=11 . 2E-9, 11 .2E-9 
AS11=0 .0,0.0 
AS22=0 .0,0.0 
MSI 1=10 . OE-6, 10 . OE-6 
MS22=5 . OE-6, 5 . OE-6 
PELS1=0 . 13E5, 0 . 13E5 
PELS2=0 . 13E5, 0 . 13E5 
PMGS1=0 .0,0.0 
PMGS2=0 .0,0.0 
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Associated MAC/GMC 4.0 Example ProblemfsV 

Example 7b, Example 7c 

Notes: 

1) In order to perform coupled thermo-electro-magneto-elasto-plastic analysis of smart composites and 
laminates, the *electromag keyword must be specified (see Section 1 of this example manual). 
Further, the coupled analysis requires the constituent materials to have certain additional 
electromagnetic properties. Currently, specification of these properties is limited in MAC/GMC 4.0 
to the case in which: 

a) The material constitutive model is the transversely isotropic elastic model (cmod= 9) 

-AND- 

b) The material properties are user-defined (matid= u) 

-AND- 

c) The material properties are read from the input file (matdb= l) 

Thus, the specification of constituent electromagnetic properties actually falls within section II. A. of 
this outline rather than in its own section III. The present format was chosen for simplicity and so as 
to make the electromagnetic property section more evident to the reader. 

2) TREF=tref is only required in the special case when the code cannot determine at what temperature 
the material properties should be taken. This special case occurs when: 

• Thermal loading is not specified (* THERM omitted) 

-AND- 

• Material properties are taken from the internal material database ( matid ^ u) 

-OR- 

Material properties are user-defined ( matid = u) - AND - temperature-dependent 

3) In all other cases TREF=£re/is optional. 

Warning : Specifying TREF =tref will override the temperature dependence and cause the code to 
employ temperature-independent material properties for the material taken at the temperature tref. 

4) The electromagnetic specifier (EM =em) is required to indicate if a particular material is 
electromagnetically active (i.e., if it has electromagnetic properties). This allows electromagnetically 
active and inactive constitutents to both be present in a composite or laminate. 

5) E=d 1 ,d 2 ,d 3 is required when the material is electromagnetically active (em= 1) because, in such a case, 
the transversely isotropic elastic model must be employed (cmod= 9). 

Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 2 (*electromag) 
MAC/GMC 4.0 Keywords Manual Section 4 (*laminate) 
MAC/GMC 4.0 Keywords Manual Section 4 (*mech) 
MAC/GMC 4.0 Theory Manual Section 2.2.2 
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Section 3 : Analysis Type and Architecture 

Keywords 


This section deals with the keywords used to select the overall analysis type and the particular 
micromechanics model(s) that will be used by MAC/GMC 4.0 to analyze the particular problem. The two 
MAC/GMC 4.0 input file keywords covered in this section are: 

*RUC — » Invokes repeating unit cell (RUC) analysis, that is, models the behavior of a composite 
material (continuum) 

* LAMINATE — » Invokes laminate analysis, therefore models the behavior of a composite laminate 
(finite in one direction) with access to individual ply, RUC, and subcell data 

RUC analysis allows simulations to be performed on a heterogeneous (composite) material that is 
represented by a repeating unit cell. There are three separate micromechanics models within MAC/GMC 
4.0 that can be employed in RUC analyses: 1) Doubly periodic generalized method of cells (GMC), 2) 
Triply periodic GMC, and 3) Doubly periodic high-fidelity GMC (HFGMC). The doubly periodic GMC 
RUC is shown in Figure 3.1, while the triply periodic GMC RUC is shown in Figure 3.2. The RUC for 
doubly periodic HFGMC is identical to that of doubly-periodic GMC with the current restriction that the 
number of subcells in each direction must be even. By assigning the different constituent materials 
(specified under * CONSTITUENTS) to the subcells within the RUC, the RUC can be constructed to 
represent a composite material with the desired architecture. As the name indicates, the RUC repeats 
infinitely in either two (doubly periodic) or three (triply periodic) Cartesian coordinate directions. The 
doubly periodic RUC is infinitely long in the third (out-of-plane) x r direction (see Figure 3.1), so both 
RUCs represent an infinite heterogeneous (composite) medium with a periodic micro structure. Thus, 
GMC and HFGMC may be thought of as models for a material point that may be part of a larger overall 
structure. These micromechanics models not only calculate the effective properties of the composite 
material located at a particular material point, but also allow the simulation of applied loading (in the 
form of global stresses, global strains, a uniform temperature change, and electromagnetic loads) on the 
composite material (see Section 4 of this manual). 

Laminate analysis, on the other hand, simulates the behavior of a structure, namely a laminated plate. 
MAC/GMC 4.0 employs classical lamination theory (Jones, 1975; Herakovich, 1998) to model the 
composite laminate as infinite in the in-plane x- and y-directions and finite in the through-thickness z- 
direction (see Figure 3.3). Simulated loading for the laminate takes the form of in-plane stress and 
moment resultants and midplane strains and curvatures (see Section 4 of this manual). As Figure 3.3 
indicates, each ply of the laminate may be a composite material, which is then represented by a doubly or 
triply periodic GMC RUC. The GMC micromechanics models are thus embedded within the laminate 
theory analysis to represent the composite material points within the laminate plies. In the presence of 
laminate bending and inelastic material behavior, the stress and strain fields in the laminate will vary 
nonlinearly in the through-thickness (z-) direction. In such a case, MAC/GMC 4.0 tracks the fields in the 
laminate at two Gauss quadrature (integration) points in each layer (see Figure 3.3). During incremental 
application of the applied laminate-level simulated loading, the appropriate GMC micromechanics model 
is called at each of the two integration points in each layer to determine the local inelastic behavior of the 
composite material. The local inelastic behavior represented by the quadrature points can then be 
employed to determine the laminate level inelastic effects, which accounted for as inelastic force and 
moment resultants (see the MAC/GMC 4.0 Theory Manual Section 3). 
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I 1 I 2 1 3 U 

hNp 


h 

h 3 

h 2 

hi i 

Figure 3.1 Doubly periodic GMC and HFGMC repeating unit cell (RUC). The RUC is composed of 
subcells, each of which contains a constituent material. The subcells are identified by the 
indices Py. Note that for HFGMC, N and N must be even. 


P=Np 



z^Xl 



Figure 3.2 Triply periodic GMC repeating unit cell (RUC). The RUC is composed of subcells, each of 
which contains a constituent material. The subcells are identified by the indices a(3y. 
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z 



Figure 3.3 Geometry employed in MAC/GMC 4.0 for laminate theory analysis. Each ply is 
represented by a doubly or triply periodic GMC repeating unit cell (RUC), which operate at 
the Gauss quadrature points. 
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*RUC : Repeating unit cell analysis and architecture 

^ Either *RUC or * LAMINATE is required in all input files ^ 

Purpose: 

Invoke repeating unit cell analysis (as opposed to laminate analysis) for the problem, which requires the 
selection of the micromechanics model used to analyze the RUC and specification of the RUC 
architecture (i.e., arrangement of the composite constituents). 

The RUC section is broken into two sub-sections: 

I. Internal RUC Architecture Library (archidi- 9 9) 

II. User-Defined RUC Architecture {archid= 9 9) 


I. Internal RUC Architecture Library (archid± 99) 


Usage: 

*RUC 

MOD =modid ARCHID =archid —> specifics : see Table 3.1 and Table 3.2 <— 


Table 3.1 Specific input data for the internal doubly periodic RUC architecture library. 


Architecture 

Specific Data Input Format 

0 

Monolithic 

Material 

*RUC 

MOD=2 ARCHID=0 M =mm 

1 

Square Fiber, 
Square Pack 

*RUC 

MOD=2 ARCHID=1 VF =vf RAD =rad RITFR =ritfr F=mf M =mm I =mi 

2 

Square Fiber, 
Hexagonal Pack 

*RUC 

MOD=2 ARCHID=2 VF =vf RAD =rad RITFR =ritfr F =mf M =mm I =mi 

3 

Square Fiber, 
Square Diagonal 
Pack 

*RUC 

MOD=2 ARCHID=3 VF =vf RAD =rad RITFR =ritfr F=mf M=mm I =mi 

4 

Cross Shaped 
Fiber, Square 
Pack 

*RUC 

MOD=2 ARCHID=4 VF =1 >f XA =xa F=mf M =mm 

6 

7x7 Circular 
Fiber Approx., 
Rectangular 
Pack 

*RUC 

MOD=2 ARCH I D= 6 VF =vf R =R RITFR =ritfr F =mf M =mm I =mi 
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*RUC : Repeating unit cell analysis and architecture 


Table 3.1 (continued) 


9 1 

Two Different 
Size Fibers, 
Rectangular 
Pack 

*RUC 

MOD=2 ARCHID=9 VF1=v/7 R=R RADl=radl RITFRl=ritfrl & 
RAD2 =rad2 RITFR2 =ritfr2 F=mfl,mf2 M =mm l=mil,mi2 

11 

Square Fiber, 
Rectangular 
Pack 

*RUC 

MOD=2 ARCHID=11 VF =1 >f R =R RAD =rad RITFR =ritfr & 
F =mf M =mm I =mi 

13 

26x26 Circular 
Fiber Approx., 
Rectangular 
Pack 

*RUC 

MOD=2 ARCHID=13 VF =vf R =R F =mf M =mm 


Table 3.2 Specific input data for the internal triply periodic RUC architecture library. 


Architecture 

Specific Data Input Format 

0 

Monolithic 

Material 

*RUC 

MOD=3 ARCHID=0 M =mm 

1 

Short Fibers in 
Square Array 

*RUC 

MOD=3 ARCHID=1 VF =1 >f AS'9=asp F=mf M =mrtl 

2 

Short Fibers in 
Diagonal 
Array 

*RUC 

MOD=3 ARCHID=2 VF =vf ASP1 =aspl ASP2 =asp2 F =mf M =mm 

3 

Offset Short 
Fibers in 
Square Array 

*RUC 

MOD=3 ARCHID=3 VF =vf ASR 1-asp 1 ASP2 =asp2 DR =dr & 
F =mf M =mm 

4 

Ellipsoidal 

Inclusion, 

Rectangular 

Array 

*RUC 

MOD=3 ARCHID=4 OPT=l A =a B =b C =C D =d H =h LL =/ & 

F=mf M =mm 
-OR — 

MOD=3 ARCHID=4 OPT=2 VF =vf A =a B =b C=C RD =rd RL =rl & 
F=mf M =mm 
-OR- 

MOD=3 ARCHID=4 OPT=3 VF =vf RA =ra RC =rc D=< i H =h LL =/ & 
F =mf M =mm 
-OR- 

MOD=3 ARCHID=4 OPT=4 VF =vf A =a B =b C =C RD =rd & 
RLC=r/c F =mf M =mm 

10 2 

Open Cell 

*RUC 

MOD=3 ARCHID=10 VF =1 >f F=mf M =mm 
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Name 

modid 3,4 


archid 5 

specifics 
ritfr 6 

yf 7 6 

rad 6 
ritfr 6 
mf 

mm 

• 6 
mi 

xa 

R 

radl 
rad2 
ritfr 1 6 
ritfr2 6 
mfl,mf2 
mil,mi2 6 
asp, aspl 
asp2 
dr 

a,b,c 

d,h,l 

rd 

rl 

ra 

rc 

rlc 


Definition 

Periodic material model indentifier: 

= 1 — » monolithic material (internally uses doubly periodic GMC with one subcell) 

= 2 — » doubly periodic GMC model 

= 3 — » triply periodic GMC model 

= 12 — > doubly periodic HFGMC model 
I.D. number of RUC architecture (see Figure 3.4, Figure 3.5, and Figure 3.6) 

(archid= 9 9 for user-defined) 

This format depends on the modid and archid selected (see Table 3.1 and Table 3.2) 

Ratio of interface thickness to fiber radius (Special Case) 

Fiber (or inclusion) volume fraction 
Fiber radius (units arbitrary) (Special Case) 

Ratio of interface thickness to fiber radius (Special Case) 

Fiber material number (one of matnum from *C0NSTITUENTS) 

Matrix material number (one of matnum from CONSTITUENTS) 

Interface material number (one of matnum from CONSTITUENTS) (Special Case) 

Length of arm of cross shaped fiber (see Figure 3.4) (assumes fiber diameter of 1.) 

Aspect ratio of the RUC = L/H (see Figure 3.1 and Figure 3.4) 

Radius of fiber number 1 (smaller fiber) for archid= 9 (units arbitrary) 

Radius of fiber number 2 (larger fiber) for archid= 9 (units arbitrary) 

Ratio of interface number 1 thickness to fiber number 1 radius for archid= 9 (Special Case) 
Ratio of interface number 2 thickness to fiber number 2 radius for archid= 9 (Special Case) 
Fiber material numbers for the two fibers in archid= 9 
Interface material numbers for the two fibers in archid= 9 (Special Case) 

Short fiber aspect ratio (fiber length / fiber diameter) 

RUC aspect ratio (d/h) (see Figure 3.2 and Figure 3.6) 

Ratio dfd 3 in archid= 3, quantifies fiber offset (see Figure 3.6) 

Ellipsoidal inclusion semi-major axes in archid= 4 

Overall dimensions of the RUC (see Figure 3.2 and Figure 3.6) 

RUC dimension aspect ratio d/h (see Figure 3.2 and Figure 3.6) 

RUC dimension aspect ratio l/h (see Figure 3.2 and Figure 3.6) 

Semi -major axis aspect ratio a/b (see Figure 3.6) 

Semi-major axis aspect ratio c/b (see Figure 3.6) 

Aspect ratio //(2c) (see Figure 3.6) 


Examples: 

*RUC 

MOD=l M=2 

*RUC 

MOD=2 ARCHID=1 VF=0 . 65 F=1 M=2 

*RUC 

MOD=3 ARCHID=1 VF=0.25 ASP=1 F=1 M=2 

*RUC 

MOD=2 ARCH I D= 6 R=1 . 0 VF=0.25 RITFR=0 . 1 F=1 M=2 1=3 
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*ruc : Repeating unit cell analysis and architecture 


*RUC 

MOD=2 ARCH I D= 9 VF1=0.10 RAD1=50. RITFR1=0 . 1 & 

VF2=0 . 15 RAD2=71 . RITFR2=0 . 1 R=1 . F=4,l M=2 1=3,3 


*RUC 

MOD=3 ARCHID=0 M=2 


Associated MAC/GMC 4.0 Example Problemfs): 

Example la-c, Example Id, Example 2a-h, Example 3a, Example 3b, Example 3d, Example 3f, 
Example 4a-c, Example 4f, Example 5a-b, Example 5d-f, Example 6a-b, Example 7a-b 


Notes: 

1) For archid= 9, the fiber volume fraction of fiber 2 is adjusted to preserve the correct packing 
arrangement and RUC aspect ratio. 

2) The open cell triply periodic RUC (archid= 10) places the matrix material (mm) in the four colored 
subcells shown in Figure 3.6 (i.e., the “arms” are the matrix material). The “open” material shown in 
the figure is occupied by the fiber material. Thus, to simulated a truly “open” material, the fiber 
elastic stiffness should be set very low so as to provide little contribution to the overall RUC 
response. 

3) modid= 1 is a special case used in MAC/GMC 4.0 to select a monolithic material. Internally, the code 
employs a double periodic GMC RUC analysis with one subcell to which the material is assigned 
using n=mm (see Example above). A monolithic material may also be selected by using archid= 0 in 
either doubly or triply periodic GMC. 

4) Doubly periodic HFGMC RUC analysis requires an even number of subcells in both the x 2 - and x 3 - 
directions and additional specialized input under the * SOLVER keyword. 

5) The RUC architectures contained within the MAC/GMC 4.0 architecture library are shown in Figure 
3.4 (doubly periodic, no interface), Figure 3.5 (doubly periodic, with interface), and Figure 3.6 (triply 
periodic). 

6) The ritfr , mi, rad , ritfrl , ritfr2 , mil , and mi2 data is associated with the special case involving an 
RUC with a separate interface material (see Figure 3.5). Only those RUCs shown in Figure 3.5 have 
the interface option available. To specify one of these RUCs with an interface, simply include the 
special case data shown in blue. 

7) Due to geometric constraints (i.e., adjacent fibers, inclusions, or interfaces touching) each 
MAC/GMC 4.0 internal RUC architecture has some restraints on its fiber volume fraction ( vj ). These 
constraints are given for each architecture in Table 3.3. 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 3 (*laminate) 
MAC/GMC 4.0 Keywords Manual Section 3 (*SOLVER) 
MAC/GMC 4.0 Theory Manual Section 2 
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Figure 3.4 MAC/GMC 4.0 internal library of doubly periodic repeating unit cells (RUCs) without 
interface. 
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Section 3: Analysis Type and Architecture Keywords 


*RUC : Repeating unit cell analysis and architecture 





Matrix Fiber or Fiber 2 Interface or Interface2 
Fiber 1 Interfacel 


x 2 



-^x 3 

ARCHID = 1 


x 2 



— x 3 


ARCHID = 2 


x 2 


- x 3 


ARCHID = 3 



I 


X 2 


1 


I 


i 


— x 3 


ARCHID = 9 


x 2 



^3 


ARCHID = 11 


Figure 3.5 MAC/GMC 4.0 internal library of doubly periodic repeating unit cells (RUCs) with 
interface. 
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ARCHID = 4 


ARCHID = 10 


Figure 3.6 MAC/GMC 4.0 internal library of triply periodic repeating unit cells (RUCs). 
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Table 3.3 Constraints on the fiber volume fraction of the MAC/GMC 4.0 internal RUC architectures. 


Doubly Periodic RUC Architectures 

Architecture 

Limit(s) on vf 

0 

Monolithic Material 


1 

Square Fiber, 
Square Pack 

CN 

^ ’ik 
+ 

VI 

V 

2 

Square Fiber, 
Hexagonal Pack 

, V3 

vf < , v> 

2 (l + ritfr) 

3 

Square Fiber, 
Square Diagonal Pack 

. f . 0.5 

(l + ritfr) 2 

4 

Cross Shaped Fiber, 
Square Pack 

vf< l-4(xa ) 2 

6 

7x7 Circular Fiber Approx., 
Rectangular Pack 

vf <— ^ R> 1 

7 ?( \ + ritfr) 

0.8125 R n , 

v/< R< 1 

( \ + ritfr) 

7 

14x14 Circular Fiber 
Approx., Rectangular Pack 

0.8148 n , 

vf < R > 1 

R 

vf< 0.81487? R< 1 

9 

Two Different Size Fibers, 
Rectangular Pack 

v f r > ^ ^ 

yj ^ ^ r ~|2 

R (l + ritfr2)+ (l + ritfr ] ) rac ^ 

rad 2] 

2jR(vf2) 

J / \2 

( rad 2 \ 

(1 + ritfr 2)+ 

y FGG 1 y 

v /1 ( rad\ Y 
1 / 2 \rad2 y 

radl\. . . 

(l + ritfrl) 

FGU 1 y 

2 

>1 

11 

Square Fiber, Rectangular 
Pack 

vf < 7 rr for 7? > 1 . 

7?(l + ritfr) 

vf < 77 for 7? < 1 . 

(l + ritfr) 2 

13 

26x26 Circular Fiber 
Approx., Rectangular Pack 

, 0.80613 r n , 

vf < for 7? > 1 . 

7? 

yf< 0.806137? for 7?<1. 
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Table 3.3 (continued) 


Triply Periodic RUC Architectures 

0 

Monolithic Material 

- 

1 

Short Fibers in Square 
Array 

v/< 1. 

2 

Short Fibers in Diagonal 
Array 

, 1 

vf < — 
2 

, 1 

vf <- 

2 

OSpl 
: ^osp2 

^ aspl' 

y aSpl j 

\ 

for asp2 > aspl 

2 

for aspl > asp 2 

3 

Offset Short Fibers in 
Square Array 

f 

vf< 4 

V 

, 1 

vf < 

2 

asp] 

aspl ^ 
asp 2 y 
^ aspl' 

y aSpl j 

, 1 

vf >~ : 

\ > aspl 

for aSpl >2 

aspl 

for ^‘<2 

aspl 

. r aspl\ 

\ ^ aspl J 

4 

Ellipsoidal Inclusion, 
Rectangular Array 

Analytical expressions unavailable 

10 

Open Cell 3 

v/< 1. 


II. User-Defined Architecture ( archid=99 ) 


Usage: 

*RUC 

MOD =modid ARCH ID =archid 
NA =na NB =nb N G=ng 
D=d(l),d(2), ...,d(na) 

R=h(l),h(2), ...,h(nb) 
L=l(l),l(2),..J(ng) 

SM=... — » specifics: see Table 3.4 <- 
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*RUC : Repeating unit cell analysis and architecture 


Table 3.4 Format of the subcell material specification for user-defined RUC architectures. 


Doubly Periodic RUC Analysis 

Triply Periodic RUC Analysis 

SM =m(nb,l),m(nb,2), ...,m(nb,ng) 
SM=m(nb-l ,l),m(nb-l ,2), ...,m(nb-l,ng) 

# -- gamma = 1 

SM=m(na,l,l),m(na,2,l), ...,m(na,nb,l) 
SM=m(na-l,l,l),m(na-l,2,l), ...,m(na-l,nb,l) 

SM=m(l,l),m(l,2), ...,m(l,ng) 
(See Figure 3.7) 

SM=m(l,l,l),m(l,2,l), ...,m(l,nb,l) 

# -- gamma = 2 

SM=m(na, l,2),m(na,2,2), ...,m(na,nb,2) 
SM=m(na-l,l,2),m(na-l,2,2), ...,m(na-l,nb,2) 


SM=m(1,1,2),m( 1,2,2), ...,m(l,nb,2) 


# -- gamma = ng 

SM=m(na,l, ng), m(na,2, ng), ...,m(na,nb, ng) 
SM=m(na-l , 1 ,ng),m(na-l ,2,ng), ...,m(na-l,nb,ng) 


SM=m(l , 1 ,ng),m(l ,2,ng), ...,m(l,nb,ng) 


(See Figure 3.8) 


Name Definition 

modid 1,2 Periodic material model indentifier: 

— 1 —> monolithic material (internally uses doubly periodic GMC with one subcell) 
= 2 — » doubly periodic GMC model 
= 3 — » triply periodic GMC model 
= 12 — » doubly periodic HFGMC model 
archid I.D. number of RUC architecture (archid= 9 9 for user-defined) 

[=99] 

na 3 Number of subcells in the Xi-direction (Special Case) (see Figure 3.2) 

nb Number of subcells in the x 2 -direction (see Figure 3.1 and Figure 3.2) 

ng Number of subcells in the x 3 -direction (see Figure 3.1 and Figure 3.2) 

d(l-na) 3 Subcell x r direction dimensions (Special Case) (see Figure 3.2) 
h(l-nb) Subcell x 2 -direction dimensions (see Figure 3.1 and Figure 3.2) 

1(1 -ng) Subcell x 3 -direction dimensions (see Figure 3.1 and Figure 3.2) 

specifics This format depends on whether the RUC is doubly or triply periodic (see Table 3.4) 
m(ib,ig) Constituent material assigned to subcell p =ib y =ig (doubly -periodic case) 
m(ia,ib,ig) Constituent material assigned to subcell a=ig P =ib y =ig (triply-periodic case) 
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Examples: 

Doubly Periodic 

*RUC 

MOD=2 ARCHID=99 
NB=3 NG=3 
H=0 . 5, 0 . 5, 1 . 0 
L=0 . 5, 1 . 0, 0 . 5 
SM=2 ,2,2 
SM=2 ,2,3 
SM=1 ,2,2 


7Wp/y Periodic 

*RUC 

MOD=3 ARCHID=99 
NA=4 NB=2 NG=2 
D=1 . , 1 . , 1 . , 1 . 
H=1 . , 1 . 

L=1 . , 1 . 

# -- gamma = 1 
SM=2 , 2 

SM=2 , 2 
SM=2 , 2 
SM=1, 2 

# -- gamma = 2 
SM=2 , 2 

SM=1 , 1 
SM=2 , 2 
SM=1, 2 


x 2 

(P) 


h, 

3,1 

3,2 

3,3 

h 2 

2,1 

2,2 

2,3 

hi 

1,1 

1,2 

1,3 


i. 

i 2 

Y dvl 

h 



matnum=1 


matnum=2 


matnum=3 


Figure 3.7 Example user-defined doubly periodic RUC. 




Figure 3.8 Example user-defined Triply periodic RUC. 


Associated MAC/GMC 4.0 Example Problem (s): 

Example 3c, Example 3e, Example 5f, Example 7d, Example 7e 
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*ruc : Repeating unit cell analysis and architecture 


Notes: 

1) modid= 1 is a special case used in MAC/GMC 4.0 to select a monolithic material. Internally, the code 
employs a double periodic GMC RUC analysis with one subcell to which the material is assigned 
using M =mm (see Example above). A monolithic material may also be selected by using archid= 0 in 
either doubly or triply periodic GMC. 

2) Doubly periodic HFGMC RUC analysis requires an even number of subcells in both the x 2 - and x 3 - 
directions and additional specialized input under the * SOLVER keyword. 

3) NA =na and D=d(l),d(2)...,d(na) are only required in the special case of triply periodic RUC analysis 


Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 3 (*laminate) 
MAC/GMC 4.0 Theory Manual Section 2 
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* LAMINATE : Laminate analysis and architecture 


^ Either *RUC or * LAMINATE is required in all input files ^ 

Purpose: 

Invoke laminate analysis (as opposed to repeating unit cell analysis) for the problem, which requires 
selection of a micromechanics model as well as specification of an RUC architecture (i.e., arrangement of 
the composite constituents) for the analysis of each each ply. Lastly, the angle and thickness of each ply 
must be specified. 

The * LAMINATE section is broken into two sub-sections: 

I. Layer RUC Architecture from Internal Library (archid^ 9 9) 

II. User-Defined Layer RUC Architecture (archid= 9 9) 


I. Layer RUC Architecture from Internal Library (archid£99) 


Usage: 

* LAMINATE 

NLY =nly 

LY =ly MOD =modid THK =thk AN G=ang ARCHID =archid & 
—> specifics : see Table 3.1 and Table 3.2 <— EM =em 


Repeated nly times 


Name Definition 


nly 

ly 

modid 


1,2 


thk 
ang 
arc hid 
specifics 
em 5 


3 


Number of layers (plies) comprising the laminate 

Layer number (data following LY=/y is associated with layer number ly) 

Periodic material model indentifier: 

= 1 — » monolithic material (internally uses doubly periodic GMC with one subcell) 
= 2 —> doubly periodic GMC model 
= 3 —> triply periodic GMC model 
= 12 — > N/A (HFGMC unavailable within laminate) 

Thickness of layer 

Off-axis angle of layer (see Figure 3.3) 

I.D. number of RUC architecture (archid= 9 9 for user-defined) 

This format depends on the modid and archid selected (see Table 3.1 and Table 3.2) 
Electromagnetic specifier for layer number ly: 

= 1 — » Layer is electromagnetic 
= 0 — » Layer is not electromagnetic 
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* laminate : Laminate analysis and architecture 


Examples: 


♦LAMINATE 

NLY=3 

LY=1 MOD=2 
LY=2 MOD=2 
LY=3 MOD=2 


THK=0 . 25 ANG=90 ARCHID=1 VF=0 . 65 F=1 M=2 
THK=0 . 50 ANG=0 ARCHID=1 VF=0 . 65 F=1 M=2 
THK=0 . 25 ANG=90 ARCHID=1 VF=0 . 65 F=1 M=2 


♦LAMINATE 

NLY=4 

LY=1 MOD=l THK=0 . 25 M=2 

LY=2 MOD=2 THK=0 . 25 ANG=60 ARCHID=6 R=1 . VF=0.25 F=1 M=2 
LY=3 MOD=3 THK=0 . 25 ANG=0 ARCHID=3 ASP1=3. ASP2=2 . DR=1 . & 
VF=0 . 25 F=1 M=2 

LY=4 MOD=3 THK=0 . 25 ANG=0 . ARCHID=1 VF=0.25 ASP=1 F=1 M=2 


Associated MAC/GMC 4.0 Example ProblemfsV 

Example lb, Example 3g, Example 3h, Example 3i, Example 4d, Example 4e, Example 4g, Example 
5c, Example 5d, Example 6c 


Notes: 

1) modid= 1 is a special case used in MAC/GMC 4.0 to select a monolithic layer. Internally, the code 
employs a double periodic GMC RUC analysis for the layer with one subcell to which the material is 
assigned using M =mm (see Example above). A monolithic material may also be selected for a layer 
by using archid= 0 in either doubly or triply periodic GMC. 

2) The ability to employ triply periodic layers in a laminate is new to version 4.0 of the code. Doubly 
periodic HFGMC RUC analysis is not currently available for laminate layers. 

3) Refer to the *RUC description earlier in this section of this manual for the specific input requirements 
for each RUC architecture. The RUC architectures contained within the MAC/GMC 4.0 architecture 
library are shown in Figure 3.4 (doubly periodic, no interface), Figure 3.5 (doubly periodic, with 
interface), and Figure 3.6 (triply periodic). 

4) The ritfr and mi data is associated with the special case involving a layer with an RUC that has a 
separate interface material (see Figure 3.5). Only those RUCs shown in Figure 3.5 have the interface 
option available. To specify one of these RUCs with an interface for the layer, simply include the 
special case data shown in blue. 

5) The electromagnetic specifier (EM =em) for the layer tells the code whether or not a particular layer is 
electromagnetically active. This is only required in the special case in which an electromagnetic 
analysis has been selected by including the *ELECROMAG keyword in the input file. 

Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 3 (*RUC) 

MAC/GMC 4.0 Theory Manual Section 3 
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II. User-Defined Layer RUC Architecture {archid— 99) 


Usage: 

* LAMINATE 

NLY =nly 

LY =ly MOD =modid THK =thk AN G=ang ARCHID =archid EM =em & 
MOD =modid ARCHID =archid 
NA =na NB =nb NG=ng 
D=d(l),d(2), ...,d(na) 

R=h(l),h(2), ...,h(nb) 

L=l(l),l(2),..J(ng) 

SM=... — > specifics: see Table 3.4 <- 


Repeated nly times 


Name 

nly 

ly 

modid 


thk 
ang 
arc hid 
[=99] 
cm 3 


4 

na 

nb 

ng 

d(l-na) 4 

h(l-nb) 

l(l-ng) 

specifics 

m(ib,ig) 

m(ia,ib,ig) 


Definition 

Number of layers (plies) comprising the laminate 

Layer number (data following LY =ly is associated with layer number ly) 

Periodic material model indentifier: 

= 1 — » monolithic material (internally uses doubly periodic GMC with one subcell) 
= 2 — » doubly periodic GMC model 
= 3 —> triply periodic GMC model 
= 12— > N/A (HFGMC unavailable within laminate) 

Thickness of layer 

Off-axis angle (Figure 3.3) of layer 

I.D. number of RUC architecture (archid= 9 9 for user-defined) 

Electromagnetic specifier for layer number ly (Special Case): 

= 1 — » Layer is electromagnetic 
= 0 — » Layer is not electromagnetic 
Number of subcells in the x r direction (Special Case) (see Figure 3.2) 

Number of subcells in the x 2 -direction (see Figure 3.1 and Figure 3.2) 

Number of subcells in the x 3 -direction (see Figure 3.1 and Figure 3.2) 

Subcell x r direction dimensions (Special Case) (see Figure 3.2) 

Subcell x 2 -direction dimensions (see Figure 3.1 and Figure 3.2) 

Subcell x 3 -direction dimensions (see Figure 3.1 and Figure 3.2) 

This format depends on whether the RUC is doubly or triply periodic (see Table 3.4) 
Constituent material assigned to subcell (3 =ib y =ig (doubly -periodic case) 

Constituent material assigned to subcell a=ig $=ib y =ig (triply-periodic case) 
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* laminate : Laminate analysis and architecture 


Examples: 

♦LAMINATE 

NLY=2 

LY=1 THK=0 . 5 ANG=0 . MOD=2 ARCHID=99 
NB=1 NG=2 
H=1 . 

L=0 .5,0.5 
SM=1 , 2 

LY=2 THK=0 . 5 ANG=90. MOD=2 ARCHID=1 VF=0.25 F=1 M=2 


♦LAMINATE 

NLY=3 

LY=1 THK=0 . 25 ANG=0 . MOD=3 ARCHID=99 EM=1 
NA=1 NB=2 NG=2 
D=1 . 

H=0 .5,0.5 
L=0 .5,0.5 
SM=1, 2 
SM=2 , 2 

LY=2 THK=0 . 5 ANG=90. MOD=2 ARCHID=1 VF=0.25 F=4 M=5 EM=0 
LY=3 THK=0 . 25 ANG=0 . MOD=3 ARCHID=99 EM=1 
NA=1 NB=2 NG=2 
D=1 . 

H=0 .5,0.5 
L=0 .5,0.5 
SM=1 , 2 
SM=2 , 2 


Associated MAC/GMC 4.0 Example Problem^: 

Example 7c 

Notes: 

1) modid= 1 is a special case used in MAC/GMC 4.0 to select a monolithic layer. Internally, the code 
employs a double periodic GMC RUC analysis for the layer with one subcell to which the material is 
assigned using M =mm (see Example above). A monolithic material may also be selected for a layer 
by using archid= 0 in either doubly or triply periodic GMC. 

2) Doubly periodic HFGMC RUC analysis is not currently available for laminate layers. 

3) The electromagnetic specifier (EM =em) for the layer tells the code whether or not a particular layer is 
electromagnetically active. This is only required in the special case in which an electromagnetic 
analysis has been selected by including the *ELECROMAG keyword. 

4) modid= NA and D=d(l),d(2)...,d(na) are only required in the special case of triply periodic RUC 
analysis. 

Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 3 (*RUC) 
MAC/GMC 4.0 Theory Manual Section 3 
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Section 4 : Loading Keywords 


This section presents the four MAC/GMC 4.0 keywords that are related to the applied loading. These are: 
*MECH — » Specify mechanical (and electromagnetic) loading profile 

* THERM — » Specify thermal loading profile 

* SOLVER — » Specify time -integration data 

*SURF — » Specify loading (and additional information) for yield surface analysis 

Because MAC/GMC 4.0 allows for inelastic material behavior, the simulated loading must be applied in 
an incremental manner. The problem solutions are (in general) path-dependent, so the overall solution (at 
the completion of the applied loading) actually involves solving the problem at each increment of the 
simulated applied loading. As such, a simulated loading profile must be specified rather than just the 
final load levels. The loading profile is essentially a load vs. time history that MAC/GMC 4.0 will apply 
to the specified RUC or laminate. 

For repeating unit cell analysis, mechanical loading takes the form of an admissible time profile for each 
of the six global (stress or strain) components present in the RUC constitutive equation. For example, 
consider the standard (non-electromagnetic) repeating unit cell effective constitutive equation, 


"^11 " 


*11 — *n -*n 

^22 


8 — 8 P — £ T 

O 22 ty 22 ty 22 

^22 

= [c] 

8 — 8 P — £ T 

c 22 ^33 ^33 

— —p —T 

^23 


y 22 ~ y 22 ~ y 22 

^13 


^1 

UJ 

1 

^1 

1 

^1 

_ G \2 _ 


_/i2 y \2 y \2 


where cr. are the global (or average) stress components on the composite, C* is the global (or effective) 
stiffness matrix, 8 tj are the global (or average) normal strain components on the composite, y tj are the 
global (or average) engineering shear strain components on the composite, e?. are the global (or average) 
normal inelastic strain components on the composite, y? are the global (or average) engineering shear 
inelastic strain components on the composite, sfj are the global (or average) normal thermal strain 

components on the composite, and yjj are the global (or average) engineering shear thermal strain 

components on the composite. Admissibility simply requires that, if a particular strain component is 
specified, the corresponding stress component cannot be specified (and vice versa). In addition to the 
arbitrary mechanical loading option (which allows the user to specify all six loading components), 
MAC/GMC 4.0 also includes several internal loading options that simplify the application of loading for 
commonly used mechanical loading types. In the special case of electromagnetic RUC analysis, the 
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number of loading components is increased from six to twelve to accommodate applied electric and 
magnetic load components (see the MAC/GMC 4.0 Theory Manual Section 2.1.2). 

The applied mechanical loading for laminate analysis is different than that of RUC analysis. Globally, 
rather than dealing with average stresses and strains, the applicable quantities are force and moment 
resultants along with midplane strains and curvatures (see the MAC/GMC 4.0 Theory Manual Section 3). 
For laminate analysis, the code’s internal mechanical loading options are numbered differently than those 
associated with RUC mechanical loading. As is the case in electromagnetic RUC analysis, 
electromagnetic laminate analysis involves additional electromagnetic loading components. 

MAC/GMC 4.0 also allows application of a simulated time vs. temperature profile for both RUC and 
laminate analysis using the * THERM keyword. It should be noted that the temperature at every point in 
time is treated as constant throughout the RUC or laminate; the code does not solve any thermal boundary 
value problems, but rather treats the temperature history as known. The * SOLVER keyword is used to 
select either a constant time step forward Euler time integration scheme, or a self-adaptive predictor- 
corrector time integration scheme. Finally, the yield surface option keyword (*SURF) is a special case in 
that the yield surface “probing” (i.e., loading along different angles in a particular stress space) controls 
the applied loading. In this case, the data specified under *MECH, * THERM, and *SOLVER are considered 
to be a “preload” that occurs before the start of the yield surface analysis. 
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*MECH : Applied mechanical loading 


Purpose: 

Specify the simulated mechanical (and electromagnetic) loading profile for the RUC or laminate. 


Usage: 

*MECH 

LOP =lop REFTIME =reftime 

NPT =npt Tl=t(l),t(2),...,t(npt) MAG=Mag(l),rnag(2),...,mag(npt) & 
MODE=mode( l ), mode(2), . ..,mode(npt-l) 


— » Repeated number of times equal to number of specified load components (see Table 4.1) 


Name Definition 

lop 1 Loading option (see Table 4.1) (lop = 99 —> arbitrary mechanical loading option) 

reftime 2 Strain reference time (Optional) 

npt Number of time vs. load magnitude points in specified profile 

t(i) 3 Time points in time vs. load magnitude profile (see Figure 4.1) 

mag(i) 4,5 Load magnitude points in time vs. load magnitude profile (see Figure 4.1) 

mode(j) 6 Mode of applied loading during loading step from t(j) to t(j+l) 

= 1 — » strain, midplane strain, electric field, magnetic field applied 
= 2 — » stress, force resultant, electric displacement, magnetic flux density applied 


Figure 4.1 
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*mech : Applied mechanical loading 


Examples: 



Time (seconds) 


*MECH 

L0P=1 

NPT=3 TI=0 . , 200 . , 300 . MAG=0 . , 0 . 02, 0 . 01 MODE=l,l 



*MECH 

LOP=7 

NPT=4 TI=0. , 100. , 200. , 300. MAG=0 . , 50 . , 50 . , 0 . MODE=2,2 
NPT=3 TI=0. , 150. , 300. MAG=0 . , -25 . , 25 . MODE=2,2 


*MECH 

LOP=99 

NPT=6 TI=0 ., 300 ., 310 ., 2310 2320 2620 . MAG=0 . , 0 . , 0 . 01, 0 . 01, 0 . , 0 . & 
MODE=l, 1, 1, 2, 2 

NPT=6 TI=0 300 310 2310 2320 2620 . MAG=0 . , 0 . , 0 . 005, 0 . 005, 0 . , 0 . & 
MODE=l, 1, 1, 2, 2 

NPT=2 TI=0 . , 2620 . MAG=0 . , 0 . MODE=2 
NPT=2 TI=0 . , 2620 . MAG=0 . , 0 . MODE=2 
NPT=2 TI=0 . , 2620 . MAG=0 . , 0 . MODE=2 
NPT=6 TI=0 ., 300 ., 310 ., 2310 ., 2320 ., 2620 . & 

MAG=0 . ,0., 0.0075, 0.0075, 0 . , 0 . MODE=l, 1, 1, 2, 2 


*MECH 

LOP=99 

NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 


MAG=0 . , 0 . 

MODE=2 

MAG=0 . , 0 . 

MODE=2 

MAG=0 . , 0 . 

MODE=2 

MAG=0 . , 0 . 

MODE=2 
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NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 
NPT=2 TI=0 . , 600 . 


MAG=0 . , 0 . MODE=2 
MAG=0 . , 0 . MODE=2 
MAG=0 . , 0 . MODE=l 
MAG=0 . , 0 . MODE=l 
MAG=0 . , 0 . MODE=l 
MAG=0 . , 0 . MODE=l 
MAG=0 . , 0 . MODE=l 
MAG=0 . , 6.E6 MODE=l 


Table 4.1 MAC/GM 

C 4.0 loading options for c 

oubly and triply periodic RUC and laminate analysis. 

lop 

number 

components 

specified 

Doubly Periodic RUC 
Analysis 

Triply Periodic RUC 
Analysis 

Laminate Analysis 



mode = 1 

mode = 2 

mode = 1 

mode = 2 

mode = 1 

mode = 2 

1 

1 

en 

an 

£11 

On 

£ ° 

^xx 

N xx 

2 

1 

£22 

022 

£22 

O22 

CO 

^ 0 

Nyy 

3 

1 

£33 

O33 

£33 

033 

y° 

/ xy 

N X y 

4 

1 

Y23 

023 

Y23 

O23 

Kxx 

M xx 

5 

1 

Yl 3 

Oi 3 

Yl 3 

O13 

K yy 

Myy 

6 

1 

Y 13 

Oi 2 

Y 13 

O12 

*xy 

M X y 

7 

2 

8 ll & 822 

CTn & O22 

£ll & £22 

Cn & O22 

N/A 

8 

2 

822 & £33 

022 & O33 

822 & £33 

G22 & O33 

N/A 

9 

2 

£ll & Y 23 

On & O23 

£ll & £33 

On & O33 

N/A 

10 

2 

£22 & Yl 3 

a22 & O13 

£ll & Y 23 

On & O23 

N/A 

11 

2 

N/A 

£22 & Yl3 

022 & O13 

N/A 

12 

2 

N/A 

£33 & Y 12 

G33 & G12 

N/A 

99 

6 

6 of 8ij or ay 

6 of 8y or Gy 

3 of s y or Ny, 
3 of Ky or My 

12 

(electromag) 

6 of 8ij or ay, 
3 of E k or D k , 
3 of H k or B k 

6 of 8ij or Gy, 
3 of E k or D k , 
3 of H k or B k 

3 of or Ny, 

3 of k ,, or Mi,, 
3 ofE k , 3 of H k 


Note: a M = stress, = normal strain, y M = engineering shear strain, N u = force resultant, My = moment 
resultant, sjj = midplane strains, Ky = midplane curvatures, D k = electrical displacements, E k = 
electric fields, B k = magnetic flux densities, H k = magnetic fields. 


Associated MAC/GMC 4.0 Example Problem fs): 

Example lc-d, Example 2a-h, Example 3a-i, Example 4a, Example 4b, Example 4c, Example 4d, 
Example 4e, Example 5b-f, Example 6a-b, Example 7a, Example 7b, Example 7c 
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*mech : Applied mechanical loading 


Notes: 

1) As shown in Table 4.1, certain values of lop are invalid in certain cases. In particular, lop= 11 and 
lop= 12 are only valid for triply periodic RUC analysis, while for laminate analysis, the only valid 
values for lop are 1-6 and 9 9 . 

2) The optional reftime =reftime specifier is used to set a reference strain value other than zero. The 
reference strain is the strain value that all the rest of the specified loading strain values are taken with 
respect to. For example, if the reference strain were 0.015, and a strain load of 0.01 were specified, 
the code would increase the strain to 0.025. Instead of specifying the reference strain directly, a time 
is given at which the strain state is taken as the six reference strain components. This feature is useful 
when simulating residual stresses prior to applying loading in the form of strains (see MAC/GMC 4.0 
Example Problem 4b). 

3) The first time specified for each mechanical (and electromagnetic) loading component must be zero, 
and the final time for each component must be the same. However, the number of intermediate times, 
as well as the intermediate times themselves, can be different for each applied loading component. 

4) The units of the applied loading should be consistent with the units employed for the constituent 
materials. That is, for example, if the units of the constituent material elastic moduli are MPa, any 
applied stress components will be interpreted by the code as having the units MPa. 

5) The first specified load magnitude should be zero. The code must begin its execution in an unloaded 
state. 

6) In the special case of an electromagnetically active laminate, the electromagnetic loading currently 
can only take the form of applied electric field components (E k ) and magnetic field (H k ) components. 
That is, laminate electric displacement ( D k ) and magnetic flux density ( B k ) components cannot 
currently be applied. 


Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 1 (*ELECTROMAG) 
MAC/GMC 4.0 Theory Manual Sections 2 and 3 
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* THERM : Applied thermal loading 


Purpose: 

Specify the simulated thermal loading profile (time vs. temperature) for the RUC or laminate. 


Usage: 

* THERM 

NPT =npt Tl=t(l),t(2),...,t(npt) TEMP =T(l),T(2),...,T(npt) 

Name Definition 

npt Number of time vs. temperature points in specified profile 

t(i) 1 Time points in time vs. temperature profile (see Figure 4.2) 

T(i) 2 Temperature points in time vs. load temperature profile (see Figure 4.2) 



Figure 4.2 Example temperature vs. time profile for applied thermal loading. 
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*therm: Applied thermal loading 


Example: 



Time (seconds) 


*THERM 

NPT=3 TI=0 6000 6600 . TEMP=900 . , 23 . , 650 . 


Associated MAC/GMC 4.0 Example Problemfsl: 

Example Id, Example 2e-h, Example 3a-i, Example 4a, Example 4b, Example 4c, Example 4d, 
Example 4e, Example 5b-f, Example 6a-b, Example 7a-c 


Note: 

1) The starting time for the time vs. temperature profile must be zero. The ending time for the time vs. 
temperature profile must correspond to the ending time for any specified mechanical loading. 

2) The temperature unit in MAC/GMC 4.0 is °C. The user may employ other temperature units if none 
of the constituent materials are taken from the code’s internal material database. 


Additional Information: 


MAC/GMC 4.0 Theory Manual Sections 2 and 3 
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* SOLVER : Time integration 

Purpose: 

Allows specification of time integration scheme and associated integration related data. 

Only two choices are currently available, thus the * SOLVER section is broken into two corresponding sub- 
sections: 

I. Forward Euler Time Integration (specified time step) (method= l) 

II. Predictor-Corrector Time Integration (self-adaptive time step) (method= 2) 


I. Forward Euler Time Integration (method= 1) 


Usage: 

* SOLVER 

METHOD =method NPT =npt Tl=t(l),t(2), ...,t(npt) STP =stp(l),stp(2),...,stp(npt-l) & 

ITMA X=itmax ERR =errtol 
NLE G=nleg NINTE G=nint 

Name Definition 

method Time integration method: 

[=1] = 1 — » Forward Euler (specified time step) 

= 2 — » Predictor-Corrector (self-adaptive time-step) 
npt Number of time vs. step size points in specified profile 
t(i) 1 Time points in time vs. step size profile 
stp(i) 2 Step sizes in time vs. step size profile 

Umax 3 Maximum number of global equilibrium iterations (Special Case) 
errtol 3 Error tolerance for convergence of global iterations (Special Case) 

nleg 4,5 Order of the Legendre polynomial used to approximate the inelastic strain field (Special Case) 

nint 4,5 Number of integration points used to represent the field variables (Special Case) 


Examples: 



stp(1) = 0.1 

stp(2) = 0.5 

stp(3) = 0.2 


to: 

) = 0 t(2) 

i = 60. 


t(3 

) = 180. t(4; 

) = 240. 


* SOLVER 

METHOD=l NPT=4 TI = 0 . , 60 . , 1 8 0 . , 2 4 0 . STP=0 . 1 , 0 . 5 , 0 . 2 


time (s) 


* SOLVER 

METHOD=l NPT=2 TI=0.,20. STP=0 . 1 ITMAX=20 ERR=0.0001 
NLEG=5 NINTEG=11 
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★solver: Time integration 


Associated MAC/GMC 4.0 Example Problem(s): 

Example lc-d, Example 2a-h, Example 3a-i, Example 4a-f, Example 5b-f, Example 6a-c, Example 7a-c 


Note: 

1) The starting time in the time vs. step size profile must be zero. The ending time must correspond to 
the ending time of any mechanical and thermal loading specified under *MECH and *THERM. 

2) The forward Euler integration scheme involves stepping through time with a specified time step size . 
At each time step, increments of the field variables throughout the RUC or laminate are calculated 
and simply added to the previous field variable values. The user specifies a number of time points 
along with the time step size to be used between each successive time point (see Example above). If 
the time step is too large, the time integration can diverge, which leads to oscillations or overflow in 
the results (see Example Problem 7a). The maximum permissible time step size that still allows 
convergence is problem dependent; users are encouraged to experiment with this value for their 
problems of interest. 

3) ITMA x=itmax and ERR =errtol are required only in the special case when one or more constituent 
material is modeled a) using classical incremental plasticity (cmod= 21) or b) implicit 
multimechanism GVIPS (cmod= 22). These two constitutive models require iterations (at each 
increment of the applied loading) to ensure that global equilibrium is satisfied. 

a) In the case of incremental plasticity, to quantify convergence of the global iterations, MAC/GMC 
4.0 employs the local effective plastic strain increment as a measure of the error. Convergence is 
achieved when the fractional change of the effective plastic strain increment from one iteration to 
the next is less than errtol throughout the RUC or laminate. 

b) In the case of implicit multimechanism GVIPS, a global energy related increment is employed as a 
measure of error. Convergence is achieved when the fractional change of this increment from one 
iteration to the next is less than errtol. The acceptable value of errtol is problem dependent. Users 
are encouraged to experiment with this value for their specific problems of interest. Note that this 
constitutive model employs a local implicit integration scheme internally. 

4) NLE G=nleg and NINTE G=nint are required only in the special case of doubly periodic RUC analysis 
using HFGMC (modid= 12). In this case, the inelastic strain field can vary arbitrarily throughout the 
RUC. In order to represent the inelastic strain field in HFGMC, MAC/GMC 4.0 employs Legendre 
polynomials, the order of which must be specified. In addition, the number of integration points at 
which all field variables are tracked must be specified, nint is the number of integration points for 
both the X 2 - and x^-directions per subcell. Thus, for example, if NINTEG=11 is specified, 11 
integration points are used per subcell in each of the two direction, for a total of 121 integration points 
for each subcell. 

5) The values nleg=5 and ninteg=11 have been shown to yield good precision in the local fields 
predicted by HFGMC. These values should be employed as default values unless the user is 
confident in using other values. 
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Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 2 (*constituents) 
MAC/GMC 4.0 Keywords Manual Section 3 (*RUC) 

MAC/GMC 4.0 Theory Manual Section 7.1 


II. Predictor-Corrector Time Integration {method— 2) 


Usage: 

* SOLVER 

METHOD =method ISTM =istm istt =istt err =errtol MINSTEP =minstp 

Name Definition 

method Time integration method: 

[= 2 ] 1 = 1 — » Forward Euler (specified time step) 

= 2 — » Predictor-Corrector (self-adaptive time-step) 
istm 2 Initial mechanical loading time step 
istt 2 Initial thermal loading time step 
errtol 3 Error tolerance for time step size changes 
minstp 4 Minimum time step size 


Example: 

* SOLVER 

METHOD=2 ISTM=0 .001 ISTT=0 . 1 ERR=0.001 MINSTEP=0 . 0 0 0 0 1 

Associated MAC/GMC 4.0 Example Problemfs): 

Example 4h 


Notes: 

1) The self-adaptive predictor-corrector integration scheme cannot currently be utilized in conjunction 
with the following features: 

a) Laminate analysis (* LAMINATE) 

b) Fatigue damage analysis (*DAMAGE) 

c) Classical incremental plasticity (cmod= 21) 

d) Implicit multimechanism GVIPS (cmod= 22). 

2) The self-adaptive predictor-corrector integration scheme employs a fourth-order Runge-Kutta starter, 
an Adams-Bashforth four-step predictor, and an Adams-Moulton four step corrector (see the 
MAC/GMC 4.0 Theory Manual Section 7.2). It starts with a small time step that is permitted to 
increase, provided the error remains within acceptable limits. ISTM =istm and ISTT =istt are used to 
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* SOLVER: Time integration 


specify the initial (starting) time step sizes for mechanical and thermal loading. These are 
independent variables because the time scale is often vastly different for mechanical and thermal 
loading. The code chooses the appropriate initial time step size from istm and istt based on whether 
the initial applied loading (specified under *MECH and * THERM) involves any change in the 
mechanical or thermal load levels. If the mechanical and thermal load levels are both initially 
changing, the code will choose the smaller of istm and istt. 

3) ERR =errtol specifies the acceptable fractional difference between predicted and corrected values of 
the field variables. If the fractional difference between these values for any field variable in any 
subcell is greater than errtol , the time step size will be cut in half. Otherwise, the time step size will 
be permitted to double, provided the history of each field variable required for prediction using the 
larger time step has been stored. An errtol value of 0.001 can be used as starting point, but an 
acceptable value can be problem dependent. Users are encouraged to experiment with this value for 
their problems of interest. 

4) MINSTEP =minstp is used to prevent the time step from becoming too small. In some cases changing 
minstep to a larger value can significantly decrease execution time (i.e., speed up the execution) 
without compromising the accuracy of the code’s results. 


Additional Information: 

MAC/GMC 4.0 Keywords Manual Section 2 (*constituents) 
MAC/GMC 4.0 Keywords Manual Section 3 (*RUC) 

MAC/GMC 4.0 Theory Manual Section 7.2 
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* SURF : Yield Surface Analysis 


Purpose: 

Generate RUC or laminate yield surfaces. 


Usage: 


*SURF 

TMA x=tmax STP =stp MMAX=mmax MODE =mode TREF =tref 
OPTION =0 j, 02 , 03 , 04,0 5, OfoO y 

ISP X=ispx(l),ispx(2),ispx(3) ISPY =ispy(l),ispy(2),ispy(3) ANGINC =anginc 
EPS =eps DR =dr !SR=isr IP =ip 


Name 

tmax 1 
stp 1 
mmax 1 
mode 1 


tref 2 

O/-O7 3 

ispx(l) 4 

ispx(2), 4 

ispx(3) 

ispy(l) 4 

ispy(2), 4 

ispy(3) 

anginc 

eps 

dr 

isr 

ip 


Definition 

Ending time for each probing excursion 
Time step size for yield surface probing 
Ending load magnitude for each probing excursion 
Mode of applied loading during yield surface probing 
= 1 —> strain or midplane strain applied 
= 2 — » stress or force resultant applied 
Reference temperature at which to generate yield surface (Special Case) 

Yield surface generation options —> up to 7 specified (see Table 4.2) (Special Case) 

Controlling component for yield surface x-axis 

Slave components to controlling x-axis component (Optional) 

Controlling component for yield surface y-axis 

Slave components to controlling y-axis component (Optional) 

Angle increment in stress space for yield surface probing 
Yield value for surface of equivalent plastic strain 
Yield value for surface of constant dissipation rate 
Yield value for surface of constant inelastic strain rate 
Yield value for surface of constant inelastic power 


Example: 

*SURF 

TMAX=400. STP=0 . 01 MMAX=0.04 MODE=l TREF=650. 

OPTION=l, 2, 3, 4 

ISPX=1 ISPY=2 , 3 ANGINC=5 . 

EPS=0 .0001 DR=0. 00025 ISR=0. 00002 IP=0.001 


Associated MAC/GMC 4.0 Example Problem (s): 

Example 4f, Example 4g 
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*surf : Yield Surface Analysis 


Notes: 

1) MAC/GMC 4.0 generates yield surfaces by “probing” in a particular stress space until all yield 
criteria are satisfied. Probing involves applying multi-axial loading such that the angle in the 
specified stress space remains constant. For example, if the Gn-o 2 2 stress space were specified, 
probing at an angle of 45° would involve incremental application of <j n = a 22 until all yield criteria 
were satisfied. Probing at an angle of 30° would involve application of g 2 2 = an tan(30°) until all 
yield criteria were satisfied. When all angles from 0° to 360° are probed and the stresses at which 
yield is detected are recorded, the yield surface can be plotted from the data files generated by the 
code. 

The loading for yield surface analysis is controlled under the *SURF keyword. This loading must be 
isothermal and always employs the forward Euler time integration scheme. In addition, the loading 
cannot be divided into separate portions with different time step sizes. Any loading information that 
is specified using *MECH, * THERM, and *SOLVER is treated as a “preload”. This preload is applied 
prior to the start of the yield surface probing and thus can be used to generate yield surfaces that 
incorporate residual stresses, tmax and mmax are the time and load magnitude at which yield surface 
probing will stop regardless of whether or not all yield surface criteria have been satisfied. Thus, if 
the code reaches this loading end point at a particular probing angle, the yield surface data for one or 
more plots will not be valid at this angle, and a warning will be written to the output file. It is 
therefore desirable to make the loading end point large so that yielding can always occur before the 
end point is reached. The code automatically stops probing at a particular angle when all yield 
criteria are satisfied, so employing large tmax and mmax values will not adversely affect execution 
time. 

Note that in order to generalize the specified loading for the multi-axial loading required for yield 
surface probing, the specified loading (mmax) is interpreted as an invariant quantity. In particular, for 

stress and force resultant control, mmax is interpreted as Jj) , while for strain control mmax is 

interpreted as \ . This allows the loading to be applied at a constant equivalent rate, regardless of 
probing angle. 

2) YREF=tref is required in the special case that no thermal preloading has been specified. If thermal 
preloading is specified (using * THERM), the code assumes that the yield surface probing occurs at the 
ending temperature of the thermal preload. In the absence of thermal preloading, the code requires 
tref to specify the temperature for yield surface probing. 

3) The yield surface options are given in Table 4.2. At least one option must be specified and up to four 
may be specified for RUC yield surface analysis. For laminate yield surfaces, up to seven options 
may be specified. Each yield surface option specified will cause a data file to be generated. This data 
file will contain the probing angle and the global stress state at which yielding was detected at that 
probing angle for each of the four yield surface types (eps, dr, isr, ip). In the case of 0PTI0N=4 or 
0PTI0N=7, all of this data will be repeated for each individual ply or each individual subcell. Note 
that elastic subcells will be skipped. 

4) The x and y stress space specifiers (ispx(j), ispy(j )) are used indicate to the code which loading 
components will be applied to generate the yield surface, ispx(l) and ispy(l jare required. These are 
the controlling components that will occupy the x- and y- axes on the plot of the yield surface. Two 
optional slave components (ispx(2), ispx(3), ispy(2), ispy(3 )) may be specified for each of the 
controlling components. The slave components will be loaded to the same value as the corresponding 
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controlling components. For example, if a n and a 2 2 were specified as the x and y controlling 
components and a 33 were specified as a slave y component (ispx=1 ispy=2 , 3), at a probing angle 
of 30°, a 33 = a 22 = an tan(30°) would be incrementally applied. 

Note that for laminate yield surfaces, ispx(i) and ispy(i) can only take the values 1, 2, or 3, 
corresponding to in-plane loading. Thus, at most one slave component can be specified. Further, the 
yield surface analysis is only valid for symmetric laminates with no bending. 


Table 4.2 Yield surface options (specified as options) and output file names for RUC and laminate 
analyses. 


Oi 

RUC analysis 

Laminate analysis 


surface generated 

file generated 

surface 

generated 

file generated 

1 

RUC level global 
yield surface 

surf global.dat 

Laminate level 
global yield surface 

surf global.dat 

2 

First subcell yield 
surface 

surf lst.dat 

First ply yield 
surface 

surf 1st ply.dat 

3 

All subcells yield 
surface 

surf all.dat 

All plies yield 
surface 

surf all plies.dat 

4 

Yield surface for each 
individual subcell 

surf local.dat 

Yield surface for 
each individual ply 

surf plies.dat 

5 

N/A 

N/A 

First subcell yield 
surface 

surf lst.dat 

6 

N/A 

N/A 

All subcells yield 
surface 

surf all.dat 

7 

N/A 

N/A 

Yield surface for 
each individual 
subcell 

surf local.dat 


Additional Information: 

MAC/GMC 4.0 Theory Manual Section 6 
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Section 5 : Damage and Failure Keywords 


This section deals with the six MAC/GMC 4.0 keywords that enable analysis of damage and failure in 
composites materials and laminates. These are: 

* allowables — » Elastic estimation of the composite allowables 

*FAILURE_SUBCELL — » Subcell level static failure analysis 

* FAI LURE_CELL — » RUC level static failure analysis 

* DAMAGE — » Cyclic fatigue damage analysis 

* DEBOND — > Internal fiber-matrix interfacial debonding 

* CURT IN — » Curtin effective fiber breakage model for longitudinal fiber failure 

Elastic allowables estimation is the simplest form of failure analysis available in MAC/GMC 4.0. 
Composite (RUC) allowable stress components are estimated based on the global level of each individual 
stress component that causes the allowable stress in a constituent material to be met. The composite 
allowable strain components are similarly estimated. These calculations are easily performed via the 
GMC concentration equations and are independent of any simulated loading applied to the composite. As 
such, the estimated allowables are like material properties for the specified RUC because altering the 
applied loading (e.g., adding simulated residual stresses) does not have an effect. 

Static failure analysis within MAC/GMC 4.0 (via the *FAILURE_SUBCELL and * FAI LURE_CELL 
keywords) is completely different from the allowables estimation. The static failure analysis monitors the 
stress and strain components in the composite or laminate (on the subcell or RUC level) and checks for 
failures using the specified failure criteria. Thus, the static failure analysis (like an actual experimental 
test) is inherently linked to the loading applied to the RUC or laminate and the predicted failure stress and 
strain are not truly material properties. That is, for example, a uni-axial tensile test on a highly strain rate 
dependent metallic alloy performed at a strain rate of 10" 2 per second would result in a completely 
different measured strength than that resulting from the same test on the same alloy performed at a strain 
rate of 10' 6 per second. Clearly, the strength measured for this material is not a material property, but is 
linked to the loading conditions to which the material is subjected. Static failure analysis within 
MAC/GMC 4.0 is the same way. The allowables estimation is a quick and simple first pass at predicting 
composite failure stresses and strains from the constituent allowables. For more realistic failure analyses, 
the static failure capabilities within MAC/GMC 4.0 should be employed. 

The final three keywords covered in this section deal with more specialized damage and failure analysis 
techniques. * DAMAGE enables use of a simple strength reduction model (typically applied to the fiber 
constituent) and a more complex stiffness reduction model (typically applied to the matrix constituent) for 
composites and laminates subjected to cyclic loading. * DEBOND allows the effect of local fiber-matrix 
debonding, which is prevalent in some composites (most notably titanium matrix composites), to be 
included within RUC and laminate analyses. Finally, the Curtin effective fiber breakage model 
(*CURTIN) is an approach based on fiber strength statistics that account for the longitudinal stiffness 
reduction due to failure of fibers in RUCs and laminates. 
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* ALLOWABLES : Estimated elastic allowables 


%. Caution % 

These estimated allowables should not be considered "design allowables", which by definition 
have been extensively validated with experimental data 


Purpose: 

Estimate the composite allowables based on the allowables of the constituent materials and linear elastic 
concentration calculations. MAC/GMC 4.0 calculates three sets of elastic stress and strain allowable 
components, listed in order from most to least conservative: 

1) First Subcell Allowables — » Global stress/strain at which the first subcell reaches its allowable 
stress/strain level 

2) Average Material Allowables —> Global stress/strain at which the average material stress/strain for 
one of the constituents reaches its allowable stress/strain level 

3) All Subcells Allowables — » Global stress/strain at which all of the subcells have reached their 
allowable stress/strain levels 

Each allowable component is estimated based purely on a linear elastic concentration calculation 
involving application of a single global stress or strain component. Three separate sets of allowables are 
provided, an upper bound (all subcells allowables), lower bound (first subcell allowables), and an 
intermediate set (average material allowables), thus enabling the user to make a judicious choice based on 
the specific problem of interest. Note that the allowables estimation is available for RUC analysis only 
(not laminate analysis). 


Usage: 

* ALLOWABLES 

NMAT =nmat 
MAT =matnum 

Sll=sll S 2 2 =s22 S33 =s33 S23 =s23 S13=sl3 S12=sl2 COMPR =comprs & 
SCll=scll SC22 =sc22 SC33 =sc33 

Ell=ell E22 =e22 E33 =e33 E23 =e23 E13=el3 E12=el2 COMPR =compre & 
ECU =ecll EC22 =ec22 EC33 =ec33 


-> Repeated nmat times 
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* allowables : Estimated elastic allowables 


Name 

nmat 
matnum 1 
sll, s22, 
s33 2 
s23, sl3, 
si 2 2 
comprs 


sell, 
sc22, 
sc33 2 ’ 3 
ell, e22, 
e33 2 
e23, el 3, 
el2 2 ’ 3 
compre 


ecll, 
ec22, 
ec33 2 ’ 4 


Definition 

Number of constituent materials for which allowables are specified 
Material number of constituent whose allowables are being specified 
Tensile normal stress allowables for material number matnum 

Shear stress allowables for material number matnum 

Compressive normal stress allowable specifier for material number matnum : 

= DIF — > Compressive stress allowables are different than tensile and must be specified 
= sam Compressive stress allowables are the same as tensile should not be specified 
= OFF — » Do not utilize compressive stress allowables, and they should not be specified 
Compressive normal stress allowables (Special Case) 


Tensile normal strain allowables for material number matnum 
Shear strain allowables for material number matnum 

Compressive normal strain allowable specifier for material number matnum : 

= DIF — > Compressive strain allowables are different than tensile and must be specified 
= SAM — » Compressive strain allowables are the same as tensile should not be specified 
= OFF — » Do not utilize compressive strain allowables, and they should not be specified 
Compressive normal strain allowables (Special Case) 


Examples: 


* ALLOWABLES 

NMAT =2 


MAT=1 
Sll=3500 . 
El 1=0 . 009 
MAT =2 
Sll=80 . 

El 1=0 . 023 


S22=91 . 2 S33 

E22=0 . 012 E33 

S22=80 . S33 

E22=0 . 023 E33 


91.2 S23=31 . 4 

0.012 E23=0 .012 

80. S23=40 . 

0.023 E23=0 .031 


S13=134 SI 2= 
E13=0 .00 9 E12 : 

S13=4 0 . S 1 2 ■ 

E13=0 .031 E12 : 


34. COMPR=SAM 
0.009 COMPR=SAM 

40. COMPR=SAM 
0.031 COMPR=SAM 


* ALLOWABLES 

NMAT=1 
MAT =3 

Sll=2500 . S12=500 . COMPR=OFF 

El 1=0 . 0 1 E22=0 . 012 E33=0.012 E23=0.012 E13=0.009 E12=0.01 COMPR=DIF & 
EC1 1=0 . 005 


Associated MAC/GMC 4.0 Example Problemfs): 

Example 5a, Example 5b 
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Notes: 

1) If a constituent material that is present in the RUC does not have allowables specified, the subcells 
that contain this material will not be considered in the allowables estimation. In such a case, the All 
Subcells Allowables will still be calculated by skipping the subcells containing the material in 
question. 

2) If a particular component is omitted from either the stress or strain allowables input for a particular 
material, that component will not be considered in the calculation for subcells containing that 
particular material. If none of the materials within the RUC have an allowable specified for a 
particular component, an allowable will not be calculated for that component. 

3) The specified shear strain allowables are interpreted by the code as engineering shear strains. 

4) The compressive stress and strain allowables {sell, sc22, sc33, ecll, ec22, ec33 ) are only required in 
the special case when the compression specifier {comprs or compre) =DIF. Otherwise, the 
compressive stress and strain allowables should be omitted. 


Additional Information: 


MAC/GMC 4.0 Theory Manual Section 5.1 
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* failure subcell : Subcell static failure 


* FAILURE SUBCELL: Subcell static failure 


Purpose: 

Perform static failure analysis based on the local stress and strain fields in the subcells. Static failure 
analysis involves checking the specified failure criteria during application of a simulated thermo- 
mechanical loading history. This is distinct from the allowables estimation, which is a linear elastic 
concentration calculation independent of the applied loading and any material history dependence. Thus, 
while the allowables may be considered predicted material properties for a composite, the static failure 
predictions are inherently linked to the specifics of the simulated loading applied to the composite. The 
maximum stress, maximum strain, and Tsai-Hill failure criteria are available within MAC/GMC 4.0. 
This availability, however, does not imply endorsement of these criteria by any parties associated with the 
code. The static failure analysis capabilities should be used with caution. 


Usage: 

* FAI LURE_SUBCELL 

NMAT =nmat 

MAT=matnum NCRIT =ncrit 

CRIT =crit Xll=xll X22 =x22 X33 =x33 X23 =x23 X13=xl3 X12=xl2 COMPR =compr & 
XCll=xcll XC22 =xc22 XC33 =xc33 ACTION =action 


Repeated ncrit times 


Repeated nmat times 


Name 

nmat 

matnum 

ncrit 

crit 


xll, x22, 
x33 2 
x23, xl3, 
xl2 2 ’ 3 
compr 


Definition 

Number of constituent materials for which failure parameters are specified 
Material number of constituent whose failure parameters are being specified 
Number of failure criteria to check for material number matnum 
Failure criterion material number matnum : 

= 1 — » Maximum stress failure criterion 
= 2 — » Maximum strain failure criterion 
= 3 — » Tsai-Hill multi-axial failure criterion 
Tensile normal failure parameters for material number matnum and criterion number crit 
(failure stresses for crit= 1, or 3, failure strains for crit= 2) 

Shear failure parameters for material number matnum and criterion number crit 
(failure stresses for crit= 1, or 3, failure strains for crit= 2) 

Compressive normal failure specifier for material number matnum and criterion number crit : 
= DIF — > Compressive parameters are different than tensile and must be specified 
= sam — » Compressive parameters are the same as tensile should not be specified 
= OFF — » Do not utilize compressive parameters, and they should not be specified 
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Name Definition 

xcll , Compressive normal failure parameters for material number matnum and criterion 

xc22, number crit (failure stresses for crit= 1, or 3, failure strains for crit= 2) 

xc33 2,4 (Special Case) 

action 5 Action for code to take upon detection of subcell failure: 

= 1 — » Zero subcell stiffness and continue execution 

= 0 — » Do nothing and continue (write notification to output file) 

= - 1 —> Stop execution 


Examples: 


* FAILURE_SUBCELL 

NMAT=2 


MAT=1 NCRIT=2 

CRIT=1 Xll=3500 . X22=91 . 2 X33=91.2 X23=31.4 X13=134. X12=134. & 

COMPR=OFF ACTION=l 

CRIT=2 X11=0 . 009 X22=0 .012 X33=0.012 X23=0.012 X13=0.009 X12=0.009 & 
COMPR=OFF ACTION=l 
MAT=2 NCRIT=2 

CRIT=1 Xll=80. X22=80 . X33=80. X23=40. X13=40. X12=40. & 

COMPR=OFF ACTION=l 

CRIT=2 X11=0 . 023 X22=0.023 X33=0.023 X23=0.031 X13=0.031 X12=0.031 & 
COMPR=OFF ACTION=l 


* FAILURE_SUBCELL 

NMAT=2 

MAT=4 NCRIT=1 
CRIT=3 Xll=3500 . 

X22=91 .2 

X33=91 . 2 

X23=31 . 4 X13=134 . 

X12=134 

ACTION=-l 
MAT=7 NCRIT=2 
CRIT=1 X22=80 . 

X33=80 . 

X23=4 0 . 

COMPR=DIF & 


XC22=100 . 

XC33=100 . 

XC23=60 . 

ACTION=0 


CRIT=2 X11=0 . 023 

COMPR=SAM 

ACTION=l 




& 


Associated MAC/GMC 4.0 Example ProblemfsU 

Example 5b, Example 5c, Example 5d 


Notes: 

1) If a constituent material that is present in the RUC does not have subcell failure parameters specified, 
the subcells that contain this material will not be considered when checking for failure (i.e., they are 
treated as infinitely strong). 

2) If a particular component (xij) is omitted from the failure parameters for either the maximum stress or 
maximum strain criterion, that component will not be considered when checking for failure of 
subcells containing that particular material. For the Tsai-Hill criterion, all of the components (xij) are 
required. 
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* failure subcell: Subcell static failure 


3) In the case of the maximum strain criterion, the specified shear failure strains are interpreted by the 
code as engineering shear strains. 

4) The compression specifier ( compr ) and the compressive failure parameters ( xcll , xc22, xc33 ) are not 
required for the Tsai-Hill criterion. The compressive failure parameters (xcll, xc22, xc33 ) are only 
required for the maximum stress and maximum strain criteria in the special case when the 
compression specifier compr= dif. Otherwise, the compressive failure parameters should be omitted. 

5) ACTION ^action is ignored in the case of fatigue damage analysis (* DAMAGE) and can be omitted in 
that case. 


Additional Information: 


MAC/GMC 4.0 Theory Manual Section 5.2 
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* FAILURE CELL : Repeating unit cell static failure 


Purpose: 

Perform static failure analysis based on the RUC level stress and strain fields. As with subcell level static 
failure analysis, RUC level static failure analysis involves checking the specified failure criteria during 
application of a simulated thermo-mechanical loading history. Thus, the RUC level static failure 
predictions are inherently linked to the specifics of the simulated loading applied to the composite. 


Usage: 

* FAI LURE_CELL 

NCRIT =ncrit 

CRIT =crit Xll=xll X22 =x22 X33 =x33 X23 =x23 X13=xl3 X12=xl2 COMPR =compr & 
XCll=xcll XC22 =xc22 XC33 =xc33 ACTION =action 


Repeated ncrit times 


Name 

ncrit 

crit 


xll, x22, 
x33 1 
x23, xl 3, 
xl2 lj2 
compr 


xcll, 
xc22 , 
xc33 u 
action 4 


Definition 

Number of RUC level failure criteria to check 
Failure criterion: 

= 1 —> Maximum stress failure criterion 
= 2 — » Maximum strain failure criterion 
= 3 — » Tsai -Hill multi-axial failure criterion 
Tensile normal failure parameters for criterion number crit 
(failure stresses for crit= 1 or 3, failure strains for crit= 2) 

Shear failure parameters for criterion number crit 
(failure stresses for crit= 1 or 3, failure strains for crit= 2) 

Compressive normal failure specifier for criterion number crit : 

= DIF — > Compressive parameters are different than tensile and must be specified 
= sam — » Compressive parameters are the same as tensile should not be specified 
= OFF — » Do not utilize compressive parameters, and they should not be specified 
Compressive normal failure parameters for criterion number crit (Special Case) 
(failure stresses for crit= 1, or 3, failure strains for crit= 2) 

Action for code to take upon detection of RUC failure: 

= 0 — » Do nothing and continue (write notification to output file) 

= - 1 —> Stop execution 


NAS A/TM— 2002-2 1 2077 /VOL2 


86 


Section 5: Damage and Failure Keywords 


* fai lure_cell : Repeating unit cell static failure 


Examples: 


*FAILURE_CELL 

NCRIT=1 

CRIT=3 XI 1=21 12 . X22=68 . 08 X33=68.08 X23=31.4 X13=24.89 X12=24.89 & 
ACTION=-l 


*FAILURE_CELL 

NCRIT=2 

CRIT=1 Xll=450 . X33=110 . X23=80. COMPR=DIF & 

XC22=300 . XC33=140 . XC23=60. ACTION=0 

CRIT=2 X12=0 . 03 COMPR=SAM ACTION=-l 


Associated MAC/GMC 4.0 Example Problem (s): 

Example 5c, Example 5d 


Notes: 

1) If a particular component {xij) is omitted from the failure parameters for either the maximum stress or 
maximum strain criterion, that component will not be considered when checking for failure of the 
RUC. For the Tsai-Hill criterion, all of the components {xij) are required. 

2) The specified shear strain allowables are interpreted by the code as engineering shear strains. 

3) The compression specifier ( compr ) and the compressive failure parameters {xcll, xc22, xc33) are not 
required for the Tsai-Hill criterion. The compressive failure parameters {xcll, xc22, xc33 ) are only 
required for the maximum stress and maximum strain criteria in the special case when the 
compression specifier compr= dif. Otherwise, the compressive failure parameters should be omitted. 

4) ACTION ^action is ignored in the case of fatigue damage analysis (when * DAMAGE is included in the 
input file) and can be omitted in that case. 


Additional Information: 


MAC/GMC 4.0 Theory Manual Section 5.2 
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*DAMAGE : Fatigue damage analysis 


Purpose: 

Perform fatigue damage analysis on an RUC or laminate. The damage related applied loading 
information and constituent material damage parameters are also specified under * DAMAGE. MAC/GMC 
4.0 contains two fatigue damage models. The first is a simple strength reduction model, which allows for 
a specified semi-logarithmic degradation of material strength, typically applied to fiber constituents. The 
second is a more involved stiffness reduction model (the Anisotropic Damage Evolution and Life Model, 
which reduces to the ONERA NonLinear Cumulative Damage Rule when the material is isotropic), which 
has a number of material parameters that are characterized based on material S-N curves (see Arnold and 
Kruch, 1994). 

In traditional fatigue tests on materials, a given loading cycle is applied repeatedly to a material until 
failure occurs. The fatigue life of the material is then the number of cycles to failure for a given cyclic 
loading profile. For composite materials, the GMC micromechanics models implemented within 
MAC/GMC 4.0 are ideal for fatigue damage analysis because of the availability of the micro scale 
(subcell level) fields enables stress redistribution during fatigue simulations. That is, for example, if a 
cyclic fatigue load of 0 to 50 ksi is imposed on a virgin composite material specimen, the fiber and matrix 
will each experience a certain multiaxial stress state at the maximum load level of 50 ksi. However, after 
a number of cycles, matrix damage may cause a local stiffness reduction. Consequently, during 
subsequent cycles, the fiber stress will increase due to the loss of load carrying capacity in the matrix. 
This type of stress redistribution cannot be accounted for using macro scale fatigue approaches (i.e., 
treating the composite as a homogenized effective material) in conjunction with Basquin’s Law or the 
Coffin-Manson Law. 

Explicitly simulating fatigue tests to account for the stress redistribution by repeatedly executing a model 
to evaluate the damage for every cycle can quickly become computationally intensive. This is 
particularly true in the present case of inelastic micromechanics where the loading for each cycle must be 
applied incrementally, and a good number of local field and state variables must be tracked for each 
subcell. The fatigue damage analysis within MAC/GMC 4.0 avoids this computationally demanding 
explicit type of analysis by solving the inverse problem. Rather than applying a loading cycle and 
determining the damage, a damage increment is applied, and the number of cycles required to cause this 
damage increment is calculated. This procedure preserves the ability to capture stress redistribution, but 
cuts the number of loading cycles that must be explicitly simulated from thousands or hundreds of 
thousands to tens. 

The MAC/GMC 4.0 fatigue damage analysis procedure is as follows for a specified loading cycle and 
damage increment (AD), where locally, the damage parameter, D, varies between an undamaged level of 
0 and a completely damaged, or zero stiffness level of 1.0: 

1) Apply specified thermo-mechanical loading cycle to the composite or laminate 

2) Determine number of cycles (AN) required to increment the local damage ( D ) in the first (most 
damaged) controlling subcell by AD. 

3) Determine damage level in all other subcells due to AN number of cycles. 

4) Add AN to the previous total number of cycles (N) to arrive at new total number of cycles. 

5) Re-apply the specified thermo-mechanical loading cycle, go to step 2. 
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* damage : Fatigue damage analysis 


Clearly, when the specified loading is applied to the damaged material in step 5, the resulting stress field 
throughout the damaged composite or laminate will be different than that experienced by the undamaged 
composite or laminate after step 1 . This redistribution of stresses then affects the evolution of damage as 
the damage increment (AD) is repeatedly applied and the total number of cycle (TV) increases. When a 
particular subcell reaches a damage level of 1.0, it has failed and (provided the RUC or laminate has not 
completely failed) a new controlling subcell will be chosen whose damage level will be incremented. 


Usage: 

Strength Reduction Model 

* DAMAGE 

MAXNB =maxnb DINC =dinc DMA X=dmax BLOCK=tl,t2 
NDMAT =ndmat 

MAT=matnum MOD =dmod SBl=sul(l 1), sul(22), sul(33), sul(23), sul(13), sul(12) & 
SU2=su2(l 1), su2(22), su2(33), su2(23), su2(13), su2(12) & 

Nl=nl(ll),nl(22), nl(33), nl(23), nl(13), nl(12) & 

N2=n2(l l),n2(22), n2(33), n2(23), n2(13), n2(12) 


Repeated ndmat times 


Stiffness Reduction Model 

* DAMAGE 

MAXNB =maxnb DINC =dinc DMA X=dmax BLOCK=tl,t2 
NDMAT =ndmat 

MNY=matnum MOD =dmod AN G=0 BN =b BP =b’ OMU =(D U OMFL = COfl OMM=&> m & 
etu=77 m etfl =/^7 E r m=r/ m be =j3 a =a sfl=o /7 xml =M su =cr u 


-» Repeated ndmat times 

Name Definition 

maxnb 1 Maximum number of load blocks to apply 

dine 2 Damage increment to apply for each load block 

dmax 3 Maximum amount of damage permitted (0 < dmax < 1 .0) 

tl,t2 4 Start and end times (from thermo-mechanical loading profile) to be associated with a load 
block (see Figure 5.1) 

ndmat Number of materials for which damage parameters are specified 
matnum Material number (from CONSTITUENTS) whose damage parameters are specified 
dmod Damage model to employ for material number matnum 

[=1] =1 — » Strength reduction fatigue damage model 

= 2 — » Stiffness reduction fatigue damage model 
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Name 

sul(ij) 

su2(ij) 

nl(ij) 


n2(ij) 


Qi.l 

b 6 ’ 1 
b ’ 6,7 


COfl 

co m 


5,7 

5,7 

5,7 


% 


5,7 


% 

1jm 


5,7 

5,7 


P 


6,7 
. 6,7 


M 

C>u 


6,7 

6,7 

6,7 


Definition 

Ultimate strength component ij point 1 from semi-log strength reduction plot 
(see Figure 5.2) 

Ultimate strength component ij point 2 from semi-log strength reduction plot 
(see Figure 5.2) 

Number of cycles point 1 for component ij from semi-log strength reduction plot (see 
Figure 5.2) 

Number of cycles point 2 for component ij from semi-log strength reduction plot (see 
Figure 5.2) 

Angle between x r axis and direction of transverse isotropy 
Quantifies effect of mean stress on M: M=M 0 { 1 - b<j mean ) 

Quantifies effect of mean stress on oft: <Jfi=(<Jfl)o (1 - bcr mean ) 

Ratio of longitudinal to transverse normal strength 

Ratio of longitudinal to transverse normal fatigue limit stress 

Ratio of longitudinal to transverse normalizing normal stress amplitude 

Ratio of longitudinal to transverse shear strength 

Ratio of longitudinal to transverse shear fatigue limit stress 

Ratio of longitudinal to transverse normalizing shear stress amplitude 

Exponent on the normalized stress amplitude 

Scaling factor for the stress dependency in the fatigue damage exponent (d=acr u /ofi) 
Normal stress fatigue limit 

Ratio of the stress amplitude constant to the ultimate strength 
Ultimate normal strength 


Examples: 

Strength Reduction Model 

♦DAMAGE 

MAXNB=100 DINC=0 . 2 DMAX=1 . 0 BLOCK=0 . , 2 0 0 . 

NDMAT=1 

MAT=1 MOD=2 SU1=3500, 91 .2, 91 .2, 31 . 4, 134 . , 134 & 

SU2=2000 . , 91. 2, 91. 2, 31. 4, 134. ,134. & 

Nl=1000, 1000, 1000, 1000, 1000, 1000 & 

N2=300000000, 300000000, 300000000, 300000000, 300000000, 300000000 


Stiffness Reduction Model 

♦DAMAGE 

MAXNB=100 DINC=0 . 2 DMAX=1 . 0 BLOCK=0 . , 2 0 0 . 

NDMAT=1 

MAT=2 MOD=l ANG=0 . BN=0 . 0 BP=0 . 0 OMU=l . OMFL=l . OMM=l . ETU=1 . & 

ETFL=1 . ETM=1. BE=9. A=0 . 05 SFL=27 . XML=150 . & 

SU=80 . 


Associated MAC/GMC 4.0 Example ProblemfsV 

Example 5d 
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♦damage : Fatigue damage analysis 


Notes: 

1) Execution of a MAC/GMC 4.0 fatigue damage analysis will end when all subcells (or RUCs) have 
failed, the code determines that the remaining fatigue life is infinite, or maxnb load blocks have been 
applied. 

2) A smaller damage increment (dine) value will generally cause the code to apply a greater number of 
load load blocks prior to failure, leading to a larger execution time. 

3) The maximum value of dmax actually employed by the code is 0.99999. This value is employed 
when a dmax value of 1.0 is specified, which is a good default value as it allows the subcells to 
damage almost completely. A D value of exactly 1.0 causes a zero stiffness for the subcell and leads 
to numerical difficulties. 

4) As shown in Figure 5.1, the load block that constitutes a single cycle of the applied fatigue loading is 
associated with a portion of the thermo-mechanical loading specified under *MECH and * THERM by 
defining tl and t2. Because of the cyclic nature of fatigue loading, the thermal loading profile 
specified under * therm must have the same temperature at tl and t2 so that the cycle starts and ends 
at the same temperature. Likewise, each component of the mechanical loading specified under 
*MECH must have the same magnitude at tl and t2. 

5) The fatigue damage model parameters, 6 , co w cty, co m , rj u , r/fl, and r/ m are applicable only to cases 
involving transversely isotropic subcell materials. For isotropic subcell materials, the specified value 
of 6 is arbitrary, while co w oty, co m , rj u , r/fl, and rj m should all be set equal to 1.0. 

6) For isotropic subcell materials, the parameters, b, b /3, &, oft, M, and cr u must be characterized for the 
materials based on experimental material S-N curves. If no mean stress dependence is desired, b and 
b’ can be set to zero. See MAC/GMC 4.0 Example Problem 5d for an example of how the stiffness 
reduction fatigue damage model parameters are characterized. 

7) For additional information on the stiffness reduction fatigue damage model material parameters, see 
the MAC/GMC 4.0 Theory Manual Section 5.3 and Arnold and Kruch (1994). 

8) If any local failures or interfacial debonding are detected during application of a particular load block, 
the load block is repeated and one cycle is added to the total number of cycles. 

9) Currently both the strength and the siffness reduction fatigue damage models are isothermal and do 
not admit temperature-dependent material properties. 


Additional Information: 


MAC/GMC 4.0 Keywords Manual Section 4 (*mech) 
MAC/GMC 4.0 Keywords Manual Section 4 (*therm) 
MAC/GMC 4.0 Theory Manual Section 5.3 
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Figure 5.1 The MAC/GMC 4.0 fatigue damage analysis requires the specification of a thermo- 
mechanical “load block”, which defines a cycle of applied loading in the fatigue damage 
simulation. The load block is chosen from the thermo-mechanical loading specified under 
*MECH and * THERM by specifying BLOCKS/, £2. (a) Simple example; (b) More involved 
example. 
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* damage : Fatigue damage analysis 



Figure 5.2 MAC/GMC 4.0 strength reduction fatigue damage model. Each component of the 
material’s ultimate strength is reduced according to the linear behavior defined by the 
points sul, su2, nl, and n2 from a semi-logarithmic plot. 
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* DEBOND : Fiber-matrix interfacial debonding 


Purpose: 

Enable interfacial debonding between the subcells within an RUC or laminate analysis. Two debonding 
models are available within MAC/GMC 4.0. Both models treat normal and shear interfacial debonding 
independently and incorporate finite shear and normal bond strengths. The constant compliant interface 
(CCI) model is formulated such that, after debonding occurs at a particular interface, the interfacial stress 
becomes a constant (for properly chosen model parameters). The evolving compliant interface (ECI) 
model, on the other hand, allows the interfacial stress to unload and be redistributed after debonding 
occurs (again, for properly chosen model parameters). While both models are made available, the ECI is 
thought to capture the physics of debonding better than the CCI model, especially when used in concert 
with HFGMC, which allows more accurate local stress redistribution than GMC. However, the CCI 
model may be useful in some instances in which one desires to maintain a certain amount of load carrying 
capacity subsequent to failure of the interface. 

Usage: 

* DEBOND 

L =ly Nil =nii 

DBCH =dbch nai =nai NBI =nbi NGI =ngi FACE =face bdn=<j db ln= n BN= B n TOLN =toln & 

BDS= x DB LS = A S BS= B s DELAY =delay 


-> Repeated nil times 


-» Repeated nly times 

Name Definition 

ly 1 Layer number (Special Case) 

nii Number of debonding interfaces for layer ly 

dbch Debonding model choice 

= 1 — » Constant compliant interface (CCI) model 
= 2 — » Evolving compliant interface (ECI) model 

nai 2 a index of subcell to which debonding model will be applied (see Figure 3.2) (Special Case) 

nbi (3 index of subcell to which debonding model will be applied (see Figure 3.1 and Figure 3.2) 

ngi y index of subcell to which debonding model will be applied (see Figure 3.1 and Figure 3.2) 

face Interface to which debonding model will be applied (see Figure 5.3 and Figure 5.4) 

<j DB 3 Normal interfacial bond strength 

A n 3,4 ECI/CCI normal debonding parameter A 

B n 3,4 ECI/CCI normal debonding parameter B 

toln 3,5 Reversal tolerance - stress below which interface is considered to be in compression 
x DB 3 Shear interfacial bond strength 

A s 3,4 ECI/CCI shear debonding parameter A 

B s 3,4 ECI/CCI shear debonding parameter B 

delay 6 Time delay - time before which interface is treated as perfectly bonded (sec.) 
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* debond : Fiber-matrix interfacial debonding 



face = 


1 

2 

3 


debonded facets) 

Top face 
Right face 
Both faces 


Figure 5.3 Definition of debonded faces for a given subcell within a doubly periodic RUC. 




face = 

debonded face(s) 

face = 

debonded face(s) 

i 

Top face 

5 

Top and Back faces 

2 

Right face 

6 

Right and Back faces 

3 

Back faces 

7 

All faces 

4 

Figure 5.4 

Top and Right faces 

Definition of debonded faces for a given subcell within a triply periodic RUC. 
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Examples: 

*DEBOND 

NII=2 

DBCH=1 NBI=1 NGI=1 FACE=2 BDN=7 . LN=0.0001 BN=60 . TOLN=0 . & 
BDS=40 LS=0 . 1 BS=100 . DELAY=64800. 

DBCH=2 NBI=3 NGI=1 FACE=2 BDN=7 . LN=0.0001 BN=8 . TOLN=0 . & 
BDS=40 LS=0 . 1 BS=100 DELAY=64800. 

*DEBOND 

LY=1 NI 1=1 

DBCH=1 NBI=1 NGI=1 FACE=2 BDN=10. LN=0.0001 BN=10. TOLN=0 . & 
BDS=40 LS=0 . 1 BS=100 . DELAY=0 . 

LY=2 NII=0 
LY=3 NI 1=1 

DBCH=1 NBI=1 NGI=1 FACE=2 BDN=10. LN=0.0001 BN=10. TOLN=0 . & 
BDS=40 LS=0 . 1 BS=100. DELAY=0 . 


Associated MAC/GMC 4.0 Example Problem (s): 

Example 5e 


Notes: 

1) L =ly is only required in the special case of laminate analysis. In this case, the debonding parameters 
for the interfaces in the RUCs representing each layer must be specified. If a layer has no debonding 
interfaces, LY =ly Nll=0, should be specified (see Example above). 

2) NAI =nai is only required in the special case of a triply periodic RUC. The triply periodic RUC may 
be part of a stand alone triply periodic RUC analysis or may represent a particular layer in a laminate 
analysis. 

3) The units of cr DB , toln , and t db must be consistent with the stress units employed throughout the 
analyzed problem. The units of B n and B s are seconds, while the units of A n and A s are l/(stress 
units). 

4) The debonding parameters, B m A S9 and B s are constants in evolution equations that affect the 
debonding behavior of the interface. The form of these evolution equations is different for the CCI 
and CCI models. See the MAC/GMC 4.0 Theory Manual Section 5.4 for details on the debonding 
model equations. 

5) The reversal tolerance, toln , can be used to force the interface to “close” when the normal stress is 
smaller than the specified value. The debonding model implementations within MAC/GMC 4.0 are 
such that, when the normal stress at a debonded interface becomes compressive, the interface “closes” 
and reacquires its original (perfectly bonded) normal stiffness. The reversal tolerance sets a small 
tensile value for the interfacial stress for the code to consider negligible. 
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6) The time delay, delay, can be used to artificially prevent debonding of the interface until after the 
specified time. This can be useful in situations where the user wished to disallow debonding during a 
portion of the applied loading history. 

Additional Information: 

MAC/GMC 4.0 Theory Manual Section 5.4 
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* CURT IN : Curtin effective fiber breakage model 


Purpose: 

Enable the Curtin effective fiber breakage model for subcells within a doubly periodic RUC (the model is 
not applicable to triply periodic RUCs). The Curtin model simulates the effects of fiber fracture and the 
unloaded portion of the fiber in the vicinity of the fiber fracture. The model uses statistical fiber strength 
data in order to simulate the effects of the breaking and subsequent load shedding of many fibers as one 
effective fiber. The effective fiber stiffness is degraded as the longitudinal stress in the fiber is increased. 
The effective longitudinal fiber stress reaches a maximum, which can be used as a failure stress for the 
effective fiber. 

Usage: 

* CURT IN 

L =ly NCURT =ncurt 

NBI =nbi NGI =ngi D =d LO =l 0 SIG0= cr 0 TAU0= x 0 M=m DELAY =delay ACTION =action 


-> Repeated ncurt times 


-» Repeated nly times 

Name Definition 

ly 1 Layer number (Special Case) 
ncurt Number of Curtin model effective fibers for layer ly 
nbi P index of subcell to which Curtin model will be applied (see Figure 3.1) 

ngi y index of subcell to which Curtin model will be applied (see Figure 3.1) 

d 2 Fiber diameter 

l 0 2 Fiber characteristic length 

g 0 3 Characteristic (average) fiber strength 
x o 3,4 Shear stress sliding resistance at fiber-matrix interface 

m Weibull modulus from fiber strength statistics 

delay 5 Time delay - time before which Curtin model is inactive 

action Action for code to take upon detection of Curtin model effective fiber failure: 

= 1 — » Zero fiber stiffness and continue execution 

= 0 — » Do nothing and continue (write notification to output file) 

= - 1 —> Stop execution 

Examples: 

*CURTIN 

NCURT=1 

NBI=1 NGI=1 D=142 . E-6 L0=0.0127 SIG0=508. TAU0=2.03 M=17.0 & 
DELAY=64800 . ACTION=0 
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* curt in : Curtin effective fiber breakage model 


* CURT IN 

LY=1 NCURT=0 
LY=2 NCURT=1 

NBI=1 NGI=1 D=142 . E-6 L0=0.0127 SIG0=610. TAU0=6.78 M=10.0 & 
DELAY=0 . ACTION=0 
LY=3 NCURT=0 

Associated MAC/GMC 4.0 Example Problemfs): 

Example 5f 


Notes: 

1) L =ly is only required in the special case of laminate analysis. In this case, the Curtin model 
parameters for the fibers in the doubly periodic RUCs representing each layer must be specified. If a 
layer has no Curtin model fibers, LY =ly NCURT=0, should be specified. 

2) The units for the fiber diameter and characteristic length, d and l 0 , must have consistent (length) units, 
but these units do not have to be consistent with the remaining material property units. 

3) The units of g 0 and must be consistent with the stress units employed throughout the rest of the 
problem. 

4) The shear stress sliding resistance at the fiber-matrix interface, x 0 , is the most difficult of the Curtin 
model parameters to obtain. Researchers have often attempted to deduce a value for this parameter 
from fiber push out and pull out tests on the composite. This parameter may also be chosen based on 
characterization vs. a known longitudinal composite stress-strain curve that is taken to failure. 

5) The time delay, delay , can be used to artificially prevent activation of the Curtin model until after the 
specified time. This can be useful in situations where the user wishes to disallow fiber stiffness 
degradation and/or breakage during a portion of the applied loading history. 


Additional Information: 


MAC/GMC 4.0 Theory Manual Section 5.5 
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Section 6 : Results and Data Output Keywords 


This section presents the four MAC/GMC 4.0 keywords that are related to the results of the analyzed 
problem. These are: 

* print — » Specify the level of results data to be written to the MAC/GMC 4.0 output file 

*XYPLOT — » Specify data to be written to files for generating x-y scatter plots 

*PATRAN — » Generate results output for post-processing with MSC/PATRAN software 

★MATLAB — > Generate results output for post-processing with MATLAB software 

The lowest level of data output in MAC/GMC 4.0 is controlled with the *PRINT keyword, which 
specifies the amount of data written to the ASCII format MAC/GMC 4.0 output file. Next, the *XYPLOT 
keyword writes time -based output to ASCII files that can be used within standard plotting software 
packages to generated x-y plots (such as predicted composite stress-strain curves). Finally, *PATRAN and 
* matlab enable the generation of data files to be read by the MSC/PATRAN and MATLAB (the 
MathWorks, Inc.) software packages for generating fringe plots that show the variation of the result fields 
over the composite architecture. Generation of fringe plots with these software packages requires the 
appropriate MAC/GMC 4.0 files (distributed with the MAC/GMC 4.0) that work in concert with 
*PATRAN and * MATLAB. 
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* print : Overall print output level for output file 


* PRINT : Overall print output level for output file 


Purpose: 

Specify the level of output data to be written to the MAC/GMC 4.0 ASCII output file 

Usage: 

* PRINT 

NPL=npl 

Name Definition 

npl Print level specifier for the MAC/GMC 4.0 ASCII output file (see Table 6.1) 


Table 6.1 Data written to the output file based on the value of npL 


npl 

Result 

-i 

Determine and write effective properties and stop code execution 

0 

Minimal output, no effective properties, no time -based output 

1 

Time based global fields output, no effective properties 

3 

Effective RUC or laminate properties and architecture information 

4 

Debonding, yield surface, and static failure analysis information 

6 

Effective stiffness matrix and CTEs or ABD matrix output at each time step 

7 

Extensive fatigue damage information 

8 

Iteration data to screen for incremental plasticity and multimechanism GVIPS 

9 

Ancillary data 

10 

Program execution trace and all array data 

% Caution - Generates large output file and will slow execution considerably % 


Example: 

♦PRINT 

NPL=6 

Associated MAC/GMC 4.0 Example Problem (s): 

All example problems 


Notes: 

1) A print level of npl is inclusive of all output information associated with all print levels lower than 
npl. 

2) If * print is not specified, the code will use a default value of npl= 1. 
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*XYPLOT : x-y plot output 


Purpose: 

Generate time -based data for x-y scatter plots. 


Usage: 

*XYPLOT 

FRE Q=freq 
LAMINATE =nlam 

NAME =name_l X=x_l Y =y_l 


-> Repeated nlam times 
MACRO=nmacro 

NAME=name_mac LYR=lyr_mac X=x_mac Y =y_mac COORD =coord 


-> Repeated nmacro times 
MICRO =nmicro 

NAME =name_mic LYR=lyr_mic I A =ia IB =ib 


IG =ig X=x_mic Y=y_mic 


-> Repeated nmicro times 


Name 

freq 
nlam 1 

rl,2 

name_i 

x_l l 

yj' 


nmacro 
name mac 2 
lyrjnac 1 
xjnac 3 

yjnac 3 


Definition 

Frequency of time steps to write output (i.e., write output every freq time step) 
Number of laminate level plot files (Special Case) 

Laminate level plot file name without path (Special Case) 

I.D. number of data to write in first column of data file for current laminate plot 
(see Table 6.2) (Special Case) 

I.D. number of data to write in second column of data file for current laminate plot 
(see Table 6.2) (Special Case) 

Number of macro (RUC level) plot files 
Macro level plot file name without path 
Layer number for current macro plot (Special Case) 

I.D. number of data to write in first column of data file for current macro plot 
(see Table 6.3) 

I.D. number of data to write in second column of data file for current macro plot 
(see Table 6.3) 
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*xyplot : X-Y plot output 


Name 

coord 1,4 


nmicro 

namejnac 2 

lyrjnic 

. 5 

ia 

ib 

ig 

x mic 3 


yjnic 


Definition 

Specifies local (principal material) or global (laminate) coordinate system for macro (RUC 
level) plots within laminate layer (Special Case) (Optional): 

= L — > Local coordinate system for layer (default if COORD =coord is omitted) 

= G — » Global (laminate) coordinate system 
Number of micro (subcell level) plot files 
Micro level plot file name without path 
Layer number for current micro plot (Special Case) 
a subcell index for current micro plot (see Figure 3.2) (Special Case) 

(3 subcell index for current micro plot (see Figure 3.1 and Figure 3.2) 
y subcell index for current micro plot (see Figure 3.1 and Figure 3.2) 

I.D. number of data to write in first column of data file for current micro plot 
(see Table 6.4) 

I.D. number of data to write in second column of data file for current micro plot 
(see Table 6.4) 


Examples: 

*XYPLOT 

FREQ=5 

MACR0=1 

NAME=example_3a X=2 Y=8 
MICRO=0 

*XYPLOT 

FREQ=1 

LAMINATE=1 

NAME = kx y_v_Mx y X=9 Y=15 
MACR0=2 

NAME =my_p lot_file LYR=1 X=2 Y=8 COORD=G 
NAME=" LYR=3 X=2 Y=14 
MICR0=1 

NAME=micro file LYR=2 IB=1 IG=1 X=2 Y=8 


Associated MAC/GMC 4.0 Example Problem (s): 

Example lc, Example 2a-h, Example 3a-i, Example 4a-d, Example 4h, Example 5b, Example 5d-f, 
Example 6a, Example 6b-c, Example 7a-c 


Notes: 

1) The x-y plot information associated with laminate level plots should only be specified in the special 
case of laminate analysis. Further, the specification of LYR =lyr_mac, LYR=lyr_mic, and 
COORD =coord for the macro and micro plots is only required in the special case of a laminate in order 
to indicate which layer the RUC or subcell information should be taken from. 
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2) Laminate level plot files have “_lam . data” appended to the specified namej. Macro level plot files 
have “_macro . data” appended to the specified namejnac. Micro level plot files have 
“_micro . data” appended to the specified namejnic. The files will be written in the same directory 
as the input file. If a name is repeated within a particular level of plots, an integer counter will be 
appended to the specified file name. Finally, the special character double quotation (“) can be 
specified as a file name to automatically repeat the file name directly preceeding this specification 
(see Example above). 

3) The macro and micro plot files associated with a laminate contain four columns of data rather than 
two. The first two columns contain the x and y data for quadrature point number 1 in the specified 
layer, while the third and fourth columns contain the x and y data for quadrature point number 2 in the 
specified layer. 

4) Even in the case of laminate analysis, COOKD=coord is optional. If it is omitted, the macro plot 
information will be written in the local coordinates of the layer. 

5) IA =ia is only required in the special case when the subcell in question is part of a triply periodic 
RUC. 
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*XYPLOT : X-Y plot Output 


Table 6.2 LAMINATE level x-y plot data file i.d. numbers. 


1 

£ ° 

° XX 

midplane normal strain 

32 

Ac 

magnetic field 

2 

£ ° 

b yy 

midplane normal strain 

33 

H y 

magnetic field 

3 

A 

average normal strain 

34 

A 

magnetic field 

4 

A 

average engineering shear strain 

35 

n e 

T XX 

electric force resultant 

5 

Txz 

average engineering shear strain 

36 

N e 

yy 

electric force resultant 

6 

r° 

/ xy 

midplane engineering shear strain 

37 

n e 

xy 

electric force resultant 

7 

K xx 

midplane curvature 

38 

m e 

XX 

electric moment resultant 

8 

K yy 

midplane curvature 

39 

m e 

yy 

electric moment resultant 

9 

K xy 

midplane curvature 

40 

m e 

xy 

electric moment resultant 

10 

Nxx 

force resultant 

41 

n m 

T XX 

magnetic force resultant 

11 

N yy 

force resultant 

42 

N m 

yy 

magnetic force resultant 

12 

A, 

force resultant 

43 

n m 

xy 

magnetic force resultant 

13 

M xx 

moment resultant 

44 

M m 

XX 

magnetic moment resultant 

14 

M yy 

moment resultant 

45 

M m 

yy 

magnetic moment resultant 

15 

M xy 

moment resultant 

46 

M m 

xy 

magnetic moment resultant 

16 

K 

inelastic force resultant 

47 

n et 

T XX 

pyroelectric force resultant 

17 

N in 

yy 

inelastic force resultant 

48 

N et 

yy 

pyroelectric force resultant 

18 

N in 

xy 

inelastic force resultant 

49 

n et 

xy 

pyroelectric force resultant 

19 

M in 

XX 

inelastic moment resultant 

50 

m et 

XX 

pyroelectric moment resultant 

20 

M in 

yy 

inelastic moment resultant 

51 

m et 

yy 

pyroelectric moment resultant 

21 

M in 

xy 

inelastic moment resultant 

52 

m et 

xy 

pyroelectric moment resultant 

22 

n t 

thermal force resultant 

53 

n mt 

T XX 

pyromagnetic force resultant 

23 

n t 

yy 

thermal force resultant 

54 

n mt 

yy 

pyromagnetic force resultant 

24 

K 

thermal force resultant 

55 

n m 

xy 

pyromagnetic force resultant 

25 

Kx 

thermal moment resultant 

56 

M mt 

XX 

pyromagnetic moment resultant 

26 

K 

thermal moment resultant 

57 

yy 

pyromagnetic moment resultant 

27 

K 

thermal moment resultant 

58 

M mt 

xy 

pyromagnetic moment resultant 

28 

N/A 

N/A 

59 - 99 

N/A 

29 

A 

electric field 

100 

T 

temperature 

30 

A 

electric field 

101 

t 

time 

31 

A 

electric field 

102 

^ creep 

creep time 
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Table 6.3 MACRO level x-y plot data i.d. numbers. 


1 


total normal strain 

38 

E x 

electric field 

2 

£ 22 

total normal strain 

39 

A 

electric field 

3 

£ 33 

total normal strain 

40 

A 

electric field 

4 

A 

total engineering shear strain 

41 

A 

magnetic field 

5 

r is 

total engineering shear strain 

42 

H 2 

magnetic field 

6 

r 12 

total engineering shear strain 

43 

A 

magnetic field 

7 

Ai 

stress 

44 

A 

electric displacement 

8 

G 22 

stress 

45 

A 

electric displacement 

9 

G 33 

stress 

46 

A 

electric displacement 

10 

G 23 

stress 

47 

B x 

magnetic flux density 

11 

A 3 

stress 

48 

A 

magnetic flux density 

12 

G \2 

stress 

49 

A 

magnetic flux density 

13 

F in 

6 11 

inelastic normal strain 

50 

w 

thermal electric field 

14 

F in 

^22 

inelastic normal strain 

51 

E t 2 

thermal electric field 

15 


inelastic normal strain 

52 

El 

thermal electric field 

16 

A” 

inelastic engr. shear strain 

53 

Hi 

thermal magnetic field 

17 

Yu 

inelastic engr. shear strain 

54 

Hi 

thermal magnetic field 

18 

n 2 

inelastic engr. shear strain 

55 

Hi 

thermal magnetic field 

19-30 

N/A 

56-99 

N/A 

31 

s T 

*11 

thermal normal strain 

100 

T 

temperature 

32 

S T 

*22 

thermal normal strain 

101 

t 

time 

33 

s T 

*33 

thermal normal strain 

102 

^ creep 

creep time 

34 

-T 

Yu 

thermal engineering shear strain 




35 

-T 

r is 

thermal engineering shear strain 




36 

-T 

Yu 

thermal engineering shear strain 




37 

N/A 

N/A 
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*XYPLOT : X-Y plot Output 


Table 6.4 MICRO level x-y plot data i.d. numbers. 


1 

£ \\ 

total normal strain 

34 

T 

r 2 3 

thermal engineering shear strain 

2 

S 22 

total normal strain 

35 

T 
Y 13 

thermal engineering shear strain 

3 

£ 33 

total normal strain 

36 

rh 

thermal engineering shear strain 

4 

Yl 3 

total engineering shear strain 

37 

K 

normal debond parameter 

5 

Y 13 

total engineering shear strain 

38 

K 

normal debond parameter 

6 

Y 12 

total engineering shear strain 

39 

R l 

normal debond parameter 

7 

Al 

stress 

40 

R] 

shear debond parameter 

8 

G 22 

stress 

41 

Rf 

shear debond parameter 

9 

G 33 

stress 

42 

R< 

shear debond parameter 

10 

G 23 

stress 

43 

e , 

electric field 

11 

G \3 

stress 

44 

E 2 

electric field 

12 

G \2 

stress 

45 

E 3 

electric field 

13 

F in 

6 11 

inelastic normal strain 

46 

H , 

magnetic field 

14 

p in 

b 22 

inelastic normal strain 

47 

H 2 

magnetic field 

15 

4 

inelastic normal strain 

48 

H 3 

magnetic field 

16 

4 

inelastic engr. shear strain 

49 

A 

electric displacement 

17 

yS 

inelastic engr. shear strain 

50 

A 

electric displacement 

18 

yZ 

inelastic engr. shear strain 

51 

A 

electric displacement 

19 

On 

inelastic model state variable 

52 

B \ 

magnetic flux density 

20 

o 22 

inelastic model state variable 

53 

A 

magnetic flux density 

21 

o 33 

inelastic model state variable 

54 

B 3 

magnetic flux density 

22 

$23 

inelastic model state variable 

55 

El 

thermal electric field 

23 

O n 

inelastic model state variable 

56 

El 

thermal electric field 

24 

Oi 2 

inelastic model state variable 

57 

El 

thermal electric field 

25 


inelastic model state variable 

58 

Hi 

thermal magnetic field 

26 

%2 

inelastic model state variable 

59 

Hi 

thermal magnetic field 

27 

^33 

inelastic model state variable 

60 

Hi 

thermal magnetic field 

28 

^23 

inelastic model state variable 

61-99 

N/A 

29 

%3 

inelastic model state variable 

100 

T 

temperature 

30 

%2 

inelastic model state variable 

101 

t 

time 

31 

S T 

b \\ 

thermal normal strain 

102 

^ creep 

creep time 

32 

s T 

*22 

thermal normal strain 




33 

e T 

*33 

thermal normal strain 
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* PATRAN : PATRAN post-processing data output 


Purpose: 

Generate data output files for post-processing using the MSC/PATRAN software in conjunction with the 
MACPOST software. The MACPOST software is an add-on to the MSC/PATRAN software package 
that enables generation of x-y plots and fringe plots. The MACPOST software add-on is distributed with 
MAC/GMC 4.0 and has its own User Manual (Goldberg et al . , 1999). The user is referred to the 
MACPOST User Manual for specific instructions on post-processing the data generated from MAC/GMC 
4.0. Note that the MSC/PATRAN output is available only for doubly and triply periodic GMC RUC 
analyses. 


Usage: 

* PATRAN 

TPRE =tpre STP =stp 

Name Definition 

tpre Preloading time before which MSC/PATRAN output data files are not written 

stp 1 Frequency for writing MSC/PATRAN output, output is written every stp time step 


Example: 

* PATRAN 

TPRE=57 600 . STP=50 


Associated MAC/GMC 4.0 Example Problemfst 

Example 6b 


Notes: 


1) A small value of stp will cause a great deal of data to be written and may slow execution 
considerably. 


2) The * PATRAN keyword causes the generation of the following 14 data files, where outfile is the 
MAC/GMC 4.0 output file name: 


outfile . macgeo 
outfile . total_pat . data 
outfile . macro l_pat . data 
outfile . macro 2_pat . data 


RUC geometry data 
Number of time steps and subcells 
RUC level strain data at each output time 
RUC level stress data at each output time 
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*PATRAN: PATRAN post-processing data output 


outfile. macro 3_pat . data 
out file . macro 4_pat . data 

outfile . microl_pat . data 
outfile . micro2_pat . data 
outfile . micro3_pat . data 
outfile . micro4_pat . data 


outfile . microl_pat . contour 
outfile . micro2_pat . contour 
outfile . micro3_pat . contour 
outfile . micro4_pat . contour 


RUC level inelastic strain data at each output time 

RUC level thermal strain, creep time, temperature, and stress 

invariant data at each output time 

Subcell level strain data at each output time 

Subcell level stress data at each output time 

Subcell level inelastic strain data at each output time 

Subcell level thermal strain, creep time, temperature, and stress 

invariant data at each output time 

Subcell level strain data for fringe plots 

Subcell level stress data for fringe plots 

Subcell level inelastic strain data for fringe plots 

Subcell level thermal strain, creep time, temperature, and stress 

invariant data for fringe plots 
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* MAT LAB : MATLAB post-processing data output 


Purpose: 

Generate data output files for post-processing using the MATLAB software in conjunction with the 
MATLAB source files distributed with MAC/GMC 4.0. The MATLAB output is available only for 
doubly periodic RUC analysis and doubly periodic RUC layers within a laminate analysis. 


Usage: 

*MATLAB 

N =ntimes TIMES =t(l),t(2), ...,t(ntimes) 

Name Definition 

ntimes Number of times at which MATLAB output will be written 

t(l-ntimes) The times at which MATLAB output will be written 


Example: 

*MATLAB 

N=3 TIMES=57600 57700 57800 . 


Associated MAC/GMC 4.0 Example Problemfst 

Example 3f, Example 6c 


Notes: 


1) The * MATLAB keyword causes the generation of the following 5 data files, where outfile is the 
MAC/GMC 4.0 output file name: 


outfile_x 2 . dat 
outfile_x 3 . dat 
outfile_s±q . dat 
outfile_e ps . dat 
outfile_e psp . dat 


x 2 subcell coordinate geometry information 
x 3 subcell coordinate geometry information 
Subcell level stress data at each output time 
Subcell level strain data at each output time 
Subcell level inelastic strain data at each output time 


To be read by the three MATLAB source files (see Note 2 below), the leading 66 outfile must be 
stripped from the outfile_x 2.dat and outfilejx 3.dat files. The remaining three MATLAB 
data files must contain only the data for a single time, without the time itself included. Thus, the user 
must copy the desired data out of the outflle_ sig.dat, outfile_e ps.dat, and outfile- 
_epsp . dat files and copy it into the corresponding file with the leading 66 outfile stripped from 
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Section 6: Results and Data Output Keywords 


*matlab : MATLAB post-processing data output 


the name. The MATLAB source files can then correctly read these three data files and create the 
appropriate fringe plots. 

2) The following MATLAB source files, which generate fringe plots from the data contained in the 
MAC/GMC 4.0 MATLAB output files, are distributed with MAC/GMC 4.0: 

stress . m Generates stress component and invariant fringe plots 

s train . m Generates strain component fringe plots 

epsp . m Generates inelastic strain component fringe plots 

Each file contains an explanation and a section for user edits that control some of the characteristics 
of the fringe plots. These MATLAB source files are printed in the Appendix. 

3) The fringe plots can be generated by placing the MATLAB source files and the MAC/GMC 4.0 
MATLAB output files in the MATLAB working directory, and then executing the MATLAB source 
files by typing the file name at the MATLAB command line. For example, typing “stress” at the 
MATLAB command line would execute the stress.m program to generate the stress component 
fringe plots. 
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Section 7 : End of File Keyword 


This section contains the keyword that is required at the end of every MAC/GMC 4.0 input file: 
*end — » Indicate that the end of the MAC/GMC 4.0 input file has been reached 
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Section 7: End of File Keyword 


*end : End of input file 


*END : End of input file 


^ *END is a required keyword R 5 ' 

Purpose: 

Define the end of the MAC/GMC 4.0 input file. 


Usage: 

*END 


Example: 

*END 


Associated MAC/GMC 4.0 Example Problem(s): 

All example problems 


Notes: 

1) Lines present in the input file that are after the *END keyword will not be read by the code. The user 
can thus hold particular lines of input file data by placing them after this keyword at the end of the 
input file. 
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Appendix 


Appendix 


This Appendix contains the following: 

• User-Defined Constituent Material Constitutive Model File 

• User-Defined Stiffness Matrix Determination File 

• User-Defined Material Property Determination File 

• Sample external material database file 

• The three MATLAB source files for generating fringe plots 

—> stress . m, strain . m, and epsp . m 
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User-Defined Constituent Material Constitutive Model File 


usrmat . F90 


!####################################################################### 
SUBROUTINE USRMAT (DSA, SA, PE , PV, D, LOCTISO, TIME , TSTEP, & 
CTEMP, DTEMPR, NIO, NE , NV, NS , MN, DMGF, NSASIZE) 


Purpose: MAC/GMC 4.0 user-defined material constitutive model for 

determination of the inelastic strain and state variable 
increments 

This subroutine is called when one or more constituent 
materials are described by user-defined constitutive 
models, that is, when: 

CMOD =99 

Reference: MAC/GMC 4.0 Keywords Manual Section 2, MAC/GMC 4.0 

Example Manual (Example 2g) 

Called From: FORMIN 
Subs Called: FATALERROR 


* On Entry 
SA 

PE (NE ) 

PV (NV) 

D(3) 

LOCTISO -- 


TIME 

TSTEP 

CTEMP 

DTEMPR -- 

NIO 

NE 

NV 

NS 

MN 

DMGF 


* 

Vector of total (integrated) quantities (see below) 

Vector of elastic constants for material MN 
(where NE = # of elastic constants --> 9 MAX) 

Vector of viscoplastic constants for material MN 
(where NV = # of viscoplastic constants --> 19 MAX) 

Vector of direction cosines (for models 3, 7, & 9) 

Flag indicating if ANY material exhibits local 
transverse isotropy (and global anisotropy) 

= 0 - all materials are at most globally transversely 
isotropic (D not used) 

= 1 - at least one material is locally transversely 
isotropic (D used) 

Current time 

Current time step 

Current temperature 

Time rate of change of temperature 

Unit number of .out file 

# of elastic constants --> 9 MAX 

# of viscoplastic constants --> 19 MAX 

Current subcell number (this routine is called from within 
a subcell loop) 

Material number 

Damage factor - if damage is included the user 
should multiply material stiffness terms by DMGF 
when using such terms in his inelastic model. 


* Expected On Exit * 

DSA -- Vector of quantity increments (see below) 


NOTE: In this subroutine, [SA] and [DSA] contain the 

micro (subcell) quantities for the micromechanics model 

Arrangement of [DSA] & [SA] arrays: 


variable location 


+ 

I total strain 

(1-6) (contains ENGINEERING shears) 

1 

| stress 

i 

(7-12) 

1 

| inelastic strain 

(13-18) (contains ENGINEERING shears) 
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Appendix 


User-Defined Constituent Material Constitutive Model File 


1 

12 "slots" 

(19-30) 



1 

for state variables 




1 

thermal strain 

(31-36) 

(contains ENGINEERING shears) 

1 

Rn & Rt debond 

(37 - Rn 

11 

face) 

1 

parameters 

a 
c L 

00 

00 

22 

face) 

1 


(39 - Rn 

33 

face) 

1 


o 

i 

W 

r+ 

11 

face) 

1 


-P 

ft 

\ — 1 

22 

face) 

1 


(42 - Rt 

33 

face) 

+- 







NOTE: The user is responsible for providing the inelastic strain and 
state variable increments, DSA slots (13-30) . The rest of DSA 
should not be disturbed 

NOTE: To use your own code within this subroutine, this subroutine 

must be compiled and linked with the MAC/GMC 4.0 library file 
mac4 . lib 

** IT IS THE USER'S RESPONSIBILITY TO ENSURE THAT THE USER'S CODE ** 
** FUNCTIONS PROPERLY AND THAT THE DESIRED CODE IS ACTUALLY BEING ** 
** EXECUTED ** 


IMPLICIT DOUBLE PRECISION (A - H, 0 - Z) 


DIMENSION SS ( 6 ) , S(6), R(6) 

DIMENSION DSA(NSASIZE) , SA(NSASIZE) 
DIMENSION PV(NV), PE(NE), D(3) 

LOGICAL LOCTISO 


NEP = 9 
NVP =19 

IF ( NE . GT. NEP ) THEN 
CALL FATALERROR (NIO) 

WRITE (NIO, *) 'TO MANY ELASTIC PROPERTIES FOR MATERIAL # ' , MN 
WRITE (NIO, *) ' # USED = ' , NE 

WRITE (NIO, *) ' # ALLOCATED = ', NEP 

STOP 
ENDIF 

IF ( NV .GT. NVP ) THEN 
CALL FATALERROR (NIO) 

WRITE (NIO, *) 'TO MANY NON-LINEAR PROPERTIES FOR MATERIAL # ', & 

MN 

WRITE (NIO, *) ' # USED = ' , NV 

WRITE (NIO, *) ' # ALLOCATED = ', NVP 

STOP 
ENDIF 


BEGIN USER EDITS * 


! -- Consider placing your own message here to ensure that your code is 
! being executed 

IF (TIME .EQ. 0 . DO ) THEN 

WRITE (NIO, *) '*** NOTE: Default USRMAT routine being executed ***' 

WRITE (*, *) '*** NOTE: Default USRMAT routine being executed ***' 

ENDIF 
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MATERIAL #1 -- Power Law Creep 


IF (MN .EQ. 1) THEN 


! -- Creep material properties 
CREEPCOEF = PV(1) 

POWER = PV (2 ) 

! -- Deviatoric stress [S] in the subcell 

TRACE3 = (SA(7) + SA(8) + SA(9)) / 3. 
S (1) = SA ( 7 ) - TRACE3 
S (2) = SA(8) - TRACE3 
S (3) = SA ( 9 ) - TRACE3 
S (4) = SA(10) 

S (5) = SA (11) 

S ( 6 ) = SA ( 12 ) 


Effective Stress 


SEFF = SQRT ( ( SA ( 7 ) - SA(8))**2 + (SA(8) 

( SA ( 7 ) - SA ( 9 ) ) **2 + & 

+ 6 * (SA (10) **2 + SA ( 11 ) **2 


SA ( 9 ) ) **2 + & 

+ SA ( 12 ) **2 ) ) / SQRT (2 . 


CREEP__F = CREEPCOEF * SEFF**POWER 

IF (SEFF .GT. 0 . DO ) THEN 

COEF = 3 * CREEP_F / (2 * SEFF) 

ELSE 

COEF = 0.D0 
ENDIF 

! -- Inelastic strain increments 


DSA ( 13 ) 

= COEF 

* 

S(l) 

* 

TSTEP 

DSA (14) 

= COEF 

* 

S (2 ) 

* 

TSTEP 

DSA ( 15 ) 

= COEF 

* 

S (3) 

* 

TSTEP 

DSA (16) 

= COEF 

* 

S ( 4 ) 

* 

TSTEP 

DSA (17) 

= COEF 

* 

S (5) 

* 

TSTEP 

DSA ( 18 ) 

= COEF 

* 

S ( 6) 

* 

TSTEP 


MATERIAL #2 — > Bodner-Partom 


ELSEIF (MN .EQ. 2) THEN 


! -- Copy appropriate viscoplastic material constants 


IF ( NV . LT. 6 ) THEN 
CALL FATALERROR (NIO) 

WRITE (NIO, *) ' NOT ENOUGH PV SPACE: NV = ', NV 

STOP 
ENDIF 


DO = PV(1) 
Z0 = PV (2 ) 
Z1 = PV ( 3 ) 
BM = PV ( 4 ) 
AN = PV ( 5 ) 
Q = PV ( 6 ) 


! -- Copy stress 
SS(1) = 
SS (2 ) = 
SS (3) = 
SS ( 4 ) = 
SS ( 5 ) = 
SS ( 6 ) = 


from [SA] to [SS] 
SA ( 7 ) 

SA ( 8 ) 

SA ( 9) 

SA(10) 

SA (11) 

SA ( 12 ) 


! -- Compute the deviatoric stress [S] in the subcell 
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User-Defined Constituent Material Constitutive Model File 


TEMP = (SS(1) + SS(2) + SS ( 3 ) ) / 3.0 
S (1) = SS (1) - TEMP 

S (2) = SS (2) - TEMP 
S (3) = SS (3) - TEMP 
S (4) = SS (4) 

S (5) = SS ( 5 ) 

S (6) = SS (6) 

AJ2 = 0.5 * (S(l)**2 + S (2 ) **2 + S(3)**2) + S(4)**2 + S (5) **2 + & 
S (6) **2 

SQ3AJ = DSQRT ( SS(1)**2 + SS(2)**2 + SS(3)**2 + 2 * (SS (4) **2 + & 
SS (5) **2 + SS (6) **2 ) ) 

SQ2 = 1.414215 

IF (SQ3AJ .EQ. 0.0) THEN 
DO I = 1, 6 

R ( I ) = 0 
END DO 

ELSE 

R ( 1 ) = SS (1) / SQ3AJ 

R (2 ) = SS (2) / SQ3AJ 
R ( 3 ) = SS (3) / SQ3AJ 

R ( 4 ) = SQ2 * SS ( 4 ) / SQ3AJ 

R ( 5 ) = SQ2 * SS (5) / SQ3AJ 

R ( 6 ) = SQ2 * SS (6) / SQ3AJ 

ENDIF 

! -- If D0=0 then assume elastic and zero-out [DSA(13-30)] (inelastic 
! strain increments and internal variable increments) , then return 
IF (DO .EQ. 0) THEN 
DO JJ = 13, 30 
DSA(JJ) =0.0 
END DO 
RETURN 

! -- Inelastic 
ELSE 

ZEF = Z0 + Q * SA (2 0 ) + (1 - Q) * (R(l) * SA(21) + R(2) * & 

SA (22 ) + R ( 3 ) * SA (23 ) + R(4) * SA(24) + R(5) & 

* SA (25 ) + R(6) * SA (2 6) ) 

IF (AJ2 .EQ. 0.0) THEN 
AL = 0.0 

ELSE 

ARG1 = ZEF* *2 . 0 / (3.0 * AJ2 ) 

IF (ARG1 .GT. 1E6) ARG1 = 1E6 
CON = .5 * (AN +1.0) /AN 
ARG = CON * (ARG1 ) **AN 

IF (ARG .GT. 50.0) ARG =50.0 
AL = DO / (DEXP (ARG) * DSQRT (AJ2)) 

ENDIF 

! -- Inelastic strain increments 


DSA ( 13 ) 

= AL 

* S (1) 

* TSTEP 

DSA (14) 

= AL 

* S(2) 

* TSTEP 

DSA ( 15 ) 

= AL 

* S(3) 

* TSTEP 

DSA (16) 

= 2 * 

AL * 

S (4) * TSTEP 

DSA (17) 

= 2 * 

AL * 

S ( 5 ) * TSTEP 

DSA (18) 

= 2 * 

AL * 

S ( 6 ) * TSTEP 


! -- Plastic work increment 

WPD = S(l) * DSA ( 13 ) + S (2 ) * DSA(14) + S(3) * DSA(15) + S(4) * & 
DSA (16) + S (5) * DSA ( 17 ) + S ( 6) * DSA(18) 

! -- State variable increments 
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DSA (19) = WPD 
ZOM = BM / ZO 

ZD = ZOM * (Z1 - ZEF) * WPD 

DSA (20 ) = ZD 

DSA (21 ) = ZD * R ( 1 ) 

DSA (22 ) = ZD * R (2 ) 

DSA (23 ) = ZD * R ( 3 ) 

DSA (24 ) = ZD * R ( 4 ) 

DSA (25 ) = ZD * R ( 5 ) 

DSA (2 6) = ZD * R ( 6 ) 

ENDIF 

ENDIF 



RETURN 

END 

EOF **' 
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User-Defined Stiffness Matrix Determination File 


User-Defined Stiffness Matrix Determination File 

usrformde . F90 


!####################################################################### 
SUBROUTINE USRFORMDE (NIO, MN, TIME , PEM, PVM, D, D_MODEL, NE, NV, & 
NEP, NVP, NMTS , DNEW) 


! Purpose: MAC/GMC 4.0 user-defined subroutine to allow formation of 

! constituent material stiffness matrices based on a user 

! constitutive model, which may use arbitrary material 

! properties (used when ncmd = 99) 

! This subroutine is called when one or more constituent 

! materials are described by user-defined constitutive 

! models (using the USRMAT subroutine), that is, when: 

! CMOD =99 

! Reference: MAC/GMC 4.0 Keywords Manual Section 2, MAC/GMC 4.0 

! Example Manual (Example 2g) 

! Called From: GETCON 
! Subs Called: 


! ** On Entry ** 

! NIO -- Unit number of .out file 

! MN -- Material number 

! TIME -- Current time 

! PEM (NE , MN) -- Vector of elastic constants for material # MN 

! P VM ( NV , MN) -- Vector of viscoplastic constants for material # MN 

! D(3) -- Vector of direction cosines 

! (required for material models 3, 7, 9, & 20) 

! D_MODEL -- Flag indicating if material # MN exhibits local 

! transverse isotropy with global anisotropy 

! = .FALSE. - material is at most globally transversely 

! isotropic (D not used) 

! = .TRUE. - material is locally transversely 

! isotropic (D used) 

! NE -- # of elastic constants --> 9 MAX 

! NV -- # of viscoplastic constants --> 19 MAX 

! NEP -- Max # of elastic constants (to dimension arrays) 

! NVP -- Max # of viscoplastic constants (to dimension arrays) 

! NMTS -- Total number of materials 

! ** Expected On Exit ** 

! DNEW (6, 6) -- Current elastic material stiffness matrix 


! NOTE: In the examples shown here, standard engineering material 

! elastic constants (E, nu) are used. However, the user is 

! free to employ any material elastic constants with his 

! constitutive model ( i . e . , bulk modulus, etc.) . Thus, the 

! user must provide the equations to determine the stiffness 

! components required by GMC from his elastic constants. 

! NOTE: To use your own code within this subroutine, this subroutine 

! must be compiled and linked with the MAC/GMC 4.0 library file 

! mac4 . lib 

! ** IT IS THE USER'S RESPONSIBILITY TO ENSURE THAT THE USER'S CODE ** 

! ** FUNCTIONS PROPERLY AND THAT THE DESIRED CODE IS ACTUALLY BEING ** 

! ** EXECUTED ** 


IMPLICIT DOUBLE PRECISION (A - H, 0 - Z) 
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DIMENSION DNEW ( 6 , 6) 

DIMENSION PEM (NEP, NMTS ) , PVM(NVP, NMTS ) 
DIMENSION D ( 3 ) 

LOGICAL D MODEL 


! * BEGIN USER EDITS 


! -- Consider placing your own message here to ensure that your code is 
! being executed 

IF (TIME .EQ. 0 . DO ) THEN 

WRITE (NIO, *) '*** NOTE: Default USRFORMDE routine being executed ***' 
WRITE (NIO, *) ' Material number', MN 

WRITE (*, *) '*** NOTE: Default USRFORMDE routine being executed ***' 
WRITE (*, *) ' Material number', MN 

ENDIF 


! -- Obtain material properties from storage 
E = PEM ( 1 , MN) 

FN = PEM (2, MN) 

GA = E / (2 . * (1 . + FN) ) 

! -- Zero the stiffness matrix 
DO I = 1, 6 

DO J = 1, 6 

DNEW (I, J) = 0 
END DO 
END DO 


! -- Determine the stiffness components 
GT = 0 . 5 * E / (1 + FN) 

FK = 0.25 * E / (0.5 * (1 - FN) - FN**2) 


DNEW ( 1 , 

1) 

= E + 

4.0 * FK ' 

DNEW ( 2 , 

1) 

= 2.0 

* FK * FN 

DNEW ( 3 , 

1) 

= 2.0 

* FK * FN 

DNEW ( 1 , 

2) 

= 2.0 

* FK * FN 

DNEW (2, 

2) 

II 

+ 

GT 

DNEW ( 3 , 

2) 

= FK - 

GT 

DNEW ( 1 , 

3) 

= 2.0 

* FK * FN 

DNEW (2, 

3) 

= FK - 

GT 

DNEW ( 3 , 

3) 

= FK + 

GT 

DNEW ( 4 , 

4) 

= GT 


DNEW (5, 

5) 

= GA 


DNEW ( 6 , 

6) 

= GA 




RETURN 

END 


j *** poF *** 
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User-Defined Material Property Determination File 


User-Defined Material Property Determination File 

usrfun. F90 


!####################################################################### 
SUBROUTINE USRFUN (NIO, MN, TIME , TSTEP, CTEMP, DTEMPR, SA, DSA, & 
DOLD, PEM, PVM, D, D_MODEL, ALPA, ALPT, & 

NE , NV, NMTS , NEP, NVP, NSASIZE, SU1, SU2 , Nl, N2 , & 

ZDAM, N Z DAM , XFAILT, XFAILC, NCRIT, CYCTOT) 


Purpose: MAC/GMC 4.0 user-defined subroutine to allow elastic and 

viscoplastic material properties to be functions of temp 
or field variables. 

This subroutine is called when one or more constituent 
materials employ user-defined functional form material 
properties, that is, when: 

(MATID .eq. 'U') .and. (MATDB = 2) 

Reference: MAC/GMC 4.0 Keywords Manual Section 2, MAC/GMC 4.0 

Example Manual (Example 2f) 

Called From: GETCON 
Subs Called: 


** On Entry ** 

NIO -- Unit number of MAC/GMC 4.0 output (.out) file 

MN -- Constituent material number 

TIME -- Current time 

TSTEP -- Current time step size 

CTEMP -- Current temperature 

DTEMPR -- Time rate of change of temperature 

SA -- Vector of total (integrated) quantities (see below) 

DSA -- Vector of quantity increments (see below) 

DOLD(6, 6) -- Previous elastic material stiffness matrix for material 

PEM (NE , MN) -- Vector of previous elastic constants for material 
# MN (where NE = # of elastic constants --> 9 MAX) 
PVM(NV, MN) -- Vector of previous viscoplastic constants for 
material # MN 

(where NV = # of viscoplastic constants --> 19 MAX) 

NMTS -- Total number of constituent materials 

NEP -- Max # of elastic constants (to dimension arrays) 

NVP -- Max # of viscoplastic constants (to dimension arrays) 

NSASIZE -- Length of the SA and DSA vectors (for dimensioning) 

NCRIT -- Number of failure criteria (for dimensioning) 

CYCTOT -- Current total number of cycles in damage analysis 

** Expected On Exit ** 

PEM(NE, MN) -- Vector of current elastic constants for material MN 
PVM(NV, MN) -- Vector of current viscoplastic constants for 
material MN 

D(3) -- Vector of direction cosines of normal to plane of 

trans. isotropy (required for models 3, 7, & 9) 

D_MODEL -- Flag indicating if material # MN exhibits local 

transverse isotropy with global anisotropy 
= .FALSE. - material is at most globally transversely 
isotropic (D not used) 

= .TRUE. - material is locally transversely 
isotropic (D used) 

ALPA(MN) -- Longitudinal cte for material MN 
ALPT (MN) -- Transverse cte for material MN 
NE -- Number of elastic properties used 

NV -- Number of viscoplastic properties used 

Strength reduction fatigue damage model 

SU1 (MN, I) -- Material damage property (component I) 
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SU2 (MN, I) 
N1 (MN, I) 
N2 (MN, I) 


Material damage property 
Material damage property 
Material damage property 


(component I) 
(component I) 
(component I) 


ZDAM(MN,20) — 


NZDAM(MN) 


Stiffness 
Material 
Z DAM ( MN , 

Z DAM ( MN , 
ZDAM (MN, 
ZDAM (MN, 

Z DAM ( MN , 

Z DAM ( MN , 
ZDAM (MN, 
ZDAM (MN, 

Z DAM ( MN , 

Z DAM ( MN , 
ZDAM (MN, 
ZDAM (MN, 

Z DAM ( MN , 

Z DAM ( MN , 
ZDAM (MN, 
ZDAM (MN, 


1 ) = ANG 

2) = BN 

3) = BP 

4 ) = OMU 

5 ) = OMFL 

6 ) = OMM 

7) = ETU 

8) = ETFL 

9) = ETM 

10) = BE 

11) = * empty * 

12) = A_hat 

13) = SFL 

14) = XML 

15) = SU 

16-20) = open space 
SK - Flag to skip damage calculations for material 
= 0 - Perform damage calculations 
= 1 - Skip damage calculations 


XFAILT 
(MN, ICRIT , I) 
XFAILC 
(MN, ICRIT, I) 


Local (subcell) failure model 

Tensile failure properties (Xi j ) , material #MN, 
criterion #ICRIT, component #1 

Compressive failure properties (XCi j ) , material #MN, 
criterion #ICRIT, component #1 


NOTE: This routine can be used in conjunction with an internal 

material constitutive model, or a constitutive model input by 
the user in the USRMAT subroutine 

NOTE: In this subroutine, [SA] and [DSA] contain the 

micro (subcell) quantities for the micromechanics model 

Arrangement of [DSA] & [SA] arrays: 


variable 

location 


total strain 

(1-6) (contains ENGINEERING shears) 

stress 

(7-12) 


inelastic strain 

(13-18) (contains ENGINEERING 

shears ) 

12 "slots" 

for state variables 

(19-30) 


thermal strain 

(31-36) (contains ENGINEERING 

shears) 

Rn & Rt debond 

(37 - Rn, 11 face) 


parameters 

(38 - Rn, 22 face) 
(39 - Rn, 33 face) 
(40 - Rt, 11 face) 
(41 - Rt, 22 face) 
(42 - Rt, 33 face) 



NOTE: quantities in [SA] and [DSA] are SUBCELL quantities - the 

values on entry are for the first subcell containing material 
# MN - the values on exit of this subroutine will be applied to 
ALL SUBCELLS containing material # MN. It is thus recommended 
that, if using the field variables, you assign the appropriate 
material # to ONE SUBCELL ONLY. Use of [SA] and [DSA] in this 
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context in conjunction with bending in laminate theory will 
result in erroneous results as field variables become dependent 
on through-thickness position (while the material # does not) . 

NOTE: To use your own code within this subroutine, this subroutine 

must be compiled and linked with the MAC/GMC 4.0 library file 
mac4 . lib 

** IT IS THE USER'S RESPONSIBILITY TO ENSURE THAT THE USER'S CODE ** 

** FUNCTIONS PROPERLY AND THAT THE DESIRED CODE IS ACTUALLY BEING ** 

** EXECUTED ** 


IMPLICIT DOUBLE PRECISION (A - H, 0 - Z) 


DIMENSION DOLD ( 6 , 6), D(3) 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


PEM (NEP, NMTS), PVM(NVP, NMTS ) 

ALPA (NMTS ) , ALPT (NMTS ) 

DSA (NSASIZE) , SA(NSASIZE) 

SU1 (NMTS , 6 ) , SU2 (NMTS , 6) 

N1 (NMTS , 6),N2(NMTS, 6) 

Z DAM ( NMT S , 20), NZDAM(NMTS) 

XFAILT (NMTS , NCRIT, 6), XFAILC (NMTS , 


NCRIT, 6) 


LOGICAL D MODEL 


BEGIN USER EDITS * 


! -- Consider placing your own message here to ensure that your code is 
! being executed 

IF (TIME .EQ. 0 . DO ) THEN 

WRITE (NIO, *) '*** NOTE: Default USRFUN routine being executed ***' 

WRITE (NIO, *) ' material number', MN 

WRITE (*, *) '*** NOTE: Default USRFUN routine being executed ***' 

WRITE (*, *) ' Material number', MN 

ENDIF 


MATERIAL # 1: Elastic, E = function of strain & previous E 


IF (MN .EQ. 1) THEN 

IF ( SA ( 1 ) .GT. 0) THEN 
EAOLD = PEM ( 1 , MN) 

EA = 700. E9 - EAOLD * SA(1) * 80.0 
ELSE 

EA = 700.E9 
EAOLD = 700.E9 
ENDIF 

IF (EA .GT. EAOLD) THEN 

WRITE (NIO, *) 'NON-PHYSICAL REGIME ENTERED' 
STOP 


ENDIF 


EA = EA 


ET = EA 
FNA = 0.41 
FNT = 0.41 
GA = EA / 
ALPA (MN) = 
ALPT (MN) = 

(2.0 * 
4.5E-6 
4.5E-6 

NE = 5 
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PEM 

(1 

f MN) 

= 

EA 

PEM 

(2 

f MN) 

= 

ET 

PEM 

(3 

, MN) 

= 

FNA 

PEM 

(4 

, MN) 

= 

FNT 

PEM 

(5 

f MN) 

= 

GA 

EA ^ 

= ! 

EA 



ET : 

= ; 

EA 



FNA 

= 

0.41 



FNT 

= 

0.41 



GA = 

= : 

EA / 

(2. 

.0 * 

ALPA(MN) = 

4 . 

. 5E- 1 

ALPT (MN) = 

4 . 

.5E-i 

NE = 

= 

5 



PEM 

(l 

, MN) 

= 

EA 

PEM 

(2 

f MN) 

= 

ET 

PEM 

(3 

, MN) 

= 

FNA 

PEM 

(4 

, MN) 

= 

FNT 

PEM 

(5 

, MN) 

= 

GA 


MATERIAL # 2 : Incremental plasticity with strain rate dependence 


ELSEIF (MN .EQ. 2) THEN 


EA = 13161. 

ET = EA 
FNA = 0.26 
FNT = 0.26 

GA = EA / (2.0* (1.0+ FNA)) 
ALPA (MN) = 9.00E-6 
ALPT (MN) = 9.00E-6 

NE = 5 

PEM ( 1 r MN) = EA 
PEM (2 , MN) = ET 
PEM (3, MN) = FNA 
PEM (4, MN) = FNT 
PEM (5, MN) = GA 


NP = 3 

NV = 2 * NP + 1 

STRAIN_RATE = 0.4714 * (DSQRT ( (DSA ( 1 ) - DSA(2))**2 + (DSA(l) - DSA(3))**2 + & 

(DSA (2 ) - DSA (3) ) **2 + 1.5 * DSA(4)**2 + & 

1.5 * DSA (5) **2 + 1.5 * DSA ( 6 ) **2 ) ) / TSTEP 


IF (STRAIN_RATE . GT . l.E-2) STRAIN_RATE = l.E-2 

IF (STRAIN_RATE . LT . l.E-7) THEN 
SY = 70. 

51 = 80. 

52 = 84. 

53 = 87. 

ELSE 


SY = 

70. 

. + 

■ 3. 

* 

LOG (1 .E7*STRAIN RATE) 

SI = 

SY 

+ 

10. 

* 

SY / 70. 

S2 = 

SY 

+ 

14 . 

* 

SY / 70. 

S3 = 
ENDIF 

SY 

+ 

17 . 

* 

SY / 70. 


E1=0 . 01 
E2=0 .015 
E3=0 . 2 

PVM ( 1 r MN) = SY 
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PVM ( 2 , 

MN) = SI 

PVM ( 3 , 

MN) = S2 

PVM ( 4 , 

MN) = S3 

PVM ( 5 / 

MN) = El 

PVM ( 6 f 

MN) = E2 

PVM ( 7 , 

MN) = E3 

MATERIAL # 3 

: properties = linear function of temperature 

for use with user constitutive model (actually b-p) 


ELSEIF (MN .EQ. 3) THEN 


IF (CTEMP . LT. 21.0) CTEMP =21.0 
IF (CTEMP .GT. 400.) CTEMP = 400. 


! -- Elastic 


E = 72 

. 4E9 

- 81. 

53E6 * 

(CTEMP - 21 . ) 

FN = 0 

.33 

+ 7.91 

6E-5 * 

(CTEMP - 21 . ) 

ALP = 

22.5E-6 + 

3. 958E- 

-9 * (CTEMP - 21 

NE = 5 





PEM ( 1 , 

MN) 

= E 



PEM (2, 

MN) 

= E 



PEM (3, 

MN) 

= FN 



PEM (4, 

MN) 

= FN 



PEM (5, 

MN) 

= E / 

(2. * 

(1 . + FN) ) 

ALPA(MN) = 

ALP 



ALPT (MN) = 

ALP 




! -- Viscoplastic (Bodner-Partom) 

DO = 1.E4 
Z0 = 340. E6 
Z1 = 435. E6 
BM = 300 . 0 

AN = 10.0 - 0.02493 * (CTEMP - 21.) 
Q =1.0 


NV = 6 


PVM ( 1 f 

MN) 

= DO 

PVM ( 2 f 

MN) 

= ZO 

PVM ( 3 , 

MN) 

= Z1 

PVM ( 4 f 

MN) 

= BM 

PVM ( 5 / 

MN) 

= AN 

PVM ( 6 / 

MN) 

= Q 


MATERIAL # 4: Fiber (NCMD = 6) with local failure = funct of temp 


ELSEIF (MN .EQ. 4) THEN 


PEM ( 1 r MN) 
PEM ( 2 , MN) 
PEM (3, MN) 
PEM (4, MN) 
PEM (5, MN) 


58 .E3 
58 .E3 
0.25 
0.25 
23.2E3 


ALPA (MN) = 3.5E-6 
ALPT (MN) = 3.5E-6 


-- Note: *SUBCELL_FAILURE must be specified in the input file to 
enable local failure checking. Any failure properties 
specified here will replace those specified in the input file 
XFAILT (MN, 1, 1) = 3500. - 1.8 * (CTEMP - 21.) 


MATERIAL #5: Matrix (NCMD = 6) with damage 


ELSEIF (MN .EQ. 5) THEN 
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PEM (1, 

MN) 

= 12770 

PEM (2, 

MN) 

= 12770 

PEM (3, 

MN) 

= 0.32 

PEM (4, 

MN) 

= 0.32 

PEM (5, 

MN) 

= PEM ( 1 


MN ) / ( 2 . * (1.+ PEM ( 3 , MN) ) ) 


ALPA(MN) = 21.06E-6 
ALPT (MN) = 21.06E-6 


-- Note: ^DAMAGE must be specified in the input file to enable fatigue 
damage analysis. Any damage properties specified here will 
replace those specified in the input file 


ZDAM (MN, 

1) 

= 

0. 

ZDAM (MN, 

2) 

= 

0. 

Z DAM ( MN , 

3) 

= 

0. 

Z DAM ( MN , 

4) 

= 

1. 

ZDAM (MN, 

5) 

= 

1. 

ZDAM (MN, 

6) 

= 

1. 

Z DAM ( MN , 

7) 

= 

1. 

Z DAM ( MN , 

8) 

= 

1. 

ZDAM (MN, 

9) 

= 

1. 

ZDAM (MN, 

10) 

= 

2.27 

Z DAM ( MN , 

11) 


= 0 

Z DAM ( MN , 

12) 


= 0.2302 

Z DAM ( MN , 

13) 


= 20.3 

ZDAM (MN, 

14) 


= 900. 

ZDAM (MN, 

15) 


= 128. 

NZDAM(MN) 

= 

0 



ENDIF 


END USER EDITS * 


RETURN 

END 


I -k -k -k EOF * * * 
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Sample external material database file 

sample_material_database .mat 


Sample Material Database Provided with MAC/GMC 4.0 

# 

# Bodner-Partom Model Material 

MAT=copper 

CMOD=l 

EL=120.E9, 120.E9, 0.33, 0.33, 45 . 11E9, 14 . 7E-6, 14 . 7E-6 & 
VI^1.E4, 63. E6, 2 50. E6, 8.19,7.5, 0.55 

# 

# Isotropic, Elastic Materials 

MAT=SCS- 6 

CMOD=6 

NTP=5 

TEM=2 1.0,316.0,427 .0,538.0,860.0 

EA=3 93 . E 9 , 382 . E 9 , 378. E9, 374 . E9 , 368. E9 

ET=3 93 . E 9 , 382 . E 9 , 378 . E 9 , 374 . E 9 , 368 . E9 

NUA=0. 25, 0.25, 0.25, 0.25, 0.25 

NUT=0. 25, 0.25, 0.25, 0.25, 0.25 

GA=157.2E9, 152.8E9, 151.2E9, 149.6E9,147.2E9 

ALPA=3 . 56E-6, 3 . 73E-6, 3 . 91E-6, 4 . 07E-6, 4 . 57E-6 

ALPT=3 . 56E-6, 3 . 73E-6, 3 . 91E-6, 4 . 07E-6, 4 . 57E-6 

# 

# User constitutive model example 

MAT=A1 (6061-0) 

CMOD=9 9 
NPE=2 NPV=6 
NTP=6 

TEM=21 . 0, 148 . 9, 204. 4, 2 60. 0,37 1.1, 400.0 
El=7 2 . 5E 9 , 69.4E9, 65. 8E9, 58.5E9, 41 . 5E9, 41 .5E9 
E2=0. 33, 0.33, 0.33, 0.33, 0.33, 0.33 

ALPA=22 . 5E-6, 22 . 5E-6, 22 . 5E-6, 22 . 5E-6, 22 . 5E-6, 22 . 5E-6 
ALPT=22 . 5E-6, 22 . 5E-6, 22 . 5E-6, 22 . 5E-6, 22 . 5E-6, 22 . 5E-6 
V1=1.E4, 1.E4, 1.E4, 1.E4, 1.E4, 1.E4 
V2=100 . E6, 100. E6, 100. E6, 100. E6, 100. E6, 100. E6 
V3=190 . E6, 190. E6, 190. E6, 190. E6, 190. E6, 190. E6 
V4=70. 0,70. 0,70. 0,70. 0,70. 0,70.0 
V5=10. 0,7. 0,4. 0,1. 6, 0.55, 0.55 
V6= 0.4, 0.4, 0.4, 0.4, 0.4, 0.4 
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stress .m 


stress. m MATLAB source file 


% Distributed with MAC/GMC 4.0 

% This program will generate fringe plots with the MATLAB software 
% from the data written to ASCII files by MAC/GMC 4.0 using the 
% * MAT LAB option 

% This program reads the file sig.dat, which contains 8 columns: 

% Sll S22 S33 S23 S13 S12 Smean Seff 

% NOTE: Smean = (Sll + S22 + S33) / 3. 

% Seff = sqrt (3*si j *si j /2 ) (sij = deviatoric stresses) 

% = [ 1/ sqrt (2 ) ] * sqrt [ (Sll - S22)**2 + (S22 - S33)**2 + 

% (Sll - S33) **2 + 6*(S23**2 + S13**2 + S12**2)] 

% While the ASCII output files from MAC/GMC 4.0 contain data for 
% all specified times (and integration points for laminate) , the 

% sig.dat file should contain just the stress field data for the 

% one particular time (and integration point) at which the fringe 

% plots will be made. 

% In addition, this program reads the files x2.dat and x3.dat, 

% which are also generated by MAC/GMC 4.0 to define the subcell 
% grid. Note that this grid will differ between GMC and HFGMC . 


% -- Clear workspace 
clear; 


USER OPTIONS 


% -- Select stress components to plot (1 - plot, 0 - don't plot) 

n s 1 1 = 1 ; 

n s 2 2 = 1 ; 

n s 3 3 = 1 ; 

n s 2 3 = 1 ; 

n s 1 3 = 1 ; 

nsl2 = 1; 

nsm = 1; 

nse = 1; 

% -- Select subplot (all on one Fig.) or individual Figs. 

% nsub = 1 --> subplot 

% nsub = 0 --> individual Figs, 

nsub = 0; 


% -- Set scale_spec = 1 to specify scale high and low values for each 
% component - set values below 
scale_spec = 0; 

if scale_spec > 0 

sll_low = -200; 
sll~high = 500; 

s22_low = -200; 
s22_high = 500; 
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s33 

low = 

-200 

s33 

high 

= 500 

s23 

low = 

-200 

s23 

high 

= 500 

sl3 

low = 

-200 

sl3 

high 

= 500 

sl2 

low = 

-200 

sl2 

high 

= 500 

sm 

low = 

-200; 

sm 

high 

500; 

se 

1 ow = 

-200; 

se 

high = 

500; 


end 


% -- Set average = 1 to perform averaging of tractions across subcell 
% boundaries. This reproduces the average sense imposition of the 
% traction continuity conditions in HFGMC . 

% -- Set average = 0 for no averaging. Use average = 0 for GMC . 
average = 0; 

% -- Set nint = number of integration points in HFGMC when average = 1 
if average > 0 
nint = 11; 

end 


% -- Note: To draw subcell boundaries, change "shading interp" 
% to "shading faceted" (GMC only) 


END USER OPTIONS 


% -- Load xl.dat and x2.dat files and determine their sizes 
load x2.dat; 
m = length (x2 ) ; 

load x3.dat; 
n = length (x3 ) ; 

% -- Create grid 

[X, Y] = meshgrid (x3 , x2 ) ; 


% -- Load sig.dat 
load sig.dat; 

f=l ; 

for i=l:m; 
for j=l:n; 

sigll (i, j ) ^ 
sig22 (i, j ) = 
sig33 (i, j ) = 
sig23 (i, j ) = 
sigl3 (i, j ) = 
sigl2 (i> j ) = 
sigmean (i, j ) 
sigef f (i, j ) 
f=f +1 ; 

end 

end 


file 


sig (f , 1) ; 
sig (f , 2) ; 
sig (f , 3) ; 
sig (f , 4 ) ; 
sig (f , 5) ; 
sig (f , 6) ; 

= sig (f , 7 ) 
= sig (f , 8) ; 


% Section to average boundary tractions 
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if average > 0 


Nbeta =length (x2 ) /nint ; 
for j=l:n 

for i=l : Nbeta-1 ; 
sig22 (nint*i , j ) = 
sig22 (nint*i+l , j ) 
sig23 (nint*i, j ) = 
sig23 (nint*i+l, j ) 
end 

end 


0.5* ( sig22 (nint*i , j ) 
= sig22 (nint*i , j ) ; 
0.5* (sig23 (nint*i, j ) 
= sig23 (nint*i , j ) ; 


+ sig22 (nint*i+l , j ) ) ; 
+ sig23 (nint*i+l , j ) ) ; 


Ngama = length (x3 ) /nint ; 
for 1=1 :m 

for j=l:Ngama-l; 

sig33 (i, nint* j ) = 0 . 5* (sig33 (i, nint* j ) 
sig33 (i, nint*j+l) = sig33(i,nint*j); 
sig23 (i, nint* j ) = 0 . 5* (sig23 (i, nint* j ) 
sig23 (i, nint* j+1) = sig23 (i, nint* j ) ; 

end 

end 


+ sig33 (i, nint* j+1) ) 
+ sig23 (i, nint* j+1) ) 


end 


if nsub > 0 

num_plots = nsll + ns22 + ns33 + ns23 + nsl2 + nsl3 + nsm + nse; 


if num 

plots == 1 


nl 

"= 2; 


n2 

= 2; 


elseif 

num plots — 

2 

nl 

2; 


n2 

= 2; 


elseif 

num plots — 

3 

nl 

2; 


n2 

= 2; 


elseif 

num plots — 

4 

nl 

- 2; 


n2 

= 2; 


elseif 

num plots — 

5 

nl 

3; 


n2 

= 3; 


elseif 

num plots — 

6 

nl 

3; 


n2 

= 3; 


elseif 

num plots — 

7 

nl 

3; 


n2 

= 3; 


elseif 

num plots — 

8 

nl 

3; 


n2 

= 3; 



end 


end 


num = 0; 

if nsub > 0 
figure; 
hold on 

end 


% -- Sll plot 
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if nsll > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sigll) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [ sll_low sll_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \ sigma \ font size { 14 } _1_1 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- S22 plot 


if ns22 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sig22 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ([ s22_low s22_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \ sigma \ font size { 14 }_2_2 ' ) ; 

xlabel ( '\bfx_3' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- S33 plot 


if ns33 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 
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colormap (jet) ; 

pcolor (X, Y, sig33 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ([ s33_low s33_high] ) ; 

end 

colorbar; 

title ( ' \bf \ fonts ize { 24 } \sigma\ font size { 14 }_3_3 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- S23 plot 


if ns23 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sig23 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ([ s23_low s23_high] ) ; 

end 

colorbar; 

title ( ' \bf\ fonts ize { 24 } \sigma\ font size { 14 }_2_3 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bf x_2 ' , ' rotation ' , 0 ) ; 


end 


% -- S13 plot 


if nsl3 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sigl3 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ([ sl3_low sl3_high] ) ; 

end 

colorbar; 
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title ( ' \bf \ fonts ize { 24 } \sigma\ font size { 14 }_1_3 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bf x_2 ' , ' rotation ' , 0 ) ; 


end 


% -- S12 plot 


if nsl2 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sigl2 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [ sl2_low sl2_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \ sigma \ font size { 14 }_1_2 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- Smean plot 


if nsm > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sigmean ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [sm_low sm_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \ sigma \ font size { 14 }_m_e_a_n ' ) ; 
xlabel (' \bfx_3' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- Seff plot 


if nse > 0 
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num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, sigef f ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [se_low se_high] ) ; 

end 

colorbar ; 

title ( ' \bf \ fonts ize { 24 } \sigma\ font size { 14 }_e_f_f ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ’ rotation ' , 0 ) ; 


end 


clear; 

break 


% -- End of file 


strain .m 


strain. m MATLAB source file 


% Distributed with MAC/GMC 4.0 

% This program will generate fringe plots with the MATLAB software 
% from the data written to ASCII files by MAC/GMC 4.0 using the 
% * MAT LAB option 

% This program reads the file epsp.dat, which contains 6 columns: 

% ell e22 e33 e23 el3 el2 

% ** NOTE: shear strains are engineering shears ** 

% While the ASCII output files from MAC/GMC 4.0 contain data for 
% all specified times (and integration points for laminate ) , the 

% eps.dat file should contain just the strain field data for the 

% one particular time (and integration point) at which the fringe 

% plots will be made. 

% In addition, this program reads the files x2.dat and x3.dat, 

% which are also generated by MAC/GMC 4.0 to define the subcell 
% grid. Note that this grid will differ between GMC and HFGMC . 


% -- Clear workspace 
clear; 


USER OPTIONS 
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% -- Select strain components to plot (1 - plot, 0 - don't plot) 

nell = X; 

ne22 = 1; 

ne33 = 1; 

ne23 = 1; 

nel3 = 1; 

nel2 = 1; 

% -- Select subplot (all on one Fig.) or individual Figs. 

% nsub = 1 --> subplot 

% nsub = 0 --> individual Figs, 

nsub = 0 

% -- Set scale_spec = 1 to specify scale high and low values for each 
% component - set values below 
scale_spec = 0; 

if scale_spec > 0 

ell_low = -0.020; 
ell_high = 0.025; 

e22_low = -0.020; 
e22_high = 0.025; 

e33_low = -0.020; 
e33_high = 0.025; 

e23_low = -0.020; 
e23_high = 0.025; 

el3_low = -0.020; 
el3_high = 0.025; 

el2_low = -0.020; 
el2_high = 0.025; 


end 


% -- Note: To draw subcell boundaries, change "shading interp" 
% to "shading faceted" (GMC only) 


END USER OPTIONS 


% -- Load xl.dat and x2.dat files and determine their sizes 
load x2.dat; 
m = length (x2 ) ; 


load x3.dat; 
n = length (x3 ) ; 

% -- Create grid 

[X, Y] = meshgrid (x3 , x2 ) ; 


% Load eps.dat file 
load eps.dat; 


f=l ; 

for i=l:m; 
for j=l:n; 

epsll (i, j ) 
eps22 (i, j ) 
eps33 (i, j ) 
eps23 (i, j ) 
epsl3 (i, j ) 
epsl2 (i, j ) 
f=f +1 ; 

end 


eps (f , 1) ; 
eps (f , 2) ; 
eps (f , 3) ; 
eps (f , 4 ) ; 
eps (f , 5) ; 
eps (f , 6) ; 
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end 


if nsub > 0 

num_plots = nell + ne22 + ne33 + ne23 + nel2 + nel3; 

if num_plots == 1 
n 1 = 2 ; 
n2 = 2; 

elseif num_plots == 2 
n 1 = 2 ; 
n2 = 2; 

elseif num_plots == 3 
nl = 2; 
n2 = 2; 

elseif num_plots == 4 
nl = 2; 
n2 = 2; 

elseif num_plots == 5 
n 1 = 3 ; 
n2 = 3; 

elseif num_plots == 6 
nl = 3; 
n2 = 3; 

end 

end 


num = 0; 

if nsub > 0 
figure; 
hold on 

end 


% -- ell plot 


if nell > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, epsll ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [ell_low ell_high] ) ; 

end 

colorbar ; 

title ( ' \bf\ font size { 24 } \epsilon\ font size { 14 }_1_1 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 1 , ' rotation ' , 0 ) ; 

end 
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% -- e22 plot 


if ne22 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, eps22 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [e22_low e22_high] ) ; 

end 

colorbar; 

title ( ' \bf \ fonts ize { 24 } \epsilon\ fonts ize { 14 }_2_2 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bf x_2 ' , ' rotation ' , 0 ) ; 


end 


% -- ep33 plot 


if ne33 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figured- 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Yd eps33 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [e33_low e33_high] ) ; 

end 

colorbar; 

title ( ' \bf\ fonts ize { 24 } \epsilon\ fonts ize { 14 }_3_3 ' ) ; 

xlabel ( '\bfx_3' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' ,0) ; 


end 


% -- e23 plot 


if ne23 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
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hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, eps23 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [e23_low e23_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \ gamma \ font size { 14 }_2_3 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ',0) ; 


end 


% -- e!3 plot 


if nel3 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, epsl3 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [el3_low el3_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \ gamma \ font size { 14 }_1_3 ' ) ; 

xlabel ( '\bfx_3' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- el2 plot 


if nel2 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, epsl2 ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [el2_low el2_high] ) ; 

end 
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colorbar ; 

title ( ' \bf \ fonts ize { 24 } \ gamma \ font size { 14 }_1_2 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 

clear; 

break 


% -- End of file 


epsp .m 


epsp.m MATLAB source file 


% Distributed with MAC/GMC 4.0 

% This program will generate fringe plots with the MATLAB software 
% from the data written to ASCII files by MAC/GMC 4.0 using the 
% * MAT LAB option 

% This program reads the file epsp.dat, which contains 8 columns: 

% epll ep22 ep33 ep23 epl3 epl2 ep_eq matnum 

% ** NOTE: inelastic shear strains are engineering shears ** 

% NOTE: In the case of incremental plasticity constit model, 

% ep_eq is the integrated equivalent plastic strain increment 

% where dep_eq = sqrt (2*dei jp*deijp/3) 

% For other constit models, ep_eq = sqrt (2*ei jp*ei jp/3) 

% see Mendelson, A. (1968) Plasticity: Theory and Application. 

% While the ASCII output files from MAC/GMC 4.0 contain data for 
% all specified times (and integration points for laminate) , the 
% epsp.dat file should contain just the inelastic strain field data 
% for the one particular time (and integration point) at which the 
% fringe plots will be made. 

% In addition, this program reads the files x2.dat and x3.dat, 

% which are also generated by MAC/GMC 4.0 to define the subcell 
% grid. Note that this grid will differ between GMC and HFGMC . 


% -- Clear workspace 
clear; 


USER OPTIONS 


% -- Select inelastic strain components to plot (1 - plot, 0 - don't plot) 

nepll = 1; 

nep22 = 1; 

nep33 = 1; 

nep23 = 1; 

nepl3 = 1; 

nepl2 = 1; 

nepeq = 1; 

nmatn = 1; 


% -- Select subplot (all on one Fig.) 
% nsub = 1 --> subplot 

% nsub = 0 --> individual Figs, 

nsub = 0; 


or individual Figs . 
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% -- Set scale_spec = 1 to specify scale high and low values for each 
% component - set values below 
scale_spec = 1; 

if scale_spec > 0 

ep 1 1_1 ow = -0.005; 
epll_high = 0.015; 

ep22_low = -0.005; 
ep22_high = 0.015; 

ep33_low = -0.005; 
ep33_high = 0.015; 

ep23_low = -0.005; 
ep23_high = 0.015; 

ep 1 3_1 o w = -0.005; 
epl3_high = 0.015; 

epl2_low = -0.005; 
epl2_high = 0.015; 

epeq_low = 0.0; 
epeq_high = 0.03; 

matn_low = 0; 
matn_high = 4; 


end 


% -- Note: To draw subcell boundaries, change "shading interp" 
% to "shading faceted" (GMC only) 


END USER OPTIONS 


% -- Load xl.dat and x2.dat files and determine their sizes 
load x2.dat; 
m = length (x2 ) ; 

load x3.dat; 
n = length (x3 ) ; 

% -- Create grid 

[X, Y] = meshgrid (x3 , x2 ) ; 


% Load epsp.dat file 
load epsp.dat; 


f=l ; 

for i=l:m; 
for j=l:n; 

epsllp (i, j ) 
eps22p (i, j ) 
eps33p (i, j ) 
eps23p (i, j ) 
epsl3p (i, j ) 
epsl2p (i, j ) 
epspq ( i , j ) 
matnum (i, j ) 
f=f+l ; 

end 

end 


epsp (f , 1) ; 
epsp (f , 2) ; 
epsp (f , 3) ; 
epsp (f , 4 ) ; 
epsp (f , 5) ; 
epsp (f , 6) ; 
epsp (f , 7) ; 

: epsp (f , 8) ; 
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if nsub > 0 

num_plots = nepll + nep22 + nep33 + nep23 + nep!2 + nep!3 + nepeq + nmatn; 


if num 

plots == 1 

nl 

"= 2; 

n2 

= 2; 

elseif 

num plots — 2 

nl 

= 2; 

n2 

= 2; 

elseif 

num plots == 3 

nl 

= 2; 

n2 

= 2; 

elseif 

num plots == 4 

nl 

= 2; 

n2 

= 2; 

elseif 

num plots == 5 

nl 

= 3; 

n2 

= 3; 

elseif 

num plots == 6 

nl 

= 3; 

n2 

= 3; 

elseif 

num plots == 7 

nl 

= 3; 

n2 

= 3; 

elseif 

num plots — 8 

nl 

= 3; 

n2 

= 3; 


end 


end 


num = 0; 

if nsub > 0 
figure; 
hold on; 

end 


% -- epll plot 


if nepll > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, epsllp ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [epll_low epll_high] ) ; 

end 

colorbar ; 

title ( ' \bf \ fonts ize { 24 } \epsilon\ fonts ize { 14 } A p_l_l ' ) ; 

xlabel ( ’ \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 
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% -- ep22 plot 


if nep22 > 0 

num = num + 1 ; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, eps22p ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [ep22_low ep22_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \epsilon\ font size { 14 } A p_2_2 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- ep33 plot 


if nep33 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, eps33p ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [ep33_low ep33_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 2 4 } \epsilon\ font size { 14 } A p_3_3 ' ) ; 
xlabel (' \bfx_3' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- ep23 plot 


if nep23 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figured- 
hold on; 

end 
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colormap (jet) ; 

pcolor (X, Y, eps23p ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [ep23_low ep23_high] ) ; 

end 

colorbar; 

title ( ' \bf \ fonts ize { 24 } \gamma\ font size { 14 } A p_2_3 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- epl3 plot 


if nepl3 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, epsl3p) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [epl3_low epl3_high] ) ; 

end 

colorbar; 

title ( ' \bf\ fonts ize { 24 } \gamma\ font size { 14 } A p_l_3 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bf x_2 ' , ' rotation ' , 0 ) ; 


end 


% -- epl2 plot 


if nepl2 > 0 

num = num + 1; 
if nsub > 0 

subplot (nl,n2, num) ; 

else 

figured- 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Yd epsl2p ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [epl2_low epl2_high] ) ; 

end 

colorbar ; 
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title ( ' \bf \ fonts ize { 24 } \gamma\ font size { 14 } A p_l_2 ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bf x_2 ' , ' rotation ' , 0 ) ; 


end 


% -- ep_eq plot 


if nepeq > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, epspq ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [epeq_low epeq_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size { 24 } \epsilon\ fonts ize { 14 } A p_e_q ' ) ; 

xlabel ( ' \bfx_3 ' ) ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 


% -- matnum plot 


if nmatn > 0 

num = num + 1; 
if nsub > 0 

subplot (nl , n2 , num) ; 

else 

figure; 
hold on; 

end 

colormap (jet) ; 

pcolor (X, Y, matnum ) , shading interp; 

axis image; 

if scale_spec > 0 

caxis ( [matn_low matn_high] ) ; 

end 

colorbar; 

title ( ' \bf\ font size {20 {materials ' ) ; 

xlabel ('\bfx_3') ; 

ylabel ( ' \bfx_2 ' , ' rotation ' , 0 ) ; 


end 

clear; 

break 


% -- End of file 
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Keywords 

♦ALLOWABLES, 80 
♦CHECK, 13 
♦CONDUCTIVITY, 14 
♦CONSTITUENTS, 21 
♦CURTIN, 98 
♦DAMAGE, 88 
♦DEBOND, 94 
♦ELECTROMAG, 16 
♦END, 113 

♦FAILU RE_CELL , 86 
♦FAILU RE_SUBCELL , 83 
♦LAMINATE, 60 
♦MATLAB, 110 
♦MDBPATH, 19 
♦MECH, 66 
♦PATRAN, 108 
♦PRINT, 101 
♦RUC, 48 
♦SOLVER, 72 
♦SURF, 76 
♦THERM , 70 
♦XYPLOT, 102 

A 

A—, 49, 89, 90 

ACTION=, 83, 84, 85, 86, 87, 98, 99 

Adams-Bashforth, 74 

ALF1=, 31, 33 

ALF2=, 31, 33 

ALF3=, 31, 33 

ALF4=, 31, 33 

ALF5=, 31, 33 

ALF6=, 31, 33 

all subcells, 80, 91 

allowables, 1, 5, 79, 80, 81, 82, 83 

ALPA=, 26, 29, 30, 34, 35, 41, 43, 129 

ALPT=, 26, 29, 30, 34, 35, 41, 43, 129 

aluminum, 39, 129 

ANG=, 60,61,62, 63,89, 90 

ANGINC=, 76 

anisotropic, 3, 5, 21, 30, 32 

arbitrary loading, 17 

ARCHID=, 11, 48, 49, 50, 51, 56, 58, 60, 61, 62, 63 
architecture, 1, 2, 45, 48, 49, 50, 51, 57, 60, 61, 62, 100, 
101 

architecture library, 48, 49, 51, 61 
array dimensions, 1 3 
AS=, 41, 43 

ASCII files, 3, 9, 10, 100, 101, 130, 136, 141 
ASP=, 49, 50, 61 
ASP1=, 49, 61 


ASP2=, 49, 61 
aspect ratio, 50, 51 
average material, 80 

B 

B=, 49 

BDN=, 94, 96 
BDS=, 94, 96 
BE=, 89, 90 
BLOCK=, 89, 90, 92 
BN=, 89, 90, 94, 96 

Bodner-Partom, 3, 21, 24, 28, 39, 118, 127, 129 

bond strength, 94 

BP=, 89, 90 

BS=, 94, 96 

build, 10 

c 

C=, 49 

Cll=, 31, 32 
C12=, 31, 32 
C13=, 31,32 
C14=, 31, 32 
C15=, 31,32 
C16=, 31,32 
C22=, 31,32 
C23=, 31,32 
C24=, 31,32 
C25=, 31,32 
C26=, 31,32 
C33=, 31,32 
C34=, 31,32 
C35=, 31,32 
C36=, 31,32 
C44=, 31,32 
C45=, 31,33 
C46=, 31,33 
C55=, 31,33 
C56=, 31,33 
C66=, 31, 33 

characteristic length, 98, 99 
characterization, 27, 28, 99 
check input file, 12 

CMOD=, 11, 16, 22, 24, 25, 26, 28, 29, 30, 31, 32, 33, 35, 
36, 37, 39,40,41,43, 129 
coefficient of thermal expansion, 27, 32, 34, 42 
Coffm-Manson Law, 88 
command prompt, 7, 10 
Compaq Visual Fortran, 9, 10 
compile, 9, 10, 36, 38 
COMPR=, 80, 81, 83, 84, 86, 87 
concentration equations, 79 
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conductivity, 14, 27, 34, 36 
consistency, 40 

constant compliant interface (CCI) model, 94, 96 
constitutive model, 2, 3, 5, 9, 16, 18, 21, 22, 23, 24, 25, 26, 
28, 30, 33, 34, 36, 37, 38, 40, 44, 73, 1 16, 121, 124, 127, 
129 

controlling component, 77 
controlling subcell, 88, 89 
COORD=, 102, 103, 104 
coordinate system, 1 03 
copper, 39, 88, 129 
creep time, 105, 106, 107, 109 
cross shaped fiber, 50 
Curtin model, 1, 4, 6, 79, 98, 99 
cyclic, 79, 88, 91 

D 

D=, 22, 23, 26, 28, 29, 36, 41, 43, 44, 49, 56, 58, 59, 62, 63, 
98,99 

damage, 1, 2, 6, 74, 79, 85, 87, 88, 89, 90, 91, 92, 93, 101, 
116, 123, 124, 127, 128 
damage increment, 88, 89, 91 
damage parameter, 88, 89, 91 
data output, 1, 100, 108, 110 
DBCH=, 94, 96 
debonded face, 95 

debonding, 1, 2, 4, 6, 79, 91, 94, 96, 97, 107, 117, 124 
debonding parameter, 94, 96, 99 
debug, 10 

DELAY=, 94, 96, 98, 99 
design allowables, 80 
DINC=, 89, 90 

directory, 7, 8, 9, 10, 19, 104, 111 
directory path, 8, 9, 19 
dissipation rate, 76 
DOS, 7, 10 
DR=, 49,61,76 

E 

El=, 33, 35, 126, 129 
El 1=, 80, 81 
E12=, 80, 81 
E13=, 80, 81 
E2=, 33, 35, 126, 129 
E22=, 80, 81 
E23=, 80, 81 
E33=, 80, 81 

EA=, 26, 29,30,41,43, 129 
ECU—, 80, 81 
EC22=, 80 
EC33=, 80 

effective constitutive equation, 64 
effective properties, 11, 45, 101 
EL=, 11, 26, 29, 30, 32, 33, 35, 39, 41, 43, 129 
electric displacement, 66, 69, 106, 107 
electric field, 66, 68, 69, 105, 106, 107 
electromagnetic, 5, 12, 16, 17, 20, 41, 42, 43, 44, 45, 60, 
61,62, 63,64, 65,66, 69 


electromagnetic specifier, 16, 17, 20, 44, 61, 63 
electromagnetically active, 44, 61, 63, 69 
EM=, 16, 17, 20, 41, 43, 44, 60, 61, 62, 63 
embedded, 45 
ending time, 71, 73 

engineering shear strain, 64, 68, 82, 85, 87, 105, 106, 107 
EPS=, 76 

equivalent plastic strain, 76, 141 

ERR=, 72, 73, 74, 75 

errors, 73, 74 

ES=, 41, 43 

ESI 1=, 41, 43 

ES12=, 41, 43 

ES26=, 41,43 

ET=, 26, 29,30,41,43, 129 

ETFL=, 89, 90 

ETM=, 89, 90 

ETU=, 89, 90 

evolution equations, 96 

evolving compliant interface (ECI) model, 94 

Example Problem Manual, 2, 3, 6, 7 

executable file, 7, 9, 10, 36, 38 

execution, 5, 8, 10, 12, 13, 16, 69, 75, 77, 84, 86, 88, 91, 

98, 101, 108, 111 
execution speed, 1 0 

external material database, 1, 5, 9, 18, 19, 20, 22, 27, 32, 
34, 37, 39, 40, 42, 115, 129 
EXTMAT=, 39 

F 

F=, 11,48,49, 50,51,61,63 
FACE=, 94, 96 

failure, 2, 5, 6, 79, 81, 83, 84, 85, 86, 87, 88, 91, 94, 98, 99, 
101, 123, 124, 127 
failure criteria, 5, 79, 83, 86, 123 
failure parameters, 83, 84, 85, 86, 87 
fatigue, 6, 74, 79, 85, 87, 88, 89, 90, 91, 92, 93, 101, 123, 
124, 128 

fatigue damage exponent, 90 
fatigue limit, 90 
FeCrAlY, 24 

fiber breakage, 1, 4, 6, 79, 98, 99 

fiber diameter, 50, 99 

fiber fracture, 6, 98 

fiber material, 5 1 

fiber radius, 50 

fiber strength, 6, 79, 98 

fiber volume fraction, 51, 55 

field variables, 38, 72, 73, 75, 123, 124, 125 

first subcell, 80, 124 

flag-type keywords, 12 

force resultant, 66, 68, 76, 77, 105 

FORTRAN, 9, 10 

forward Euler, 65, 73, 77 

Freed-Walker viscoplastic model, 3, 5, 21, 25, 28 

FREQ=, 102, 103 

fringe plots, 1,4, 6, 100, 108, 109, 111, 115, 130, 136, 141 
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G 

GA=, 26, 29, 30,41,43, 129 
Gauss quadrature, 45, 47 
Getting Started, 1,6,7 
global equilibrium, 72, 73 
graphite, 11 

GVIPS, 3, 5, 21, 24, 25, 28, 73, 74, 101 

H 

H=, 49, 56, 58, 62, 63 

HFGMC, 2, 3, 4, 5, 6, 45, 46, 50, 51, 57, 59, 60, 61, 62, 63, 
73,94, 130, 131, 136, 141 

I 

I=, 48, 49, 50,51 
IA=, 102, 104 
IB=, 102, 103 
IG=, 102, 103 

incremental plasticity model, 3, 5, 36, 73, 74, 101, 141 

indices, 46 

ineffective length, 98 

inelastic material property, 27, 28, 34 

inelastic power, 76 

inelastic strain rate, 76 

inelasticity, 2 

initial time step, 75 

input file, 1,3,4, 6, 7, 8, 9, 10, 11, 12, 13, 18, 19, 20,21, 
26, 27, 32, 34, 37, 39, 40, 42, 44, 45, 48, 60, 61, 87, 104, 
112, 113, 127, 128 

integration, 1, 45, 64, 65, 72, 73, 74, 77, 130, 131, 136, 141 
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K=, 11,26, 29, 33,35, 36,39 
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KS1 1=, 41, 43 
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MATDB=, 11, 26, 29, 30, 31, 32, 33, 35, 37, 38, 39, 41, 43 

material assignment, 57, 62 

material number, 22, 27, 31, 34, 37, 38, 39, 42, 50, 81, 83, 
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matrix material, 5 1 
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particulate, 5 
PATH=, 19 
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PDF, 7 
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PELS1=, 42, 43 
PELS2=, 42, 43 
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piezoelectric, 42 
piezomagnetic, 42 

plane of transverse isotropy, 22, 27, 42 
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96, 98, 99, 103, 104 
PMGS=, 41, 43 
PMGS1=, 42, 43 
PMGS2=, 42, 43 
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post-processing, 1, 6, 7, 100, 108, 110 
preload, 65, 77 
print level, 101 
probing, 65, 76, 77, 78 
project, 9, 10 
pyroelectric, 42, 43, 105 
pyromagnetic, 42, 43, 105 
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QS=, 41, 43 
QS1 1=, 41, 43 
QS12-, 41, 43 
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R=, 48, 49, 50,51,61 
RA=, 49 
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redistribution, 88, 89, 94 
references, 7 
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reversal tolerance, 96 
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SCS-6, 24, 39, 129 
self-adaptive, 65, 72, 74 
semi-logarithmic, 88, 93 
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SFL=, 89, 90 

shape memory alloy, 3, 5, 21 
shear modulus, 27, 42 
shear strength, 90 
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sliding resistance, 98, 99 
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step size, 72, 73, 74, 75, 76, 77, 123 
stiffness matrix, 36, 64, 101, 121, 122, 123 
stiffness reduction, 6, 79, 88, 91 


STP=, 72, 76, 108 
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strength reduction, 4, 6, 79, 88, 90, 93 

stress allowables, 8 1 

stress-strain curve, 5, 27, 28, 99, 100 

SU=, 89, 90 

Sum, 89, 90 

SU2=, 89, 90 
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surf_lst . dat, 78 
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surf_global . dat, 78 

surf_local . dat, 78 
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T300, 24 
T50, 24 
TAU0=, 98, 99 

TEEN, 26, 29, 30, 31, 32, 33, 35, 41, 43, 129 
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thermal loading, 1, 64, 70, 73, 74, 75, 91 
THK=, 60, 61, 62, 63 
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TI=, 66, 67, 68, 70,71,72 
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Ti-6-4, 24 
time delay, 97, 99 
time integration, 65, 72, 73, 77 
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TIMES=, 110 
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144, 145, 146 
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ultimate strength, 90, 93 
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usrmat, 1, 3, 9, 21, 36, 38, 116, 141 
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Vl=, 26, 29,30, 34,35, 129 
V2=, 26, 29,30, 34,35, 129 
VF=, 11,48,49, 50,61,63 
VF1=, 49, 51 

VI=, 26, 29,33,35,36, 39, 129 
volume fraction, 50, 51, 55 
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Weibull modulus, 98 
Windows, 7, 9 
workspace, 130, 136, 141 
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X=, 102, 103 
XI 1= 83,84, 86, 87 
X12=, 83,84, 86, 87 
X13=, 83, 84, 86, 87 
X22=, 83, 84, 86, 87 
X23=, 83, 84, 86, 87 
X33=, 83, 84, 86, 87 
XA=, 48 
XC11=, 83, 86 
XC22=, 83, 84, 86, 87 
XC33=, 83, 84, 86, 87 
XML=, 89, 90 

x-y plot, 4, 100, 103, 105, 106, 107, 108 
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Y=, 102, 103 
yield criteria, 77 

yield surface, 2, 5, 27, 28, 64, 65, 73, 76, 77, 78, 101 
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zero stiffness, 88, 91 
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